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Section  1 


NTROOUCnON 


Delay  lines  employing  microwave  frequency  acoustic  waves  are  being 
used  in  radar  signal  processing  applications.  These  devices,  however,  have 
generally  been  found  practical  only  for  delay  times  up  to  a few  tens  of  micro- 
seconds. The  principal  factor  limiting  the  development  of  long  memory  time 
(milliseconds)  applications  has  been  the  attenuation  of  the  acoustic  waves  in 
the  crystals  used  as  the  delay  line  media.  In  general  this  attenuation  is  sub- 
stantially reduced  at  low  temperatures,  and  thus  longer  delay  time  applications 
may  become  feasible  if  the  additional  cryogenic  equipment  can  be  justified. 

As  advances  in  cryogenic  equipment  occur,  tradeoff  calculations  for  cryogenic 
delay  lines  in  practical  systems  will  have  to  be  updated.  Ambient  temperature 
operation  is  always  preferred,  however,  and  is  imperative  in  most  cases.  It 
has  been  the  objective,  therefore,  of  a continuing  program  at  the  General 
Electric  Research  and  Development  Center,  supported  by  U.  S.  Air  Force 
contracts,  to  develop  crystals  with  low  microwave  acoustic  attenuation  at 
ambient  temperature. 

Throughout  the  prior  work  at  General  Electric  emphasis  has  been  given 
to  a systematic  approach  to  the  development  of  useful  crystals.  Very  briefly, 
the  theory  of  microwave  frequency  sound  attenuation  was  put  in  a form  which 
permitted  correlation  of  attenuation  data  near  room  temperature  for  various 
crystals,  making  it  possible  to  define  a figure  of  merit  for  candidate  crystals. 
The  value  of  tnis  figure  of  merit  could  then  be  calculated  from  known  properties 
or  from  the  results  of  relatively  simple  thermal  measurements.  The  figure 
of  merit  is  applicable  in  the  high  temperature  regime,  including  room  tem- 
perature, where  the  attenuation  is  .uused  by  the  scattering  of  sound  waves  by 
thermal  phonons.  Theory  also  exit-.  3 for  the  temperature  region  in  which  a 
rapid  falloff  in  attenuation  occurs.  This  theory  can  be  put  in  a form  which 
will  provide  a figure  of  merit  for  materials  at  cryogenic  temperatures.  At 
very  low  temperatures,  observed  sound  attenuation  is  dominated  by  scattering 
due  to  crystal  imperfections  of  various  kinds.  In  highly  perfect  specimens  the 
attenuation  can  be  very  small  at  sufficiently  low  temperatures,  and  is  deter- 
mined mainly  by  diffraction  losses. 

Many  crystal  systems  have  been  investigated,  both  as  potential  practical 
materials  and  to  verify  the  figure  of  merit  theory.  The  lowest  known  room 
temperature  attenuation  coefficients  which  have  been  reliably  measured  are 
for  shear  waves  in  MgAl£>4  (spinel)  and  Y3  AlBOxi  (garnet), and  for  longitudinal 
waves  in  0- rhombohedral  boron  and  AlgOg  (sapphire).  At  present,  delay  lines 
are  constructed  principally  from  sapphire,  using  longitudinal  acoustic  waves, 
or  from  spinel,  using  shear  waves. 

The  need  for  broad  bandwidth  in  delay  lines  does  not  significantly  limit 
the  choice  of  crystal  media  because  sound  dispersion  is  negligible  in  the  fre- 
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qucncy  range  of  interest  for  almost  all  crystals.  Bandwidth  enters  most  im- 
portantly in  the  design  of  transducers  and  the  transducer  coupling  circuits. 

The  transducer  design,  electrode  materials,  and  the  coupling  circuits  also 
have  very  important  effects  on  the  performance  ol'  delay  lines  made  from  mate- 
rials such  as  spinel.  Shear  mode  transducers  require  improvements  in  the 
electrode  materials  used,  in  the  types  of  bonds  employed  with  bonded  crystal 
plates,  and  in  reproducible  fabrication  of  deposited  thin  film  transducers. 

In  Section  2 of  this  report,  "Design  of  Delay  Lines,  " the  practical  design 
is  discussed  thoroughly.  The  treatment  includes  those  factors  wMch  are 
under  the  designers  control:  attenuation  in  transducer,  electrodes,  and  propa- 
gating medium:  transducer  efficiency  and  bandwidth;  diffraction  of  the  acoustic 
wave:  and  beam  steering  due  to  elastic  anisotropy.  The  results  give  appropriate 
criteria  for  transducer  design,  show  the  effects  of  transducer  design  upon 
performance,  and  show  the  relative  performance  capabilities  for  sapphire, 
yttrium  aluminum  garnet,  and  spinel.  The  calculations  give  a reasonable  up- 
per limit  to  performance,  since  such  factors  as  crystalline  imperfections  and 
the  limits  of  oolishing  and  fabrication  methods  cannot  be  predicted  theoretically. 

The  implications  of  matched  filtering  by  means  of  silicon  clocked  analog 
delay  devices  (bucket  brigade  or  charge  coupled  devices)  and  surface  acoustic 
wave  devices  are  examined  and  related  to  the  capabilities  of  Bilicon  clocked 
analog  delay  (CLAD)  devices,  passive  surface  acoustic  wave  (SAW)  delay 
lines,  and  bulk  acoustic  wave  delay  lines  used  as  memory  elements.  With 
this  information  as  a background,  the  effect  of  acoustic  impedance  mismatch 
between  delay  line  and  transducer  is  examined.  The  problem  is  treated  as  a 
multipath  problem  in  a form  applicable  to  arbitrary  coding  of  the  input  signal 
which  is  to  be  stored  in  the  delay  line.  Simple  criteria  are  developed  to 
describe  the  effects  of  acoustic  impedance  mismatch  on  both  the  zoded  signal 
and  the  matched  filter  output  for  that  signal. 

In  Section  3,  "Delay  Line  Fabrication,  " particular  attention  is  paid  to 
the  problem  of  attaching  piezoelectric  crystals  to  delay  line  crystals  by  means 
of  an  intermediate  metallic  bond.  The  properties  of  such  bonds  are  discussed 
and  those  metals  appropriate  to  spinel  and  lithium  niobate  are  indicat’d.  A 
scheme  for  fabricating  these  bonds  has  been  developed,  and  methods  and 
equipment  for  the  subsequent  thinning  of  the  transducer  crystals  are  discussed. 

In  Section  4,  "Growth  and  Characterization.  " the  whole  development  of 
YBm  crystals  for  delay  lines  is  discussed.  Particular  attention  is  paid  to 
crystal  purity,  stoichiometry,  and  perfection.  Stoichiometric  crystals  appear 
to  have  a composition  of  YBW,  while  the  congruent  melting  composition  is 
ybm.7.  Material  with  low  acoustic  loss  may  be  YBa4,  This  knowledge  com- 
bined with  acoustic  loss  measurements  on  actual  samples  Bhows  that,  at  pre- 
sent, these  crystals  are  not  suitable  for  low  loss  delay  lines. 

In  Section  5,  "iron  Doped  Spinel  Crystals.  " synthetic  crystals  of  iron 
doped  spinel  are  analyzed  for  their  suitability  for  delay  lines.  In  Section  6, 


"Sound  Attenuation  and  Diffraction  Loss,  " measurements  of  acoustic  attenuation 
at  low  temperature  and  low  frequency  are  compared  with  the  calculated  diffrac- 
tion loss.  Good  agreement  is  found. 
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Section  2 


DESIGN  OF  DELAY  UNES 


A.  BACKGROUND 

In  this  section  the  design  of  acoustic  delay  lines  for  long  reverberation 
times  is  considered.  Delay  time  is  limited  by  material  choice.  Calculations 
are  made  here  for  practical  delay  devices  in  order  to  establish  which  com- 
binations of  transducer,  electrode  material,  and  delay  medium  will  provide 
the  longest  delay  time  and  the  longest  delay  time- bandwidth  product.  The 
calculations  include  the  major  factors  which  affect  delay  time:  attenuation  in 
the  delay  medium,  attenuation  in  the  metal  electrodes,  conversion  of  acoustic 
energy  to  electrical  energy  by  the  transducer,  diffraction,  and  acoustic  beam 
steering  due  to  the  elastic  anisotropy  of  a crystal  medium.  The  calculations 
are  a performance  upper  limit  in  the  sense  that  a number  of  small  and  unpre- 
dictable effects  (nonparallelism  of  end  surfaces,  surface  damage  from  polishing, 
internal  strain  due  to  crystal  growth,  local  imperfections)  are  not  included 
in  the  calculations.  The  results  for  sapphire,  however,  correspond  reason- 
ably well  with  experiment.  It  is  shown  that  other  materials,  particularly 
MgAljD4  spinel,  can  provide  substantial  improvement  in  performance  over 
AlfOj,. 


Attenuation  in  the  electrodes  and  the  delay  medium  increases  quadratically 
with  frequency,  and  diffraction  losses  vary  inversely  with  frequency.  At  low 
frequencies  delay  time  is  limited  by  diffraction;  at  high  frequencies,  by  at- 
tenuation. Delay  time  and  delay-time  bandwidth  product  each  go  through  a 
maximum  at  intermediate  frequencies.  For  delay  line  approximately  3/8  of 
an  inch  in  diameter  these  maxima  occur  between  200  MHz  and  1 GHz.  The 
bandwidth  of  the  devices  is  approximately  one-half  the  operating  frequency, 
or  from  100  MHz  to  500  MHz  in  the  0.  2 to  1.0  GHz  range.  Since  these  band- 
widths  are  adequate  for  most  applications  and  signals  can  be  readily  hetero- 
dyned to  this  frequency  range,  calculations  have  been  carried  out  in  the  fre- 
quency interval  of  0.  2 to  1.  0 GHz.  Delay  times  can  be  increased  by  considering 
delay  lines  of  larger  diameter,  but  the  improvement  comes  at  lower  frequencies, 
where  diffraction  limits  the  delay  lime.  At  the  delay  time  maximum,  the  op- 
erating frequency  is  proportional  to  D”l/a  and  the  maximum  delay  time  is 
proportional  to  D“/#  where  D is  the  transducer  diameter. 


To  calculate  a delay  time,  assumptions  must  be  made  about  the  initial 
signal  level  and  receiver  sensitivity.  These  assumptions  do  not  affect  the 
relative  performance  of  delay  lines  fabricated  from  different  combinations  of 
materials.  The  relative  performance  shown  by  these  calculations  will  be  valid 
for  a variety  of  applications.  To  show  the  magnitude  of  delay  which  may  be 
expected,  a set  of  typical  system  parameters  --  drive  power,  receiver  sen- 
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sitivity,  intcrpulse  period  --  were  chosen.  The  total  delay  time  in  a specific 
application  will  differ  from  these  calculated  values  if  the  system  parameters 
are  different.  That  is,  this  calculation  gives  the  relative  performance  under 
a variety  of  system  parameters  and  several  delay  line  material  combinations, 
but  the  calculated  total  delay  times  are  typical  values  valid  for  a specific  set 
of  system  parameters.  A simple  system  will  be  assumed,  with  an  input  of 
+40  dBm  (10  watts  peak)  into  the  delay  line  and  a low -noise  amplifier  on  the 
output  with  a sensitivity  of  -70  dBm.  The  system  thus  has  a dynamic  range 
of  110  dB  with  a good  signal -to-noise  ratio  which  will  remain  constant  for 
all  cases  discussed.  Three  delay  line  configurations  will  be  analyzed:  the 
double-end  (DE)  delay  line,  with  transducers  on  both  input  and  output;  the 
single-end  (SE)  delay  line,  which  requires  a circulator;  and  the  switched  im- 
pedance delay  line,  which  requires  an  electronic  switch  and  two  impedance 
transformers.  The  delay  line  materials  will  be  the  lowest  loss  materials 
known  that  are  readily  available  in  large  size  — longitudinal  propagation  in 
sapphire  (Al3Oa)  and  shear  pr.  pagation  in  YAG  (Y^AlgO^  and  spinel  (MgAlg04). 

The  transducer  thickness  in  the  frequency  range  0.2  to  1 GHz  varies  from 
5 to  11  pm,  depending  on  the  transducer  material  usee  Either  bonded  plate 
transducers,  thinned  after  bonding,  or  deposited  transducers  can  be  UBed. 

The  low  loss  acoustic  modes  are  shear  waves  and  require  shear  transducers. 
Now  bonding  methods  are  required  for  the  use  of  bonded  plate  transducers  with 
low  acoustic  loss  metals  for  electrodes.  Deposited  shear  mode  transducers 
have  never  been  deposited  reproducibly.  Transducer  development  will  be  re- 
quired in  Implementing  long  delay  time  devices  with  either  type  of  transducer. 

The  properties  of  transducer  materials,  electrode  materials,  and  delay 
line  materials  needed  for  the  loss  calculations  are  listed  in  the  following  sub- 
section. Next,  the  expression  for  the  total  loss  in  a reverberating  delay  line 
is  derived,  as  well  as  the  basic  equations  necessary  to  calculate  the  acoustic 
conversion  loss,  the  electrode  loss,  the  material  loss,  and  the  acoustic  diffrac- 
tion loss.  A series  of  graphs  compares  the  total  delay  time  versus  insertion 
loss  (1L)  and  frequency  and  the  IL  versus  frequency  for  several  delay  line  con- 
figurations and  materials. 

B.i  TRANSDUCER.  ELECTRODE.  AND  DELAY  MffiiA  PROPERTIES 

The  material  properties  of  transducers,  electrode,  and  propagation  media 
needed  to  calculate  the  acoustic  loss  are  listed  in  Tables  1 through  3. 

The  data  in  Table  1 are  from  Reeder  and  Winslow  (Ref.  1).  All  these 
materials  can  be  cut  into  thin  transducer  plates.  The  longitudinal  mode 
ZnO  transducer  can  also  be  sputtered  with  reduced  coupling  coefficient 
kt  = 0.22  to  0.25. 

All  the  data  in  Table  2 were  experimentally  determined  by  Larson  and 
Winslow  (Ref.  2)  by  fitting  parameters  to  insertion  loss  versus  frequency  data 
for  delay  lines,  using  thin  films  of  these  metals  as  electrodes.  The  propagation 
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Table  1 


PROPERTIES  OP  TRANSDUCER  MATERIALS 


Material 

Cut 

Mode 

Coupling 

Coefficient 

k 

Dielectric 

(.'■instant 

% 

Velocity 
V,  (10*  cm/sec) 

Acoustic 
In>  pedance 
Z,  (I0*g/s  cm9! 

Quart*,  SIO, 

mm 

Shear 

k„  • 0.  14 

4.5 

3.80 

10.  1 

Cadmium  sulfide,  CdS 

■9 

Long 

k,  = 0. 1 5 

8.  5 

4.46 

215 

Zinc  oxide,  /nO 

mm 

Long 

k,  • 0.  28 

8.  8 

6.  33 

38.0 

Zinc  oxide,  ZnO 

Shear 

- 0.32 

8.  3 

2.  n 

15.  5 

Lithium  niohate,  I.iNMI, 

KS 

Shear 

k > 0.68 

44.0 

4.  80 

22.6 

Table  2 

PROPERTIES  OF  ELECTRODE  MATERIALS 


M:it*»rta  l 

Propagation 
[ >t t iM:in»n 

V4  4 1 O'"  i'tii  1 i»er) 

Acoustic  Impedance 
Z,  <10»  g/s  cm*) 

AUnnuai  um 
tllJ/otu  at  1 fill/ 

Indium 

■ 101 1 

i...ug. 

m 

Hi.  4 

HC.000 

tmiluM 

iuii 

Shear 

r».  4 

r.n.ooo 

Silver 

' s 1 i i 

long. 

41. 5 

250 

silver 

1 1 1 1 

Shear 

jfl  1 

Hi.  5 

Cicitfl 

■mi 

Long. 

65.  f. 

200 

fluid 

ini! 

Slice  r 

n.4 

1.000 

Table  3 

PROPERTIES  OF  DELAY  MEDIA 


Mali'l'inl 

Propa  nation 
fNrrt'tfun 

Mode 

Velocity 
V„  ( IUf  cm/ace) 

A'-iuis’h*  t<tip<»diinri* 
/. . (10*  IJ/.M  onr) 

Atlrmiut.on 
IdB/p*  nt  1 till/.) 

Sapphire  AigO, 

/. 

Long, 

1 1.  1 

14.  i 

o.  :• 

YAfl  YV\l,Ou 

riooi 

Shear 

5,  o:i 

22. !' 

U.  OH 

Spinel  MgAljOg 

fiooi 

Shear 

6.  54 

23.  4 

0.  04 

directions  listed  are  the  orientations  usually  obtained  with  vapor -deposited 
metal  films.  The  velocities  are  slightly  different  from  the  values  published 
for  bulk  materials.  The  attenuation  values  are  normalized  ;o  1 GHz  and  have 
a frequency  squared  dependence;  thus  the  attenuation  valuer  Rated  at  1 GHz 
will  be  reduced  by  a factor  of  25  at  a frequency  of  200  MHz,  The  shear  at- 
tenuation value h are  always  higher  than  the  longitudinal  attenuation  values. 
Indium  is  not  a good  bond  material  for  reverberating  delay  lines  because  of 
its  excessively  high  loss.  Delay  line  performance  is  calculated  for  devices 
made  with  gold  electrodes  and  for  delay  lines  in  which  the  attenuation  in  the 
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electrode  is  negligible.  These  calculations  will  show  the  desirability  of  re- 
placing gold  with  a lower  loss  metal  in  the  case  of  the  best  combination  of  trans- 
ducers and  delay  media. 

The  attenuation  values  listed  in  Table  3 are  taken  from  the  authors'  acous- 
tic attenuation  measurements  and  those  in  the  literature  at  1 GHz.  The  at- 
tenuation constants  for  these  materials  are  proportional  to  frequency  squared 
over  the  range  of  measurement  of  10  GHz  to  500  MHz.  From  500  MHz  to  200 
MHz  the  attenuation  values  are  extrapolated  with  the  quadratic  frequency  de- 
pendence. The  three  materials  considered  have  the  lowest  known  losses  for 
commercially  available  crystals  prepared  in  large  size.  The  attenuation  data 
for  spinel  were  obtained  on  small  flux  grown  samples.  Large,  high-quality 
Czochralski  boules  are  now  available,  but  it  has  not  yet  been  shown  that  this 
material  has  the  perfection  required  for  low  acoustic  losses. 

C.  TOTAL  ACOUSTIC  LOSS 

The  delay  time  of  a delay  line  is  determined  by  the  loss  characteristics 
of  the  line  and  the  dynamic  range  of  the  system  in  which  it  is  used. 

Dynamic  Range  = Transducer  Insertion  Loss  (loss  in  converting  input 

electrical  signal  into  acoustic  energy  and  back  into 
an  electrical  signal) 

4- 

Acoustic  Conversion  Loss  (acoustic  energy  converted 
to  an  electrical  signal  every  time  the  acoustic  pulse 
strikes  the  transducer) 

+ 

Electrode  Attenuation  (acoustic  energy  loss  in  the  metal 
electrodes) 


+ 

Propagation  Medium  Attenuation 

+ 

Diffraction  Loss  (diffraction  causes  the  acoustic  pulse 
to  spread  in  size  as  it  propagates;  thus  a decreasing 
fraction  of  it  strikes  the  transducer) 

If  the  transducer  is  made  100  percent  efficient  in  converting  the  in- 
put energy  into  acoustic  energy,  insertion  loss  = 0;  if  it  is  not  electrically 
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switched  to  change  its  efficiency,  then  the  first  acoustic  "echo"  to  strike  the 
transducer  is  completely  converted  hack  to  an  electrical  signal  and  the  total 
delay  time  is  very  small.  If  the  transducer  is  made  very  inefficient  the  acous- 
tic conversion  loss  will  be  small  and  very  little  of  the  acoustic  energy  will  be 
converted  to  electrical  energy  by  the  transducer;  however,  very  little  acoustic 
energy  will  be  obtained  from  the  input  electrical  signal  and  total  delay  time  will 
suffer.  A broad  optimum  of  the  insertion  loss  exists,  so  that  a balance  is  struck 
between  converting  the  initial  electrical  signal  to  acoustic  energy  and  the  con- 
versio  of  acoustic  energy  when  the  acoustic  pulse  strikes  the  transducer.  In 
the  calculations  which  follow,  the  optimum  insertion  loss  is  found,  and  the  total 
delay  time  is  then  calculated  from  the  dynamic  range  of  the  system,  the  in- 
sertion loss  value,  and  the  attenuation  and  diffraction  constants  of  the  materials 
used  in  the  delay  line. 


The  origin  of  these  losses  will  become  apparent  as  the  DE  delay  line 
shown  in  Figure  1 is  analyzed. 
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Figure  1.  Output  Power  of  the  First  Reflection  of  a DE  Delay  Line 

Assume  an  electrical  input  power  to  the  transducer  P0.  The  transducer  con- 
verts a fraction  r of  P0  into  acoustic  power,  which  propagates  through  the 
delay  material.  At  the  second  transducer  a fraction  r of  P0r  is  converted  to 
electrical  power  and  (1-r)  of  P0r  is  reflected. 


The  electrical  power  output  at  the  first  reflection,  Pl#  is  P0r8  exp 
(-2«ui  -o-ML-tt^L),  where  2 is  the  acoustic  loss  of  two  passes  through 
the  electrodes,  <?ML  is  the  acoustic  propagation  loss  through  the  delay  mate- 
rial, and  a „L  is  the  acoustic  diffraction  loss.  The  next  series  of  reflections 
is  shown  in  Figure  2. 
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Figure  2.  Output  Power  of  the  Second  and  Third  Reflections 
of  a DE  Delay  Line 
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The  general  expression  for  the  Nth  reflection  can  be  written  and  the  terms 
identified} 

P0r2(l-r)N“l  exp(-2Na,  l -No*  L-Na0  L)  = PH  <N=  1,  2.  3. . . ) (1) 

The  r*term  is  the  transducer  insertion  loss  IL  * 10  loglfira(dB){  the  (l-r)*“l 
term  is  the  acoustic  conversion  loss  due  to  the  transducer,  rc;  and  the  three 
exponential  terms  are  the  electrode  loss  IJ,  the  material  loss  n*»  and  the  dif- 
fraction loss  r0.  The  calculation  procedure  is  now  straightforward;  assumed 
values  for  IL  will  be  subtracted  from  the  constant  dynamic  range  (110  dB), 
and  the  result  will  be  called  the  available  power.  All  the  loss  terms  will  be 
calculated  on  a decibel  per  reflection  basis  and  summed.  This  number  will  be 
divided  into  the  available  power  to  calculate  the  total  number  of  reflections 
and  then  multiplied  by  the  travel  time  per  reflection  to  give  the  total  delay 
time.  Since  loss  will  be  calculated  for  a large  number  of  reflections  (N  = 

100  to  1000),  the  (1-r)  Iosb  will  be  included  in  the  first  reflection  and  the  cal- 
culation will  be  slightly  on  the  conservative  side  because  (l-r)*iB  being  used 
rather  than  the  (l-r)""1  in  Equation  1.  The  total  acoustic  loss  in  a SE  delay 
line  is  similar  to  Equation  1 with  one  important  difference  (see  Figure  3). 
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Figure  3.  Output  Power  of  the  First  and  Second  Reflections 
of  a SE  Delay  Line 
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(2) 


The  general  expression  for  the  Nth  reflection  of  a SE  delay  line  1st 
Po^U-r)**-1  exp(-2Nas  <-2Na*  L-2NaB  L)  - PM  (N*  1.  2,  3. . . ) 

The  SE  and  DE  cases  must  be  compared  at  the  same  delay  time;  N detected 
echos  and  2N,  respectively.  The  material  loss  and  the  diffraction  loss  are 
the  same.  However,  since  the  number  of  reflections  is  less  than  in  the  DE 
case  the  acoustic  conversion  loss  and  the  electrode  loss  are  reduced  in  the 
SE  line.  Thus,  as  will  be  seen  later,  the  number  of  reflections  for  the  SE 
case  is  always  greater  than  one -half  the  reflections  for  the  DE  case,  and  hence 
the  total  delay  time  is  greater  for  the  SE  delay  line.  Notice  the  way  in  which 
the  reflections  nre  counted:  a reflection  is  counted  every  time  it  passes  through 
a transducer.  This  must  be  done  in  the  DE  case  because  the  even-numbered 
transducers  have  acoustic  conversion  loss  and  electrode  loss  even  though  only 
the  odd-numbered  reflections  are  read  out.  All  reflections  are  read  out  in  the 
SE  case. 

In  the  electrically  switched  SE  delay  line  configuration,  the  driving  re- 
sistance is  switched  from  a low  value  with  an  efficient  transducer  conversion 
fraction  of  r1  to  a high  value  with  an  inefficient  transducer  conversion  fraction 
of  ra  for  readout.  The  general  expression  for  the  Nth reflection  is: 

Porlra(l-ra)N-1  exp<-2NaB  l -2Nol«  L-2Na„  L)  • H,  (N=  1,  2,  3.  . , ) (3) 

The  advantage  of  the  switched  configuration  is  that  the  input  electrical  pulse 
is  efficiently  converted  to  acoustic  power  but  the  transducer  is  then  switched 
so  that  it  abstracts  little  power  from  the  acoustic  pulse  each  time  it  Btrlkes 
the  transducer.  Thus  there  is  more  available  power  to  compensate  for  the 
acoustic  conversion,  electrode,  material,  and  diffraotion  losses. 

D.  BASIC  LOSS  EQUATIONS 

The  basic  equations  necessary  to  calculate  the  insertion  loss,  acoustic 
conversion  loss,  electrode  loss,  material  loss,  and  diffraction  loss  will  now 
be  derived. 

Insertion  Loss 


In  pulse-echo  measurements  transducer  IL  is  defined  as  the  ratio  of 
output  electric  power  (P *)  to  the  input  electrical  power  (P0)  in  the  absence  of 
attenuation  and  diffraction  losses  (see  Figure  1).  Thus  IL  is: 

p 

IL  * 10  log — I « 10  log  r8 
Po 

with  corrections  for  electrode  loss,  material  loss,  and  diffraction  of  the 
first  acoustic  echo.  The  IL  is  then 

IL  = 20  log  r decibels  (4) 
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Equation  4 will  be  used  to  calculate  r for  assumed  values  of  IL 
Acoustic  Conversion  Loss 


The  method  of  DeKlerk  (Ref*  3)  will  be  followed  in  calculating  the  acous- 
tic conversion  loss*  which  depends  on  the  electromechanical  coupling  coeffi- 
cient of  the  transducer.  This  acoustic  conversion  loss  is  the  ratio  of  the  re- 
flected power  to  the  incident  power  either  on  the  electrical  side  of  the  trans- 
ducer or  the  acoustic  side.  Referring  to  Figure  1,  this  ratio  is  Pe(l-r)/P0 
for  the  electrical  side  or  P0r(l-r)/P0r  for  the  acoustic  side.  Since  this  loss 
occurs  at  every  reflection,  the  acoustic  conversion  loss  per  reflection  is 

rt  s 10  log(l-r)  dB/ reflection  (5) 


Acoustic  Electrode  Loss 


The  electrode  losses  are  calculated  from  the  measured  attenuation  con- 
stants at  1 GHz  with  a frequency -squared  frequency  dependence  (see  values  In 
Table  2).  The  electrode  attenuation  constant  is: 

«• " A (nnnr)  dB/cm 

where 

A ■ 200  dB/cm  for  longitudinal  propagation  through  a gold  electrode 
■ 1000  dB/  cm  for  shear  propagation  through  a gold  electrode 
f ■ frequency  (MHz) 

Since  the  electrode  losses  are  high,  particularly  for  shear  waves,  only  1-jum- 
thick  gold  electrodes  will  be  considered.  There  are  two  passes  through  the 
electrode  per  reflection;  hence  the  acoustic  electrode  loss  per  reflection  is 

T,  * 2 X ICT*  * A dB/ reflection  (6) 

Acoustic  Material  Loss 


The  material  loss  is  calculated  in  a manner  similar  to  that  used  for  the 
electrode  loss.  Thus  the  material  loss  per  reflection  is 

T*B  (k^o)*  dB/reflectlon  (7) 

where 

r <=  travel  time  per  reflection  (L/Vd  for  DE  delay  lines; 

2L/V0  for  SE  delay  lines) 

B » attenuation  constant  (dB/usec)  from  Table  3 
f ° frequency  (MHz) 
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Acoustic  Diffraction  Loss 


Acoustic  diffraction  loss  occurs  because  the  cross-sectional  area  of  the 
acoustic  beam  does  not  remain  collimated  for  long  travel  distances  but  spreads 
out,  becoming  larger  than  the  receiving  transducer,  and  this  represents  an 
acoustic  energy  loss.  This  loss  is  difficult  to  calculate.  Until  recently  authors 
(Ref.  4)  considered  only  diffraction  from  circular -piston  pressure  sources 
radiating  into  fluid  or  elastic  isotropic  half-spaces.  Papadakis  (Ref.  5)  ex- 
tended this  work  to  include  longitudinal  wave  propagation  in  anisotropic  mate- 
rials. To  the  best  of  the  authors’  knowledge,  no  one  has  yet  calculated  diffrac- 
tion of  shear  waves  in  anisotropic  materials.  Papadakis'  results  will  be  used 
here  to  estimate  the  diffraction  loss  for  shear  waves. 

The  calculation  of  diffraction  for  longitudinal  wave  in  anisotropic  materials 
is  done  in  two  steps!  1)  finding  the  pressure  and  phase  profiles  at  several  dis- 
tances from  the  input  transducer,  and  2)  finding  the  response  of  the  receiver 
to  these  profiles.  Papadakis  Bhows  that  the  pressure  at  any  point  from  the  in- 
put transducer  is  given  by 

juipV  f exp  j {u>t-0o  |R|[l+b(l-2b)  9a]}  d o 

P ‘ "W, w 

where 

w ■ angular  frequency  of  the  acoustic  wave 
d ■ density  of  the  delay  medium 
VB  = velocity  along  the  direction  of  propagation 
0O  « propagation  vector  which  makes  an  angle  with 
the  transducer  normal  N 

R ■ a vector  lying  along  the  Poynting  vector  of  acoustic 
energy  flux 

o " area  of  the  transducer 

b a a coefficient  which  measures  the  elastic  anisotropy 
of  the  material 

The  second  Integration  is  performed  over  the  receiver  area  o',  using  the 
radial  coordinate  on  the  receiving  transducer  p'.  Because  the  receiver  is  a 
pressure-  and  phase- sensitive  detector  Its  signal  Is  proportional  to  the  max- 
imum over  one  cycle  of  the  integral 

H - f Ctp'lcostt-elp'ljdo'  I 

J0»  max 

where  C(p')  and  6(p')  are  the  amplitude  and  phase  of  the  pressure  profile 
given  above.  The  diffraction  loss  is  then  computed  from 

dB  - 20  log  ( |P.i  | / |1>,  |> 
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where  IV  is  the  pressure  at  some  initial  location.  Papadakis  has  calculated 
the  diffraction  loss  in  terms  of  normalized  distance  S so  that  the  data  for  all 
wave  lengths  and  transducer  radii  (r)  can  be  reduced  to  one  curve  for  any 
given  anisotropy.  His  curves  for  b * -1.0  to  b = 0.4  are  shown  in  Figure  4. 
The  diffraction  loss  curves  are  all  similar  in  shape  for  different  values  of  b. 
The  loss  fluctuates  with  peaks  (A  = peak)  and  valleys  In  the  Fresnel  region 
before  monotonicaliy  increasing  for  large  values  of  normalized  distance 
S.  The  diffraction  loss  will  be  calculated  from  the  slope  of  these  curves  for 
large  values  of  S.  The  magnitude  of  the  diffraction  loss  is  a strong  function  of 
material  anisotropy  factor  b.  To  minimize  the  diffraction  loss  the  value  of  b 
should  be  positive  (0  to  0.4),  because  these  curves  have  the  minimum  slope. 


8s  z X /r  * 


Figure  4.  Acoustic  Diffraction  Loss  Versus  Normalized  Distance 
for  Values  of  Anisotropy  Parameter  b 

The  b paremater  is  calculated  from 


where  V is  the  velocity  of  the  propagating  acoustic  wave  and  Q is  the  angle 
between  the  propagation  vector  and  the  transducer  normal  N,  If  b is  neg- 
ative, the  energy  propagation  vector  Rv  (Poynting  vector)  makes  a larger 
angle  with  the  transducer  normal  than  does  the  propagation  vector.  The  acous- 
tic energy  surface  spreads  more  rapidly  than  does  the  phase  surface,  the 
acoustic  beam  spreading  is  greater,  and  the  diffraction  loss  is  higher.  If 
b = 0,  the  material  is  elastically  isotropic  and  vectors  /3,  R,  and  N are  colinear. 
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If  b is  positive,  vector  R lies  between  vector  ft  and  the  transducer  normal  N, 
the  acoustic  energy  is  focused  toward  the*  transducer  normal  direction,  and 
beam  spreading  and  hence  diffraction  loss  is  reduced.  In  surface  wave  de- 
vices this  effect  has  been  called  beam  steering.  Papadakis  has  shown  that 
physical  reasons  exist  (Ref.  6)  to  limit  b to  less  than  0.5. 

The  b parameter  is  calculated  from  AV/V,  which  can  be  written  in  terms 
of  the  elastic  constants  of  the  delay  material.  Papadakis  has  calculated  b for 
longitudinal  propagation  along  the  threefold  axis  for  trigonal  crystals  (Z  prop- 
agation in  Ala03)  as 


b 


(C33-Ci3-2C4i)(C33  +C13) 
2C3j(Cj3  - C44) 


(8) 


where  Cw  = 1130,  CM  1 4902,  and  = 1454  for  Ala03.  Using  these  values, 
b is  calculated  for  AlaOa  as  b *=  0. 154.  Since  no  one  has  calculated  the  acoustic 
diffraction  for  shear  waves,  the  diffraction  loss  will  be  estimated  by  calculating 
the  b parameter  from  the  A V/V  for  shear  waves  given  by  Waterman  (Ref.  7) 
and  the  diffraction  loss  curves  shown  in  Figure  4.  In  this  case  the  beam 
spreading  is  not  conical,  and  a complete  numerical  calculation  would  have  to 
include  asymmetric  diffraction  with  N,  R,  and  [i  no  longer  coplanar.  For  shear 
waves  along  the  [1001  in  cubic  materials  (MgAl?04  and  YAG)  the  shear  waves 
are  degenerate  with  a velocity  (C^  p)V*.  The  three  vectors  ft  R,  and  N are 
not  coplanar  but  depend  on  the  angles  9 and  cp(Ref.  7).  The  b parameter  in 
terms  of  elastic  constants  is 


b(0)  * j Ka  i*Kl  -sin^teKa-l)]1/*  j (9) 

where 

Kj  * Cji  - Cia  -2C44 


Ka 


C44-C 


For  MgAla04,  Cu  = 2.79,  C18  = 1.  53.  Cw  = 1.  53;  for  YAG,  Cn  * 3.  33,  Cia  = 

1.  11,  C,^  - 1.  15.  Since  the  angle  cp  is  not  specified,  b was  calculated  for  sin8 
2cp  = 0 and  sin8  2cp  5 1 and  the  results  were  averaged.  Equation  9 gives  two 
values  of  b.  Using  the  elastic  constants  for  MgAl804,  b = 2.  02  and  b = 0 for 
sin8  2cp  = 0,  and  b = 1.  72  and  b = 0.  294  for  sin8  2cp  = 1.  Since  b values  greater 
than  0.  5 are  not  allowed,  b varies  between  0 and  0.  294.  An  average  value  of 
b = 0,  15  was  used  for  shear  waves  along  the  [100]  axis  in  MgAlB04.  Similarly 
for  YAG,  b * 0.  071  and  b = 0 for  sin9  2(fi  = 0 and  b = 0.  054  and  0.  174  for  sin9 
2cp  - 1 . A value  of  b - 0 was  used  for  shear  waves  along  the  [100]  axis  in 
YAG,  which  is  almost  an  elastically  isotropic  material. 
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From  the  b parameter  values  for  longitudinal  and  shear  waves  in  the  de- 
lay materials  of  interest,  one  can  calculate  the  acoustic  diffraction  from  the 
curves  given  by  Papadakis  and  snown  in  Figure  4. 

Relative  loss,  Rt,  is  plotted  there  as  a function  of  normalized  propaga- 
tion distance  S for  several  values  of  the  anisotropy  parameter  b.  At  the  or- 
igin successive  curves  have  been  offset  by  1 decibel  for  clarity;  all  curves 
have  an  actual  value  of  0 decibels  at  the  origin.  The  normalized  propagation 
distance  S is  related  to  the  actual  propagation  distance  z by 

S * z x/ra  (10) 

where  Xis  the  wavelength  of  the  sound  wave  and  r is  the  transducer  radius. 

In  dealing  with  delay  lines  of  length  L,  one  can  observe  signal  outputs  only 
for  propagation  distances  which  are  integral  multiples  of  L.  For  a double 
ended  delay  line  the  nth  observable  pulse  occurs  for  a sound  wave  which  has 
traveled 

z ■ nL,  or  S„  « nLX/r8  (11) 

and  for  a single-ended  delay  line  the  nth  observable  pulse  has  traveled 

z 3 2nL,  or  SB  ■ 2nL\/r8  (12) 


The  propagation  distance  X between  observable  pulses  is  x 9 2L  for 
single -ended  lines  and  X » L for  double-ended  lines.  In  this  study  only  long 
delay  times  which  correspond  to  large  S values  are  of  interest.  The  atten- 
uation due  to  diffraction,  <?„,  is  estimated  by  fitting  a str  aight  line  to  the  ap- 
propriate diffraction  curve  in  Figure  4 as  is  shown  for  the  curve  b 9 0.0. 


at> 


« Rt  fa)  “ Rl  (*n) 

z(n)  - z(m) 


(13) 


The  attenuation  per  reflection  due  to  diffraction,  r„,  is 

^ , Rt  (n)  - Rt  (m) 

r>  ^ — (i4) 

Using  Equations  12  and  14,  r6  may  be  rewritten 

r„  9 [ — s '."sR|>  (m)  ] 6 " 1 for  DE  line;  6 « 2 for  SE  line 

(15) 


where  the  constant  in  brackets  is  the  slope  of  a line  fit  at  large  S to  the  ap- 
propriate diffraction  curve  of  Figure  4 . The  linear  fit  to  the  actual  diffraction 
loss  ignores  the  fluctuations  which  occur  in  the  Fresnel  region  and  extrapolates 
at  the  origin  to  a value  which  is  in  error  by  a few  decibels.  This  intercept 
error  is  small  in  comparison  to  the  total  diffraction  loss  of  the  last  observable 
pulses  and  has  little  effect  on  the  calculation.  Even  so,  it  is  taken  into  con- 


16 


sideration  as  a small  correction  in  these  calculations  by  adding  it,  with  ap- 
propriate sign,  to  the  insertion  loss.  For  low  diffraction  loss  a delay  material 
should  have  low  sound  velocity,  a large,  positive  b value,  and  should  be  avail- 
able in  large  size  so  that  transducers  of  as  large  a diameter  as  possible  can 
be  used.  The  diffraction  loss  will  decrease  as  the  frequency  increases.  When 
the  delay  material  is  specified  and  the  transducer  diameter  is  chosen,  Equa- 
tion 15  can  be  normalized  to  a convenient  frequency  (200  MHz).  Then  the 
diffraction  loss  per  reflection  can  be  written  as 

Fe  a C dB/  reflection  (16) 

where  f is  the  frequency  in  MHz  and  C is  a constant. 

E.  DELAY  LINE  DESIGN 

The  design  of  long  reverberation  time  delay  lines  can  be  approached  in 
two  ways.  One  way  is  to  select  the  materials  from  Tables  1,  2,  and  3 and  to 
decide  on  the  thickness  of  the  top  electrode,  transducer,  and  bond  electrode; 
the  clamped  capacitance  of  the  transducer;  and  the  driving  and  load  resistance. 
Then  using  the  parameters  from  the  tables  with  an  equivalent  circuit,  as  de- 
scribed in  Reference  1,  calculate  the  IL  versus  frequency  response.  Having  the 
IL,  calculate  the  losses  as  described  above  under  "Basic  Lobb  Equations" 
and  divide  the  losses  per  reflection  into  the  available  power,  to  obtain  the 
number  of  reflections  and  hence  the  total  delay  time. 

The  other  approach  is  to  select  only  the  delay  line  material  and  the  elec- 
trode material,  assume  values  for  the  IL,  and  calculate  the  losses  and  total 
delay  time  at  one  frequency.  Plot  the  results  and  find  the  IL  that  gives  the 
maximum  delay  time;  then  design  the  transducer  to  have  that  IL  using  the  equiv- 
alent circuit  computer  program.  Tho  second  approach  will  be  followed  here, 
since  it  appears  to  be  more  straightforward.  Because  the  delay  line  materials 
to  be  compared  have  different  velocities  of  propagation,  the  only  fair  compar- 
ison fs  on  the  basis  of  a fixed  interpulse  period.  It  is  assumed  that  all  the 
delay  lines  have  an  interpulse  period  of  6 ps. 

Longitudinal  Waves  in  AUQ»  at  200  MHz 

Consider  a DE  delay  line  with  Al#Oa  as  the  delay  material  and  gold  elec  - 
trodes.  The  travel  time  per  reflection  is  3 ps  even  though  every  other  echo  is 
read  out.  The  delay  line  length  is  L * t»  V0  * 3.33  cm  long.  Assuming  values 
of  IL  from  30  to  60  decibels,  calculating  r from  Equation  4 and  using  these 
values  for  r,  one  can  calculate  the  acoustic  conversion  loss  per  reflection 
from  Equation  5.  These  acoustic  conversion  loss  values  will  be  used  for  all 
the  delay  lines  to  be  considered. 

The  electrical  skin  depth  for  gold  at  200  MHz  is  4.5  pm.  For  good  elec- 
trical isolation  a gold  bond  electrode  about  5 pm  thick  should  be  used,  How- 
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ever,  the  acoustic  loss  is  excessive  for  shear  waves  in  5-pm -thick  gold  elec- 
trodes and  some  IL  variation  is  evident  in  the  transducer  pass  band  due  to  this 
thick  gold  electrode.  The  electrical  resistance  of  a l-/Ltm-thick  gold  electrode 
0.4  inch  in  diameter  is  only  0.00078  ohm,  so  extra  series  resistance  is  not 
a detrimental  factor.  Transducers  can  be  bonded  with  electrodes  of  1 pm  thick- 
ness; therefore  only  1 -pm -thick  gold  electrodes  will  be  considered  for  all  the 
delay  line  designs.  The  acoustic  loss  per  reflection  in  a 1-pm -thick  gold  elee- 
trode  is  given  by  Equation  6,  with  A 2 200  dB/cm  at  1 GHz  for  longitudinal 
propagation; 

IS  » 0.04  or  0.0016  dB/ reflection  at  200  MHz 

The  transducer  material  loss  is  so  low  that  its  contribution  to  the  total 
loss  can  be  neglected,  as  this  calculation  will  show.  The  transducer  material 
with  the  highest  loss  is  quartz  (3  dB/cm  at  1 GHz).  The  maximum  thickness 
of  the  transducer  is  10  and  there  are  two  passes  through  the  transducer 
per  reflection.  Thus  the  transducer  material  loss  is  2.4  x 10“4  dB  per  re- 
flection at  200  MHz.  Even  with  1000  reflections  thiB  loss  Is  only  0,24  dB, 
which  is  insignificant  compared  to  the  other  losses.  (See  Table  4 for  conver- 
sion loss  in  transducers.) 


Table  4 

ACOUSTIC  CONVERSION  LOSS  PER  REFLECTION 
IN  TRANSDUCERS 


IL  dB 

r 

f t dB/ reflection 

-30 

0. 03162 

0.  1305 

-35 

0. 01778 

0. 0760 

0.  01000 

0.0436 

-4  5 

0. 005623 

0.0244 

-50 

0.  003162 

0. 0139 

-55 

0. 001776 

0. 0078 

-60 

0.  001000 

0.0043 

The  material  loss  of  AlaOa  is  calculated  from  Equation  7,  with  r'  *3ps 
per  reflection  and  B = 0.2  dB/ps.  Then  the  material  loss  per  reflection  is 

r*  ■ 0.6  ) * or  0,024  dB/ ruction  at  200  MHz 

The  acoustic  diffraction  loss  is  expected  to  dominate  at  200  MHz  for  long 
reverberation  time  delay  lines.  It  is  therefore  important  to  decide  what  is 
the  maximum  possible  crystal  diameter  for  the  delay  material.  Good  acoustic 
quality  crystals  of  AlsO„  of  YAG,  and  possibly  of  MgAla04  can  be  obtained  up 
to  a diameter  of  0,4  inch.  The  maximum  diameter  of  0.328  inch  was  chosen 
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for  the  top  transducer  electrode  for  these  calculations.  The  h parameter  was 
previously  calculated  for  Ala03  by  Equation  8 as  b * 0*154.  The  slope  of  the 
diffraction  loss  curve  for  b = 0. 15  in  Figure  4 gives  the  following  values: 

Rt(n)  = 6.31  decibels,  Rt(m)  =4.87  decibels.  SB  = 5,  and  S.  = 3.  Inserting 
these  values  in  Equation  10  along  with  L = 3.33  cm,  V0  - 11.1  xl08cm/s, 
aa  * 0. 694  cm8,  and  f = 200  MHz  the  diffraction  loss  is  calculated  to  be 

rB  = 0.0767  ~2.  ( or  0.07  67  dB/ reflection  at  200  MHz 

This  loss  is  rather  high  because  the  crystal  is  long  and  the  velocity  is  high  for 
A1„03  even  though  the  b parameter  is  in  the  favorable  range. 

With  all  the  losses  now  calculated,  the  total  delay  time  can  be  determined. 

At  anIL  of  -30  decibels  the  available  power  is  80  decibels  (110  dB  dynamic 
range  less  the  IL)  to  compensate  for  the  losses.  The  total  losses  at  200  MHz 
are  acoustic  conversion  loss  (0.1395)  plus  electrode  loss  (0.0016)  plus  material 
loss  (0.0240)  plus  diffraction  loss  (0.767)  which  equals  0.2418  decibel  per  re- 
flection. Dividing  the  total  loss  into  the  available  power  gives  the  total  number 
of  reflections:  80  * 0.2418  ■ 331  reflections.  Multiplying  the  total  number  of 
reflections  by  the  travel  time  per  reflection  (3  ps  per  reflection)  gives  the  total 
delay  time  of  993  ps.  For  the  next  IL  (-35  dB)  the  available  power  is  75  decibels 
and  all  the  losses  remain  the  same  except  the  acoustic  conversion  loss,  which 
decreases  to  0.078  decibel  per  reflection.  The  total  loss  becomes  0.1803  deci- 
bel per  reflection,  and  the  total  delay  time  is  1248  p seconds.  These  calcu- 
lations are  repeated  for  the  other  IL  values  and  the  results  plotted  as  curve  1 
of  Figure  5. 

This  curve  for  a DE  line  shows  several  interesting  points.  The  maximum 
total  delay  time  (1.  55  ms)  occurs  at  an  IL  of  50  decibels.  Thus,  to  design  long 
reverberation  time  delay  lines  the  IL  must  be  rather  high.  At  50  decibels  IL 
the  acoustic  loss  is  dominated  by  diffraction  and  material  loss,  which  limits 
the  total  delay  time.  The  curve  is  rather  flat  from  40  to  55  decibels  IL. 

Next  consider  a SE  delay  line  with  A1803  as  the  delay  material  and  a gold 
electrode.  The  travel  time  per  reflection  is  6 ps,  and  every  echo  is  read  out. 

The  acoustic  conversion  loss  and  the  elect  'e  loss  are  the  same  as  in  the 
DE  case.  The  material  loss  and  the  diffrac.. loss  are  both  doubled.  At  an 
IL  of  -30  decibels  the  available  power  is  80  decibels.  The  total  losses  at 
200  MHz  are  acoustic  conversion  loss  (0.1395)  plus  electrode  loss  (0.0016) 
plus  material  loss  (0.0480)  plus  diffraction  loss  (0.1534),  which  equals  0.3425 
decibel  per  reflection:  80  -r  0.3425  = 233  reflections,  and  multiplying  by  6 ps 
per  reflection  gives  a total  delay  time  of  1398  ps.  The  total  delay  time  is  cal- 
culated for  the  other  IL  values  and  the  results  plotted  as  curve  2 in  Figure  5. 

The  curve  is  similar  to  curve  1 except  that  the  maximum  total  delay  time  is 
1.72  ms  and  occurs  at  an  IL  of  45  decibels.  The  SE  delay  line  always  has  a 
greater  total  delay  time  than  the  DE  delay  line,  because  the  total  number  of 
reflections  is  less  and  thus  the  acoustic  conversion  and  the  electrode  losses 
are  less  compared  to  the  DE  delay  line. 
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Figure  5,  Total  Delay  Time  Versus  Insertion  Loss  at  200  MHz 
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It  ciui  now  be  seen  that  the  transducer  should  be  designed  to  have  an  IL  of 
from  40  to  55  decibels  for  maximum  total  delay  time.  The  equivalent  circuit 
computer  program  will  be  used  to  design  two  transducers  (longitudinal  mode 
ZnO  and  CdS)  on  AlaOa  to  show  how  the  IL  varies  with  frequency  and  how  the 
designer  can  change  the  response  to  achieve  the  desired  IL.  Consider  longi- 
tudinal mode  ZnO,  which  has  a good  acoustic  impedance  match  to  A1803  (ratio 
Zj Zt  = 1.23).  The  transducer  thickness  is  given  by  t * Vj/2f0,  where  Vt  is 
the  longitudinal  mode  velocity  for  ZnO  from  Table  1 and  f0  is  the  desired 
center  frequency.  For  f0  = 200  MHz  the  transducer  thickness  is  15.8  pm. 

The  maximum  transducer  bandwidth  is  obtained  when  the  transducer  is  elec- 
trically matched  to  the  driving  and  load  resistance  (R0  s 50  ohms)  (Ref.  8). 

This  occurs  when  wC0R0  = lj  thus  the  clamped  capacitance  of  the  transducer  is 
C0  = 15.9  pf  at  200  MHz. 

The  top  electrode  thickness  was  selected  to  be  lOOOA  of  gold,  the  bond  electrode 
thickness  1 pm  of  gold,  and  the  velocity,  acoustic  impedance,  and  coupling  coef- 
ficient were  taken  fromTables  1,  2,  and  3.  The  two-way  IL  was  computed  for  these 
parameters  from  60  to  340  MHz  and  the  results  are  shown  as  curve  1 of  Fig- 
ure 6.  A single  minimum  IL  occurs  at  140  MIIz,  which  is  characteristic  of 
acoustically  matched  transducer-propagation  material.  The  minimum  IL  does 
not  occur  at  200  MHz  because  the  transducer  is  mechanically  loaded  by  the  de- 
lay line  material.  The  minimum  IL  can  be  moved  to  200  MHz  by  decreasing 
the  transducer  thickness  to  10.5  pm,  which  increases  the  transducer  f„  to 
300  MHz.  The  results  are  shown  in  curve  2 of  Figure  6.  The  minimum  IL 
has  shifted  to  200  MHz  and  has  a value  of  16,5  decibels,  which  is  less  than 
before  because  this  is  also  the  electrical  matching  frequency  for  a C0  » 15.9  pf. 
This  is  the  minimum  IL  that  can  be  expected  for  this  transducer  with  a kt  «=  0.22. 
The  3-decibel  bandwidth  is  47  percent. 


Figure  6.  Insertion  Loss  Versus  Frequency  for  ZnO-L  on  A1203 
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The  transducer  response  of  curve  2 is  not  a good  design  to  use  for  long 
reverberation  time  delay  lines.  The  IL  is  too  low;  this  would  limit  the  total 
delay  time  because  of  high  acoustic  conversion  loss.  Furthermore,  for  a 
10. 5 -pm -thick  ZnO  transducer  the  top  electrode  diameter  is  0.065  inch  for  a 
clamped  capacitance  of  15.9  pf;  this  would  result  in  an  enormous  acoustic 
diffraction  loss.  The  only  parameter  the  designer  can  change  at  this  point 
is  the  clamped  capacitance  C0.  With  an  increase  in  the  top  electrode  diam- 
eter to  0.164  inch,  the  C0  becomes  101.5  pf  and  the  IL  versus  frequency 
response  is  shown  as  curve  3 in  Figure  6.  The  minimum  IL  at  200  MH2  has 
increased  10  decibels  but  so  has  the  VSWE  of  the  device.  These  devices  in 
general  will  all  have  a high  VSWR.  Increasing  the  top  electrode  diameter  to 
0.328  inch,  the  value  used  for  the  loss  calculations,  raises  the  C0  to  403  pf; 
the  results  are  shown  as  curve  4 in  Figure  6.  The  minimum  IL  at  200  MHz  is 
38  decibels  and  the  IL  varies  4. 5 decibels  over  the  desired  bandwidth  from 
150  to  250  MHz.  By  varying  the  transducer  thickness,  the  IL  can  be  made  to 
decrease,  remain  flat  (as  shown  in  curve  4),  or  Increase  as  a function  of  fre- 
quency over  the  desired  bandwidth.  The  authors  consider  the  IL  response  of 
curve  4 to  be  an  adequate  transducer  design  for  a long  reverberation  time 
delay  line. 

The  second  transducer  design  is  longitudinal  mode  CdS,  which  has  an 
acoustic  impedance  mismatch  to  Al#Os  (ratio  Z0/ Zl  ■ 2.06).  The  CdS  trans- 
ducer thickness  for  f0  B 200  MHz  is  11.15  pm.  The  same  parameters  were 
used  as  in  the  previous  transducer  design:  C0  B 15.8  pf,  R0  * 50  ohms  for  uC0R0  »1 
at  200  MHz,  top  electrode  thickness  1000A  of  gold,  the  bond  electrode  thickness 
lpm  of  gold,  and  the  appropriate  velocity,  acoustic  impedance,  and  coupling 
coefficient  from  Tables  1,  2,  and  3,  The  results  are  shown  in  curve  1 of  Fig- 
ure 7.  The  IL  response  now  has  a double  minimum  at  110  MHz  and  290  MHz, 
which  is  characteristic  of  acoustically  mismatched  transducer -propagation 


Figure  7,  Insertion  Loss  Versus  Frequency  For  CdS-L  On  AlpOa 
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material.  The  II-  of  CdS  transducers  is  higher  than  ZnO  transducers  because 
of  the  lower  value  of  kt  for  CdS.  Increasing  the  top  electrode  diameter  to 
0.164  inch,  C becomes  103  pf  and  again  the  IL  has  increased  10  decibels, 
as  shown  in  curve  2 of  Figure  7.  By  increasing  the  top  electrode  diameter  to 
0.328  inch,  the  value  used  for  the  loss  calculations,  C0  * 411  pf  and  the  IL 
is  55  decibels  at  200  MHz,  as  shown  in  curve  3 of  Figure  7.  The  IL  variation 
of  curve  3 is  12  decibels  from  150  to  250  MHz.  This  IL  variation  could  be 
reduced  by  changing  the  transducer  thickness,  but  that  will  not  be  necessary, 
as  will  be  seen  when  the  total  delay  time  variation  across  the  transducer  band- 
width is  calculated. 

The  next  step  is  to  go  back  and  recalculate  the  SE  delay  line  losses,  using 
the  actual  values  for  the  IL  at  several  frequencies  across  the  transducer  band 
width,  to  recalculate  the  acoustic  conversion  loss  using  the  computer  calculated 
values  of  IL,  and  to  account  for  the  frequency  squared  variation  of  electrode 
and  material  losses  and  the  reciprocal  frequency  variation  of  the  diffraction 
loss.  At  the  frequencies  marked  on  curve  4 in  Figure  6 for  the  ZnO-L  AlgOs 
SE  delay  line,  the  total  delay  time  variation  is  as  shown  in  curve  1 of  Figure  8. 
The  total  delay  time  variation  is  13.8  percent.  At  the  frequencies  marked  on 
curve  3 in  Figure  7 for  the  CdS-L  Alg03  SE  delay  line,  the  total  delay  time 
variation  is  shown  as  curve  2 in  Figure  8.  There  the  total  delay  time  vari- 
ation is  only  4.4  percent.  The  IL  variation  in  this  case  has  conveniently  com- 
pensated for  the  frequency  variation  of  the  electrode  material  and  diffraction 
losses.  Other  methods  of  total  delay  time  compensation  such  as  variable 
system  gain  could  be  considered,  but  that  involves  system  details  which  are 
better  left  to  the  system  designer.  The  point  to  be  stressed  here  is  that 
there  is  considerable  flexibility  in  the  transducer  design,  depending  on  what 
one  wants  to  design  for. 

Shear  Waves  in  MgAl,QA  at  200  MHz 

Consider  a DE  delay  line  with  MgAlg04  as  the  delay  material  and  gold  elec  - 
trodes.  The  travel  time  per  reflection  is  3 ps,  even  though  every  other  echo 
is  read  out.  The  delay  line  length  is  r'VD  ■ 1.98  cm.  Values  are  again  as- 
sumed for  IL,  and  r calculated  from  Equation  4 and  the  acoustic  conversion 
loss  from  Equation  5 as  before. 

The  acoustic  loss  per  reflection  in  a 1 -pm -thick  gold  electrode  is  given 
by  Equation  6 with  A ■ lOOOdB/cm  at  1 GHz  for  shear  propagation 

H ■ 0.2  * or  dB/ reflection  at  200  MHz 

The  material  loss  of  MgAl^O-  is  calculated  from  Equation  7,  with  t'  ■ 3 ps 
per  reflection  and  B = 0,04  dB/ps.  The  material  loss  per  reflection  is 

r«  = 0.12  or  0.0048  dB/ reflection  at  200  MHz 
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Figure  8.  Total  Delay  Time  Variation  Across  Transducer  Bandwidth 
in  SE  Delay  Lines 

The  acoustic  diffraction  loss  parameter,  b,  was  previously  calculated  for 
MflAljO*  using  Equation  0 as  b - 0 to  0.  204.  An  average  value  of  b ■ 0.  IB  was 
estimated  for  shear  wave  propagation  along  the  [100]  direction  In  Mg  a 4* 
The  slope  of  the  diffraction  loss  curve  for  b • 0.  IB  in  Figure  4 18  “amuaB 
for  A1  CLs  inserting  those  values  in  Equation  10  along  with  L « 1.  00  tm,  Vp 
6 54  cm/.™  a*  . 0.  604  cm',  end  f ■ 200  Mil*,  the  diffraction  lo.a  1. 

calculated  to  be 

rp  ■ 0.0266 1222\,  or  0.0286  dB/ reflection  at  200  MHz 

This  loss  is  a factor  of  three  less  than  the  diffraction  loss  for  A1.A, 
though  both  materials  have  the  same  b parameter,  because  (1)  the  crystal  is 
shorter  for  the  same  travel  time  and  (2)  the  velocity  is  lower. 
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With  all  the  losses  calculated,  the  four  losses  are  summed  and  the  total 
delay  time  calculated  as  before.  The  results  are  plotted  as  curve  3 in  Figure  5. 
The  maximum  total  delay  time  is  3.5  ms  and  ooours  at  an  IL  of  55  decibels. 

At  55  decibels  1L  the  acoustic  loss  la  dominated  entirely  by  diffraction  loss. 

The  rapid  increase  in  the  total  delay  time  at  low  values  of  IL  is  caused  by  a 
similar  rapid  decrease  in  the  acoustic  conversion  loss  as  the  IL  increases. 

The  SE  delay  line  with  MgAla04  delay  material  is  calculated  as  previously 
described,  with  a travel  time  of  6 mb  per  reflection  and  the  acoustlo  conversion 
and  electrode  losses  remaining  constant  while  the  material  and  diffraction 
losses  are  both  doubled.  The  results  are  shown  as  curve  4 in  Figure  5,  Again 
it  is  seen  that  the  SE  delay  line  has  a greater  total  delay  time  than  the  DE  case, 
The  maximum  total  delay  time  is  4.25  mB  and  occurs  at  an  IL  of  50  decibels. 
The  difference  in  total  delay  time  is  greater  between  curves  3 and  4 than  be- 
tween curves  1 and  2 because  the  material  and  diffraction  losses  are  so  much 
lower  in  shear  propagating  MgAla04  than  in  longitudinal  propagating  AlaOa. 

Again  the  transducer  IL  should  be  45  to  55  decibels  for  maximum  total  delay 
time. 

Since  a high  IL  is  desired,  a transducer  material  with  only  s modest 
coupling  coefficient  will  be  selected.  Y-cut  S10a  has  an  acoustlo  Impedance 
mismatch  for  shear  wave  propagation  in  MgAl,04  (ratio  Zj Zl  ■ 2.32)  and  will 
be  entirely  satisfactory  as  a transducer  material.  The  shear  wave  polariza- 
tion direction  in  Y-cut  SiOa  must  be  aligned  with  the  shear  wave  polarization 
in  MgAla04,  but  this  can  be  done  with  careful  x-ray  alignment  toohnlques. 

The  transducer  fe  ■ 250  MHz  has  been  set  to  flatten  the  IL  response  and  the 
transducer  thickness  is  7.6  pm.  The  same  parameters  were  used  as  in  the 
previous  transducer  designs)  C0  ■ 15.0  pf,  RB  * 50  ohms  for  uC0Rp  ■ 1 at 
200  MHz,  top  electrode  thickness  1000A  of  gold,  the  bond  electrode  thiokness 
1 pm  of  gold,  and  the  appropriate  velocity,  acoustic  impedance,  and  coupling 
coefficient  from  Tables  1,  2,  and  3.  The  computer  results  are  shown  in 
curve  1 of  Figure  9.  The  IL  is  28  decibels  and  1b  flat,  with  a 3-deoibel  band- 
width from  08  to  262  MHz.  Increasing  the  top  eleotrode  diameter  to  0.328  inch, 
the  value  used  for  the  loss  calculations,  C9  « 2BB  pf  and  the  IL  is  47  decibels  at 
200  MHz,  as  shown  in  curve  2 of  Figure  B.  The  Cs  of  the  SlOa  transducer  is 
less  than  the  C0  of  the  ZnO  or  CdS  transducers,  because  the  dielectric  constant 
of  SiOa  is  4.5  while  it  is  approximately  9 for  ZnO  and  CdS.  The  IL  variation 
of  curve  2 1b  B, 5 deoibels  from  150  to  250  MHz. 

At  the  frequency  points  marked  in  curve  2 of  Figure  0 the  IL  values  were 
used  to  recalculate  the  acoustic  conversion  loss.  The  frequency  variation  of 
the  bond,  material,  and  diffraction  losses  were  accounted  for  and  the  total 
delay  time  calculated  as  before.  The  results  are  shown  as  curve  3 of  Figure  8. 
The  total  delay  time  variation  is  15.7  percent  across  the  transducer  bandwidth 
for  a shear  wave  St09  transducer  on  MgAl#04.  The  lower  total  delay  time  at 
150  MHz  is  caused  by  the  dominating  diffraction  loss,  which  is  not  reduced  by 
the  decrease  in  material  and  electrode  losseB  because  they  are  only  a small 
fraction  of  the  total  losses.  The  IL  response  should  have  a negative  slope 
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Figure  0,  Insertion  Loss  Versus  Frequency  For  SiO,-S  On  MgAlg04 

across  the  transducer  bandwidth  to  help  flatten  the  total  delay  time  variation. 

No  further  attempt  was  made  to  modify  the  transducer  response  In  order  to 
Improve  the  total  delay  time  variation. 

The  IL  response  with  a SiO,-S  transducer  on  YAG  is  almost  identical  to 
the  IL  response  with  MgAlp04  delay  material,  because  the  acoustlo  impedances 
of  YAG  and  MgAl,04  nre  so  similar.  Other  transducer  materials  were  also 
tried  (ZnO-S  and  LiNbOj-S),  but  because  of  the  high  dielectric  oonstant  of 
these  materials  the  C0  was  very  high,  and  thus  they  did  not  match  well  with  the 
50 -ohm  driving  and  load  resistance, 

Shear  Waves  In  MflAliOtSE  Delay  Line  with  a Swltohrd  Impedance 
Transducer  at  200  MHz  ~ 

All  the  acoustic  losses  for  shear  waves  in  a MgAl„04  SE  delay  line  having 
boon  calculated,  the  problem  remaining  is  the  transducer  design  for  a switched 
impedance  delay  line.  The  idea  here  is  to  have  a transducer  with  a low  IL  to 
drive  the  delay  line  (two-way  IL  of  10  to  20  dB),  then  switch  the  transducer  to 
a high  IL  for  readout  (two-way  IL  of  45  to  55  dll).  A transducer  material  with 
a high  ooupling  coefficient  is  needed  to  drive  the  delay  line.  Consider  a shear 
mode  ZnO  on  MgAlt04  as  the  SE  switched  delay  line.  The  transducer  f„  ia 
increased  to  250  MHz  to  flatten  the  IL  response  at  200  MHz.  The  transducer 
thickness  is  given  by  t - V J2ta,  where  V,  is  thr  shear  wave  velocity  in  ZnO 
from  Table  1 and  f0  is  the  transducer  resonant  frequency  250  MIIz  in  this  case. 
The  transducer  thickness  is  5.44  pm.  The  clamped  capacitance  is  calculated  from 


Co  = 


8.854  x 10“u  x c x n x a* 
r 

4 x t 


where 


c r * relative  dielectric  constant  of  ZnO 
a ■ top  electrode  diameter  (0.833  cm) 
t B transducer  thickness  (cm) 

The  clamped  capacitance  C„  ■ 736  pf  and  the  capacity  reaotanoe  at  200  MHz 
is  1.08  ohms.  The  transducer  will  have  a minimum  IL  when  uC„Ra  ■ J at 
200  MHz,  which  is  the  electrical  matching  condition.  Thus  R0  should  be  1.08 
ohms  for  minimum  IL.  A driving  resistance  of  2 ohms  can  be  used  with  only 
a small  Increase  in  IL.  The  stepdown  impedance  transformer  haB  an  imped- 
ance ratio  of  25:1. 

The  remaining  parameters  are  the  same  as  before:  top  electrode  thick- 
ness, 1000A  of  goldj  the  bond  electrode  thickness,  1 pm  of  gold;  and  the  ap- 
propriate velocity,  aooustic  impedance,  and  coupling  coefficient  from  Tables 
1,  2,  and  3.  The  computed  two-way  IL  is  shown  as  curve  1 of  Figure  10. 


Figure  10.  Insertion  Lobh  Versus  Frequency  For  ZnO-S  On 
MgAlB04  Switched  Impedance  Delay  Line 
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The  IL  is  13.5  decibels  at  200  MHz,  with  a 3-decibel  bandwidth  of  45  percent. 

If  the  driving  and  load  resistance  is  switched  to  200  ohms,  the  IL  response  is 
as  shown  in  curve  2 of  Figure  10.  The  IL  is  now  50.5  decibels  at  200  MHz, 
with  a very  flat  IL  response.  The  stepup  Impedance  transformer  required 
has  an  impedance  ratio  of  4tl  for  a 50-ohm  system.  The  IL  for  the  switched 
SE  delay  line  was  shown  in  Equation  3 to  be  r^r,,  which  can  be  calculated  from 
the  computer  two-way  IL  curves  as  1/ 2 <IL)*.,+  1/2  (IL)Ml00  ■ 32  decibels 
at  200  MHz.  Thus  the  available  power  is  78  decibels  for  the  SE  switched  delay 
line,  as  compared  to  60  decibels  for  the  nonswitched  SE  delay  line.  With  18 
decibels  more  power  to  compensate  for  the  losses,  the  total  delay  time  for  the 
SE  switched  delay  line  Is  5.5  ms,  as  shown  at  point  5 In  Figure  5,  an  Increase 
of  2.3  ms  over  a similar  nonswitched  SE  delay  line.  If  the  bond  electrode  loss 
could  be  eliminated,  the  total  delay  time  would  increase  to  6,1  ms  (point  6 in 
Figure  5).  The  advantage  of  the  switched  SE  delay  line  is  Its  lower  overall  IL, 
and  thus  higher  available  power  for  acouBtic  Iobs  compensation.  The  required 
Impedance  transformers  are  realizable,  and  it  should  be  possible  to  build  a 
practical  switched  SE  delay  line  using  this  concept.  A S-LiNbO,  transducer  on 
MgAla04  was  also  simulated  in  the  switched  configuration,  but  the  C0  value  Is 
so  high  that  electrical  Impedance  matching  becomes  impractical. 

Figure  11  shows  a comparison  of  all  the  materials  considered.  In  the  form 
of  curves  of  total  dolay  time  versus  frequency  for  SE  delay  lines.  Only  the 
results  for  SE  delay  lines  were  plotted,  because  that  configuration  gives  the 
maximum  delay  time.  These  curves  were  computed  on  the  assumption  of 
constant  IL  for  all  frequencies,  using  the  IL  value  that  gave  the  maximum 
total  delay  time  at  200  MHz  (46  dB  IL  for  longitudinal  and  50  dB  IL  for  shear 
waves).  Transducers  could  be  designed  to  obtain  these  IL  values  at  the  higher 
frequencies t it  Is  felt,  however,  that  the  Iosbbs  change  with  frequency  and  it  is 
not  certain  that  the  total  delay  time  will  peak  at  these  IL  values  for  all  fre- 
quencies. The  electrode  material  and  diffraotion  losses  were  all  corrected 
for  their  frequency  variationj  only  the  acoustic  conversion  loss  was  assumed 
to  be  constant. 

All  the  curves  have  a similar  shape,  with  the  total  delay  time  decreasing 
as  reciprocal  frequency  squared  at  high  frequencies.  This  is  not  surprising, 
since  electrode  and  material  losses  are  proportional  to  frequency  squared  and 
these  losses  dominate  at  high  frequencies.  The  diffraction  loss  can  be  ne- 
glected above  800  MHz  for  the  transducer  top  electrode  diameter  used  in 
these  calculations.  The  maximum  total  delay  times  peak  at  low  frequencies 
(200  to  300  MHz)  and  have  a range  of  values  from  1.7  ms  for  L-AltOa  to 
6.6  ms  for  S-MgAlt04  (switched).  The  effect  of  electrode  loss  is  clearly 
shown  in  the  curves.  If  the  electrode  loss  is  eliminated,  the  total  delay 
time  will  increase  about  0,5  ms  for  YAG,  1.0  ms  for  MgAlaOt,  and  1.5  ms 
for  MgAl,04  (switched).  These  increases  are  not  constant  but  depend  on  fre- 
quency becoming  less  at  high  frequencies  and  greater  at  low  frequencies.  The 
ranking  of  the  delay  line  materials  is  just  what  one  would  expect  it  to  be.  At 
1 GHz  the  higher -material -Iobs  L-Al,Oa  has  the  lowest  total  delay  time;  lower- 
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Figure  11.  Total  Delay  Time  Versus  Frequency 
For  SE  Delay  Lines 
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loss  S-YAG  has  a greater  total  delay  time;  and  finally,  lowest-loss  S-MgAl#04 
has  the  highest  total  delay  time. 

In  summary,  the  switched  impedance  SE  delay  line  with  a low -loss  metal 
electrode,  using  shear  waves  propagating  in  MgAlt04,  will  give  the  longest 
possible  delay  time.  It  appears  possible  to  build  a 1-ms  delay  line  up  to  1 GHz 
with  the  switched  impedance  configuration.  If  the  system  requires  less  total 
delay  time,  then  the  SE  delay  line  can  be  used  with  less  circuit  complexity  or 
the  DE  delay  line  can  be  used  for  greater  electrical  isolation.  The  acoustic 
loss  calculations  can  be  readily  repeated  for  lnterpulse  periods  other  than 
the  one  considered  here  (6jus). 


The  optimum  IL  for  long  reverberation  time  delay  lines  has  been  shown 
to  be  45  to  55  decibels.  The  design  of  transducers  with  adequate  bandwidth  to 
obtain  these  IL's  was  discussed.  The  report  also  showed  how  the  IL  response 
of  the  transducer  could  be  modified  to  obtain  flat  IL  versus  frequency  or  flat 
total  delay  time  versus  frequency. 

Finally,  it  was  shown  that  greater  than  1 ms  delay  lines  are  possible  with 
L-A1,03  up  to  500  MHz,  with  S-YAG  up  to  670  MHz,  with  S-MgAlg04  up  to 
870  MHz,  and  with  switched  Impedance  S-MgAlg04  up  to  1000  MHz.  Bond  loss 
can  be  reduced  by  using  other  metals  for  electrodes  so  that  these  results  can 
be  obtained  at  even  slightly  higher  frequencies.  It  will  be  necessary  to  de- 
termine if  these  results  are  experimentally  achievable  by  constructing  several 
types  of  delay  lines. 


It  is  Interesting  to  note  that  the  information  handling  capacity  of  a device 
is  related  to  its  delay  time  -bandwidth  product.  Since  the  bandwidths  are 
nearly  50  percent,  this  parameter  can  be  estimated  by  locating  the  frequency 
at  which  a line  of  slope  unity,  one  decade  decrease  in  t for  one  decade  in- 
crease in  frequency,  is  tangent  to  the  curves  of  Figure  11.  This  is  the  fre- 
quency of  greatest  delay  time -bandwidth  product.  The  frequency  is  near  400 
MHz  and  the  delay  time -bandwidth  product  is  near  10®.  The  value  of  10®  is 
considerably  higher  than  that  of  most  other  devices. 

F.  MEMORY  DEVICES  AND  MATCHED  FILTERS 

In  communication,  radar,  or  sonar  systems,  it  has  been  found  useful  to 
code  signals  which  are  transmitted  and  then  to  decode  them  with  a matched 
filter  (Refs.  9 and  10)  that  haB  the  property  of  providing  maximum  signal -to- 
noise  ratio  in  the  presence  of  white  noise.  The  signal  may  be  coded  to  pro- 
vide secure  communication  channels  or  to  provide  simultaneous  range  and 
velocity  information.  Signals  may  be  coded  by  frequency,  amplitude,  or 
phase  modulation  or  a combination  of  these.  Examples  of  coded  waveforms 
are  described  below  mathematically. 
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Frequency  Modulated  Signal: 


e(t)  s (u(to)  - u(to+A)]cos ^n^t+^kat*  + -Ikst*  + * * 
for  which  the  instantaneous  frequency  f is  time  dependent 

f * ^-(f0t  + + |kst3  + . . .^ 

f » f0  + M + kjt3  + . . . 

The  function  u(t)  is  a unit  step  function  at  time  tD: 


u(t0)  5 0 t < t0 

u(t0)  =1  t > t„ 


and  the  quantity  fu  (t0)  - u(t0  +A) "]  makes  the  frequency  modulated  Bignal  a 
pulse  of  duration  A beginning  at  time  t = te. 


Discrete  Coded  Signal: 


e<t)  ■ &n 


(u(t+(n-lH)  - u(tfnA)Jcos[(uj0+uun)t  + 9n] 


which  consists  of  N signal  parts  that  may  vary  in  amplitude,  An;  frequency, 

i-  <wc+wB); 

2ir 

and  phase,  Qt,  Signals  such  as  these  are  designed  to  meet  specific  system 
requirements, and  this  is  generally  carried  out  by  examining  the  ambiguity 
or  autocorrelation  function  of  the  signal  (Refs.  9 and  10).  A familiar  example 
of  the  use  of  a coded  waveform  is  that  of  linearly  frequency  modulated  radar, 
chirp  radar.  The  signal  is  given  by  Equation  17,  with  value  kj  ■ 0 for  j >2. 
Upon  reception,  the  modulated  signal  is  sent  through  its  matched  filter,  which 
has  the  property  of  being  linearly  dispersive  and  of  delaying  the  frequency  com- 
ponents of  the  signal  by  different  amounts  so  that  they  catch  up  with  one  another 
to  produce  a short  pulse  of  large  amplitude.  The  signals  must  be  designed  to 
meet  system  requirements,  since  the  performance  of  the  system  will  be  de- 
termined by  the  signal  coding.  Once  the  signal  is  specified,  the  matched 
filter  is  determined.  If  the  signal  is  e (t)  with  a spectrum  E(u), 


eftje^dt 


The  the  matched  filter  is  a filter  network  with  impulse  response  h(t)  propor- 
tional to  the  time  inverse  of  the  signal 


h(t)  a ke(-t) 


where  U is  a constant. 


An  equivalent  definition  of  the  matched  filter  is  made  from  the  matched 
filter  transfer  function,  H(w): 

H(u)  = kE*  <w)e  T4  (2 

where  Td  is  a delay  constant,  and  where 


H(w)  = J h(t)e3Wtdt 


Design  of  the  signal  determines  the  system  performance  and  specifies  the 
matched  filter.  In  practice  one  does  not  attempt  to  build  a filter  with  suffi- 
cient bandwidth  to  accommodate  all  of  the  sidelmnds  of  E(w)  as  specified  by 
Equation  24.  Filters  actually  used  are  approximations  to  Equations  23  and 
24,  with  a finite  number  of  stages  and  a finite  bandwidth.  In  some  applications 
it  is  more  important  to  reduce  the  side,  lobes  of  the  filter  output  than  to  main- 
tain maximum  signal-to-nolse  ratio,  and  the  filter  characteristics  are  slightly 
modified  for  this  purpose. 

Summary  of  Matched  Filters 


With  this  brief  introduction  to  matched  filters  as  background,  the  purposes 
of  this  portion  of  this  report  can  be  summarized.  The  practical  advantages 
of  coded  waveforms  and  matched  filters  imply  that  acoustic  delay  devices  must 
be  examined  for  their  usefulness  as  memory  elements  for  coded  waveforms. 
Specification  of  delay  time,  bandwidth,  and  insertion  loss  is  not  sufficient. 
Distortion  of  a coded  waveform  must  also  be  considered.  In  the  remainder  of 
this  section  the  practical  advantages  of  matched  filtering  will  be  listed.  A 
parameter  will  be  introduced  to  specify  the  quality  or  information  handling 
capability  of  a matched  filter,  A description  will  be  given  of  the  ranges  in 
frequency  and  bandwidth  in  which  matched  filters  can  be  constructed  with 
available  technologies.  An  examination  will  be  made  of  the  applicability  of 
bulk-wave  acoustic  delay  lines  and  other  devices  as  memory  elements  for 
coded  signals.  The  relationship  will  be  presented  between  the  choice  of 
materials  for  delay  medium,  bond,  and  transducer  and  the  performance  of 
bulk-wave  delay  lines  as  memory  elements  for  coded  signals. 
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The  advantages  of  matched  filtering  include: 

• Efficient  use  of  available  transmitter  power. 

• Improved  system  capability  --  range  and  velocity  resolution  in  radar 
or  sonar,  information  density,  and  security  in  communications. 

• Reduced  interference  by  unwanted  signals  that  do  not  have  the  proper- 
ties of  the  coded  waveform. 

• Real  time  extraction  of  important  parameters  from  the  received 
signal. 

Grant  et  al  (Ref.  11)  have  recently  described  the  value  of  matched  filtering 
techniques  in  meeting  the  increasingly  stringent  requirements  for  communi- 
cation, identification,  and  surveillance  in  air  traffic  control  systems. 

A useful  parameter  for  specifying  the  information  handling  capability  of 
a matched  filter  is  the  product  of  filter  delay  time  and  bandwidth,  D,  frequently 
called  the  pulse  compression  ratio  in  radar.  It  may  be  understood  by  reference 
to  Figure  12.  A sharp  pulse  of  duration  T,  contains  a wide  spectrum  of  fre- 
quencies. When  such  a pulse  is  applied  to  a filter  that  delays  the  high-fre- 
quency comp6nents  with  respect  to  the  tow-frequency  components,  a time 
stretched  signal,  duration  Ta,  with  frequency  modulation  is  the  output  from 
the  filter.  This  signal  may  be  amplified  and  transmitted.  On  reception  it  can 
be  compressed  again  by  passing  through  a conjugate  filter  which  delays  low- 
frequency  components  with  respect  to  high-frequency  compone  nts.  Bandwidth 
is  required  if  a short  pulse  is  to  be  handled,  and  time  delay  is  required  if  the 
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pulse  is  to  be  stretched  and  recomprcssed.  The  product  of  delay  and  band- 
width, D,  is  equal  to  the  ratio  of  the  durations  of  the  stretched  and  compressed 
pulses,  and  is  a measure  of  the  information  handling  capability  of  the  system. 

An  alternative  way  of  thinking  about  D is  to  note  that  the  range  resolution 
of  a radar  system  is  inversely  proportional  to  D.  Early  radar  systems  used  a 
continuous  signal  for  accurate  velocity  estimation,  with  no  range  information; 
and  a sharp  pulse  for  range  estimation,  with  little  velocity  information.  Such 
a system  can  be  optimized  to  provide  simultaneous  range  and  velocity  resolu- 
tion by  proper  design  of  a coded  signal  and  the  fabrication  of  matched  filters 
for  that  signal.  The  delay  time -bandwidth  product,  D,  is  the  parameter  which 
relates  to  the  product  of  range  resolution  and  velocity  resolution. 

In  recent  years,  two  technologies  have  produced  a new  flexibility  and  op- 
erating range  for  matched  filters  that  may  be  readily  fabricated.  These  are 
the  surface  acoustic  wave  (SAW)  technology  and  silicon  technology.  Figure  13 
summarizes  the  capabilities  of  SAW  technology  today.  The  boundaries  in- 


Figure  13.  Delay  Time  Versus  Bandwidth  for  Surface  Acoustic  Wave 
Devices 
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dicated  will  move  as  the  technology  improves.  Plots  similar  to  Figure  13 
have  been  independently  developed  by  Bell  et  al.  (Ref.  12)  and  by  Martin 
(Ref.  13).  The  horizontal  and  vertical  axes  of  Figure  13  are  bandwidth  and 
delay  time;  consequently  the  dashed  diagonal  lines  are  lines  of  constant  pulse 
compression  ratio  D,  and  several  such  lines  are  drawn  in  and  labeled.  The 
shaded  region  bounded  by  a solid  line  indicates  the  region  in  which  SAW  devices 
can  be  fabricated  today.  The  bandwidth  limitation  is  set  by  the  resolution  capa- 
bilities of  photolithography  or  electron  beam  lithography  in  defining  very  fine 
transducer  finger  pairs.  The  vertical  limits  are  set  by  substrate  size,  uni- 
formity of  sound  velocity  over  the  substrate,  the  information  handling  power 
of  lenses  used  in  processing,  etc.  For  the  very  long  delay  times  shown  at 
lower  frequencies,  closed  paths  or  spiral  paths  around  a substrate  must  be 
used. 

Silicon  technology  has  provided  several  types  of  devices  that  store  elec- 
trical charge  capacitatlvely  and  move  it  from  site  to  site  with  a clock  pulse. 

The  quantity  of  charge  can  be  varied  continuously,  so  these  devices  can  be 
operated  in  either  digital  or  analog  modes.  For  the  purposes  of  this  report 
the  authors  will  not  distinguish  between  junction  devices  (Ref.  14)  (bucket 
brigades,  BBD)  and  MOS  surface  charge  devices  ( (CCD  (Ilrf.  15),  SCT 
(Ref.  16)  ).  To  avoid  the  multiplicity  of  names  associated  with  the  silicon 
devices,  they  will  be  collectively  referred  to  here  as  clocked  annlog  delay 
devices  (CLAD),  In  shifting  charge  from  site  to  site,  some  charge  is  left 
behind  and  the  efficiency  of  the  transfer  degrades  as  the  frequency  increases. 
Consequently,  the  speed  or  bandwidth  of  the  device  is  related  to  the  number  of 
sites,  and  this  is  a limitation  on  the  transfer  of  information  into  or  out  of  the 
device. 

Reasonable  limits  for  the  input  and  output  capabilities  of  CLAD  devices 
are  shown  in  Figure  14  by  the  two  diagonal  lines  that  were  calculated  for  P 
and  N channel  MOS  realizations.  Bucket  brigade  devices  can  be  fabricated  by 
conventional  techniques  with  capabilities  comparable  or  slightly  superior  to 
the  N channel  MOS  devices.  Further  developments  in  silicon  technology  with 
buried  channels  and  very  fine  structures  can  be  expected  to  increase  CLAD 
bandwidths  by  about  a factor  of  two.  In  the  time  delay -bandwidth  region  in 
which  they  function,  CLAD  devices  have  considerable  advantage  over  SAW  de- 
vices in  being  compatible  with  other  silicon  devices,  in  both  signal  level  and 
fabrication  procedures,  and  in  being  clocked  devices  which  can  store  informa- 
tion for  significant  times  until  called  for.  A reasonable  estimate  of  storage 
time  is  given  by  the  horizontal  line  labeled  "CLAD  Memory  Time.”  Tills 
value  is  considerably  below  the  theoretical  limit  and  is  affected  by  differen- 
tial charge  leakage  at  the  storage  sites  which  must  be  controlled  by  processing 
techniques.  CLAD  devices  may  be  used  as  matched  filters  and  as  memory 
elements.  The  dark  circles  show  a reasonable  set  of  goals  for  CLAD  memory 
(Ref.  17). 

The  memory  time  of  a MgAl,04,  spinel  delay  line  is  shown  on  Figure  14 
for  two  different  transducer  geometries.  The  calculations  are  reported  in 
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Figure  14.  Delay  Time  Versus  Bandwidth  for  Silicon  Clocked  Analog  Delay 
(CLAD)  Devices  and  for  Bulk  Acoustic  Wave  Devices 
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detail  in  this  section  of  the  report,  under  "Delay  Line  Design. " Inspection 
of  Figure  14  shows  that  CLAD  devices  should  be  considered  for  applications 
requiring  matched  filtering  or  memory  with  bandwidth  requirements  of  about 
10  MHz  or  less.  Where  greater  bandwidths  are  required,  SAW  devices  provide 
a unique  technology  for  realizing  matched  filters.  Bulk -wave  delay  lines  are 
seen  to  be  suitable  as  memory  elements  in  this  bandwidth  region. 

Another  potentially  useful  memox'y  device  for  bandwidths  above  10  MHz 
is  the  wraparound  surface  acoustic  wave  delay  line.  This  device  has  been 
limited  to  the  cubic  piezoelectric  material  Bi^GeO^  (BGO)  up  to  the  present, 
because  of  mode  conversion  in  other  piezoelectric  materials  bb  the  surface 
acoustic  wave  propagates  around  the  curved  edges  of  the  device.  The  intrinsic 
attenuation  in  BGO  is  greater  than  in  MgAl,04  and  limits  delay  time  for  pas- 
sive devices  as  shown  in  Figure  15, 


Delay  Time  Versus  Band  width  for  Surface 
Acoustic  Wave  Wraparound  Passive  Devices 


Figure  15. 
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It  is  possible  to  compensate  for  these  losses  by  acoustic  amplification,  but 
only  at  the  cost  of  introducing  noise  and  dispersion.  To  the  authors'  knowledge 
the  limitations  of  active  SAW  wraparound  delay  lines  due  to  noise  and  dispersion 
have  not  been  characterized  for  memory  of  coded  signals.  Similarly,  the  per- 
formance of  bulk-wave  acoustic  delay  lines  as  memory  elements  for  coded  sig- 
nals is  markedly  affected  by  the  impedance  match  between  the  delay  medium 
bond  and  transducer.  This  problem  is  considered  in  detail  below.  The  present 
status  of  the  several  types  of  devices  just  discussed  is  summarized  in  Table  5. 

Table  3 

COMPARISON  OP  MATCHED  FILTER  AND  MEMORY  DEVICES 


Function 

Device 

llandwidth  (MHz) 

Problems 

Matched  filter 

CHAD 

£ 10 

. ■ 

Matched  filter 

RAW 

>10 

-- 

Memory 

CLAD 

£ 10 

Differential  charge  leakage 

Memory 

SAW 

> 10 

Noise  buildup,  dispersion 

Memory 

Hulk  acoustic 

> 10 

Impedance  matching 

0.  AND  IKE., MULTIPATH  PftQfllXM 

Application  of  a memory  element  in  a generalized  system  with  coded  sig- 
nals is  shown  in  Figure  10.  A coded  signal  is  generated  and  put  on  a trans- 


Figure  16.  General  System  Using  Coded  Signals  and  Memory 
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The*  function  of  the  memory  element  is  to  sample  the  coded  signal,  store  It, 
and  return  it  to  the  receiver  with  sufficiently  low  distortion  so  that  the  output 
signal  from  the  matched  filtor  is  not  appreciably  distorted,  The  distortion 
allowable  on  the  output  signal  is  a system  requirement  which  will  vary  from 
system  to  system.  In  this  section  the  effect  of  the  acoustic  impedance  mis* 
match  between  transducer  and  delay  lino  is  considered  as  a source  of  distor- 
tion, and  a simple  relationship  Is  derived  between  the  impedance  mismatch 
and  the  distortion  of  the  output  signal  for  arbitrary  coding.  If  the  allowable 
distortion  la  known  from  system  requirements,  the  required  acouracy  of  acoustic 
matching  can  be  calculated. 
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Flguro  17,  Transducer  and  Delay  Medium  Schematic 


Whatever  the  coding  of  a signal  waveform,  it  has  a finite  bandwidth  and 
a oenter  frequency  fn.  The  delay  line,  shown  schematically  in  Figure  17,  is 
designed  to  have  a bandpass  characteristic  appropriate  to  the  signals  this  im- 
plies that  the  transducer  thickness  would  be  approximately  1/2  wavelength 
thick  at  fn.  Thus, 

d»v,/2f0  (20) 


where  d is  the  transducer  thickness  and  vt  is  the  velocity  of  sound  in  the  trans- 
ducer material. 


The  delay  line  functions  ns  a memory  element  by  conversion  of  the  elec- 
tronic signal  into  a sound  wave  by  the  piezoelectric  transducer.  The  sound 
pulse  travels  back  and  forth  through  the  delay  line  many  times.  Each  time 
that  it  reaches  the  transducer,  a small  fraction  of  the  sound  pulse  is  recon- 
verted into  electrical  energy  by  the  transducer  to  provide  a replica  of  the  orig- 
inal signal.  If  there  Is  an  acoustic  impedance  mismatch  between  the  transducer, 
delay  medium,  top  electrode,  or  the  bond,  a partial  reflection  occurs  at  that 
interface,  as  Indicated  in  Figure  17.  Because  the  top  electrode  and  the  bond  a 
are  thin,  the  most  important  Impedance  match  1h  that  between  the  transducer 
and  the  delay  media.  The  portion  of  the  signal  which  is  transmitted  Into  the 
transducer  is  delayed  with  respect  to  the  directly  reflected  portion  by  a time,  T: 


HO 


T « 2u/vt  ■ l/f0 


(27) 


I 

li 


j. 


which  la  juat  the  period  of  the  signal  center  frequency.  When  an  aoouetlo  pulse 
reachea  a mismatched  transducer  an  Infinite  series  of  pulses  is  generated, 
eaeh  succeeding  pulse  delayed  In  time  by  an  amount  T.  The  first  pulse  is 
the  reflected  one.  The  second  pulse  traverses  the  transducer  twice  and  Is 
delayed  by  time  T.  The  n'h  pulse  has  bounced  baok  and  forth  within  the  trans- 
ducer (n-1)  times.  If  the  acoustic  Impedances  are  nearly  matohed,  then  the 
second  pulse  Is  much  larger  than  all  the  others  because  the  reflection  coef- 
ficient at  the  transducer-delay  medium  interface  Is  very  small.  If  the  elec- 
trical output  from  the  delay  line  is  fed  into  a matched  filter,  the  filter  output 
will  be  an  Infinite  series  of  the  outputs  for  the  original  signal,  eaoh  output  de- 
layed by  time  T.  For  sufficiently  aocurate  impedance  matching  only  a single 
output  from  the  matched  filter  will  be  appreciable.  Depending  upon  system 
requirements,  it  may  be  possible  to  allow  the  output  from  the  matched  filter 
to  bo  broadened  by  several  factors  of  T,  and  the  Impedanoe  match  required 
can  be  relaxed  somewhat. 

The  acoustic  Impedances  of  the  bond  and  the  top  electrode,  and  their  at- 
tenuation coefficients,  are  also  Important.  If  the  bond  has  a thickness  that 
in  an  appreciable  fraction  of  an  acouatlc  wavelength,  It  will  act  aa  an  Imped- 
ance transformer  between  delay  medium  and  transducer  (Ref.  18)  end  will 
also  contribute  slgnlfioantly  to  the  delay  time,  T.  Theae  effects  oan  be  readily 
incorporated  into  the  calculation  whioh  follows,  but  they  hsve  been  omitted  for 
simplicity  and  clarity,  and  to  focua  attention  on  the  dominant  fabrication  con- 
straint! the  impedance  match  of  the  transduoer  to  the  delay  medium, 

The  calculation  of  signal  distortion  by  multiple  reflection  is  carried  out 
in  the  manner  of  the  discussion  given  above.  The  distorted  waveform  from 
the  delay  lino  is  not  calculated  directly,  because  its  form  will  depend  on  the 
particular  manner  of  coding  as  well  as  on  the  Impedanoe  mismatch.  Rather, 
the  signal  amplitude  is  calculated  for  eaoh  possible  time  delay  (time  delay 
■ 0,  T,  2T,  3T  ..,).  These  amplitudes  and  delays  give  the  response  of  the 
matched  filter  to  the  stored  signal,  since  the  response  to  the  initial  signal  Is 
known,  The  discussion  above  showed  that  a single  acoustic  pulse  became  un 
infinite  series  of  pulses  after  one  Impingement  on  the  transducer.  On  the 
next  encounter  with  the  transducer,  each  of  these  pulses  becomes  an  Infinite 
series.  In  Figure  18  the  situation  is  Illustrated  for  a sound  pulBe  which  has 
tmveled  through  the  delay  line  three  times  and  has  encountered  the  transduoer 
three  times.  The  quantity  N gives  the  number  of  times  that  the  pulse  has 
traveled  through  the  delay  line  and  Impinged  on  the  transducer.  The  quantity 
M denotes  the  time  delay,  M x T,  for  a group  of  pulses  with  the  same  total 
delay.  The  quantity  L accounts  for  the  number  of  reflections  at  the  transducer- 
delay  media  interface. 

For  a given  time  delay,  or  value  of  M,  there  Is  more  than  one  way  of 
generating  that  delay.  The  first  received  pulse  of  the  infinite  series  arrives 
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Figure  IB.  Pulses  Derived  from  n Single  Signal  After  Three 
Trips  Through  the  Delay  Line 

with  no  delay,  M * 0,  and  corresponds  to  three  reflections  from  the  Interface, 
as  shown.  The  next  pulHc  Is  delayed  by  T,  M « 1,  and  onrresponda  to  two 
reflections  from  the  interface  and  one  transmission  through  the  transducer. 

This  can  occur  in  throe  waysi  on  the  first,  second,  or  third  pass  through  the 
delay  lino.  Its  multiplicity  Is  3.  The  minus  sign  shows  that  the  transducer  - 
air  inter  face  has  refloated  the  wave  onoe,  causing  a 180°  phase  shift.  The 
third  pulse  has  a delay  of  3T,  TVI  « 3,  and  corresponds  to  two  passages  through 
the  transducer.  Note  that  this  can  happen  in  two  wayss  by  two  passages  through 
the  transducer  during  one  trip  through  the  drlny  line,  L * lj  or  by  one  passage 
through  the  transducer  on  each  of  two  trips  through  the  delay  line,  L « 2,  Each 
of  those  terms  has  a multiplicity  of  3. 

At  a time  delay  of  IV!  " N " 3 one  encounters  what  may  be  colled  the  main 
term}  hore  impedances  are  nearly  matched  so  that  the  reflection  coefficient  Is 
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small . The  main  term  lstheLBM»N*3  term,  which  consists  of  three 
passages  through  the  transducer,  one  on  each  trip  through  the  delay  line. 

There  Is  only  one  such  term.  The  two  other  terms  at  a time  delay  of  3T  cor- 
respond to  three  passages  through  the  transducer  for  one  trip  through  the  de- 
lay line,  L ■ U or  for  two  passages  through  the  transducer  during  one  trip  and 
one  passage  during  another  trip.  The  terms  for  successively  longer  delays  in- 
volve more  passages  of  the  sound  pulse  back  and  forth  within  the  transducer. 
For  M • 4,  time  delay  of  4T,  the  L ■ 2 case  shows  two  different  ways  of  getting 
the  same  time  delay  with  the  same  number  of  reflections  from  the  Interface. 
For  large  numbers  of  trips  through  the  delay  line,  N large,  or  for  long  time 
delays,  large  M,  the  drawing  out  of  all  diagrams  and  counting  becomes  tedious. 
The  problem  is  considered  mathematically  below. 


The  tranamission  coefficient  DtJ  is  defined  (Ref.  19)  as  the  ratio  of  the 
amplitude  of  the  transmitted  displaooment  wave  in  medium  j to  the  amplitude 
of  the  incident  ditplacement  wave  In  medium  1.  The  reflection  coefficient 
Vn  is  defined  as  the  ratio  of  the  amplitude  of  the  reflected  wave  in  medium  j 
to  that  of  the  Incident  wave  in  medium  j for  a wave  lncidient  on  the  boundary 
with  medium  i.  In  terms  of  the  impedances,  Z,  of  the  media 
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vt< 

Dji 


VZJ 

(28) 

Vzj 

-V„ 

(28) 

2Zi 

(30) 

zi+zj 

(31) 

z1+z. 

Cj  VJ 

(32) 

where  Oj  is  the  density.  One  notes  that  the  transmission  coefficients  for  the 
displacement  waves  can  be  greater  than  unity.  Energy  conservation  may  be 
verified  by  noting  that  the  energy  flux,  P,  is  given  by 

PJ  a 2 <ejNj)u,*s*  <33) 

where  S is  the  displacement.  The  subscripts  T and  D are  used  to  denote  the 
transducer  end  the  delay  medium.  A sound  pulse  incident  on  the  transducer 
may  be  written  as 

S„  (t,  z)  * j)a<tH  z/v)  -u(t-z/v+A  )]  S“  (t,  (34) 


where  the  function  S°0(t,z)  includes  the  form  of  coding,  s\ich  as  Equation  17 
or  21.  The  coordinates  are  chosen  with  z ■ 0 at  the  interface  between  the 
delay  medium  and  transducer,  and  the  expression  S\  (t,  z)  is  written  for  the 
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waves  leaving  this  interface  when  SB  from  Equation  34  is  incident  upon  it  once. 
(The  transducer  thickness  Is  d. ) : 

Sj  <t,  z)  « {[u(0)-u(A)]  S*  (0)VTB 

- fu(2d/vt)-u(2d/vt+A)]  Dre  DBT  SJ  (■2d/vt)ei^td 
+ [u(4d/vt)-u(4d/vt+A)l  DT9  V6T  DbT  8\  <-4d/vt)e*JM 


- [u(8d/vt)-u(6d/vt+ A)]  DTb  Vjr  DoTS*  (-8d/v»)e»jM 


(35) 


The  alternation  In  sign  accounts  for  the  190s  phase  shift  at  the  transducer-air 
interface.  A time  shift  operator  U(M)  is  defined  by 


and  an  operator  G(M) 


U(M)  f(t)  « f(t-2Md/vi) 
G(M)  ■ (-1)*VdmtU(M) 


From  Equation  35j 

SJ  - VtBS„+  DTbDbTS  <-l)Xr»U(M)Se 


M«1 


SJ  - VtDSB  + DTbD#t  t G(M)SB 

*■  1 


or  in  terms  of  a reflection  operator  R, 


(38) 

(37) 

(38) 

(39) 


SJ  • rsb 


r . vrB  + Eu2h  ! g(m) 

V-T 


(40) 


(41) 


8t 


Equation  40  expresses  the  fact  that  a single  incident  pulse  is  transformed  into 
an  infinite  series  of  pulses  when  it  strikes  the  delay  medium-transducer  inter- 
face. The  usefulness  of  Equation  40  arises  from  an  additive  property  of  the 
operator  G(N)  and  the  ability  to  account  for  N encounters  with  the  interface. 
Letting  Sjf  be  the  infinite  train  of  pulses  after  N encounters 

(42) 


S£  = r*sd 


and  applying  the  binomial  theorem 


R 


N1 


i=o  (N-L)iLi 


DTPDgT  * G(M) 


bt  m=i 


(4  3) 
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The  operator  G(M)  has  the  useful  property 

GOvl^GOvg  = GtMj-t-  , (44) 

which  can  be  used  to  reduce  the  product  of  L infinite  summations  to  a single 
summation 


n L f.  G(M)  ■ r PjiG(M)  (45) 

HOl 

where  P„  Is  the  number  of  ways  the  number  M can  be  constructed  by  adding  to- 
gether L positive  integers  not  including  zero.  The  P„  are  binomial  coefficients 
(Ref.  20) 

(M  -1)1  (46) 

P"  " (L-l)l  (M-L)  1 

which  gives  upon  substitution  into  Equations  43  and  45: 


RN 


V?„  + 


L<=1 


r 

L 


N1 


(N-L)  1 LI 


(M-l)l  / * vm -tTN—L—t.  r\L  \ i rv-L  (47) 
(L-l)  1 (M-Ul  ( VtpdtdDotVdt 


It  is  noted  that  VoT  - -VT0  and  that  the  summations  may  be  interchanged  if  Nl 
is  denoted  as  the  smaller  of  M or  Nj 


R".  V’"»  + ' WM)  L tS-tTl  M iL-lCT-U  <'1>  MB) 


Equation  4B  can  be  readily  interpreted  in  terms  of  the  diagrams  of  Fig- 
ure 7,  for  the  case  of  N B 3 trips  through  the  delay  line.  Since  U(M)  is  the 
operator  which  shifts  time  be  a factor  of  2M/Vt,  a given  value  of  M describes 
the  collection  of  terms  for  a time  shift  of  2M/Vt  and  that  collection  of  terms  is 
given  by  the  finite  summation  within  the  curly  brackets.  The  binomial  co- 
efficients give  the  correct  number  of  ways  that  reflections  and  transmissions 
can  occur.  The  quantity  (-1)L  accounts  for  phase  change  on  reflection  at  the 
transducer-air  interface.  The  attenuation  and  transmission  coefficients 
properly  account  for  changes  in  amplitude  on  transmission  and  reflection. 

The  advantage  of  Equation  48  is  that  is  may  be  readily  evaluated  numeri- 
cally for  hundreds  of  trips  through  the  delay  line.  Further,  the  results  of 
Equation  48  are  directly  applicable  to  the  output  of  a matched  filter.  If  the 
exact  form  of  the  initial  signal  is  specified,  the  results  of  Equation  48  can  be 
used  to  construct  the  output  waveform  after  N trips  through  the  delay  line. 

The  terms  of  Equation  48  are  written  out  below  for  N = 3,  three  trips  through 
the  delay  line,  and  for  the  first  four  pulses  of  the  infinite  series. 
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N •-  3 


M = 0 
L “ 1 

M-l  -3VTWI>6T 
M - 2 +3VfiTV?6nToDeT 
M = 3 “3V*T  V‘?nD-ienBT 

M“4  *8V“tV1bDt#DdT 


2 3 

+ 3VlBDf„DjT 

-flV#TVT#Df0Dft  -D?dD»t 
+OV*TVfBDfcDB\  + 3V6T  D^b  DbT 


Comparison  oi  these  expressions  with  Figure  18  shows  the  one-to-one  cor- 
respondence between  the  diagrams  and  the  terms  of  Equation  48. 


,/hen  impedances  arc  nearly  matched,  VftT  “ "^TB  1 and  Dtp  and  r)oT  ^ 1 . 
The  term  for  which  L = M = N has  been  called  the  main  term  of  the  nearly 
matched  case  because  it  is  the  pulse  which  has  been  transmitted  through  the 
transducer  once  for  each  passage  through  the  delay  line.  It  is  the  only  term 
completely  independent  of  VTo  and  VeT,  having  the  form  D?0DbT,  For  small  N 
it  will  be  the  dominant  term  in  the  series.  However,  as  N becomes  very  large 
the  multiplicative  factorial  expressions  become  large,  because  they  account 
for  the  number  of  ways  the  other  terms  can  occur.  The  ratio  of  the  next  largest 
term  in  the  sequence  of  terms  for  M » N to  the  main  term  is 


NVt„ 

DTbDbT 


These  ratios  suggest  that  if 


nvT0 


<1 


(49) 


(50) 


then  the  main  term  will  dominate  others  in  the  train  of  pulses.  The  impedance 
mismatch,  I,  is  twice  the  reflection  coefficient: 

T „ 2 I Ze~Zil  - 7 VT  (fill 

1 z8’'+zy  TD  ' ' 

A useful  criterion  for  dominance  by  the  main  term  is  that 


NI  = 1 (52) 

or  in  words,  a 1 -percent  impedance  mismatch  results  in  one  dominant  term  in 
the  pulse  train  and  little  distortion  for  one  hundred  or  less  trips  through  the  de- 
lay line.  In  Figure  19,  results  are  plotted  for  the  case  of  N = 100  trips  with  a 
1 -percent  impedance  mismatch.  The  horizontal  time  axis  is  in  units  of  the  delay 
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time,  T;  because  oT  transit  through  tfie  transducer,  T - 2tl/VT.  Time  is  taken 
as  zero  for  the  main  term,  L a M «*  N,  The  Hum  of  all  terms  for  a given  time 
delay,  M value,  Is  plotted  vertically  ns  an  amplitude  for  that  time  delay.  The 
incident  signed  is  assumed  to  have  unity  amplitude,  A verbal  Interpretation  of 
Figure  10  follows  in  the  next  paragraph. 


N * 100  TRIPS 

1%  IMPEDANCE 
MISMATCH,  t >0.01 


TIME  OF  OUTPUT  FROM  THL  DELAY  LINE  IN  UNITS  OF 
DELAY  TIME , T , OF  THE  TRANSDUCER 


1!),  Amplitude  Versus  Time  of  Output  for  100  Trips  Through 
the  Delay  l.lm*.  Calculated  for  l-Peroent  Mismatch  of 
Impedance  lletween  Transducer  nnd  Delny  Medium 


A coded  signal  of  unity  amplitude  is  stored  in  a delay  line  memory  device 
with  a 1 -percent  impedance  mismatch  between  the  delay  medium  and  the  trans- 
ducer. Because  of  reflections  at  the  interface  of  the  transducer  and  delay 
medium,  the  single  signal  becomes  an  infinite  train  of  duplicate  signals,  sep- 
arated from  one  another  by  the  transit  time  T through  the  transducer.  The 
largest  of  these  signals  has  an  amplitude  of  0.76.  It  is  preceded  and  followed 
at  a time  interval  T by  signals  of  amplitude  0.44  and  at  a time  interval  of  2T 
by  signals  of  amplitude  0.1;  at  a time  interval  3T  it  is  preceded  by  a signal 
of  amplitude  0.019  and  followed  by  a signal  of  amplitude  0.  020;  and  so  on. 

Note  that  only  the  points  at  integral  multiples  of  T can  have  any  meaning.  In 
subsequent  graphs  the  points  are  joined  with  continuous  curves  to  aid  the  eye 
in  following  from  point  to  point,  but  only  the  points  have  any  meaning. 

In  Figures  20  and  21,  amplitude  is  plotted  versus  delay  time,  T,  for  sig- 
nal trains  that  result  from  100  trips  through  the  delay  line  for  impedance  mis- 
matches of  0.1,  1,  2,  3,  4,  and  5 percent.  One  sees  that  for  0. 1 -percent 
and  1-percent  mismatch,  the  M = N (time  0)  signal  is  largest  and  the  other 
signals  decrease  in  amplitude  monotonically  with  time  difference.  For  a mis- 
match of  2 percent  or  greater,  the  M = N signal  is  not  the  largest  and  the  am- 
plitudes are  appreciable,  yet  variable  in  magnitude,  for  an  appreciable  number 
of  time  delays.  The  criterion  of  Equation  52  which  states  NI  « 1 for  dominance 
of  the  M = N signal  is  verified.  In  addition,  one  can  define  a useful  width  par- 
ameter, W,  as  the  number  of  delay  times,  T,  encountered  before  the  signal 
level  drops  by  a factor  of  3 and  continues  to  decline  monotonically.  For  small 
values  of  mismatch  it  is  found  empirically  that 

W = 3IN  (53) 

Beyond  this  width  the  signals  in  the  train  drop  in  power  by  about  10  decibels  per 
delay  period  T . 

The  form  of  Equation  52,  NI  = 1,  suggests  that  the  shape  of  the  amplitude 
versus  T curves  should  be  dependent  principally  upon  NI.  As  one  varies  N and 
I,  keeping  the  product  fixed,  one  would  not  expect  the  details  of  amplitude  vari- 
ation near  T = 0 to  be  reproduced,  but  would  expect  the  regions  of  the  curves 
that  decrease  monotonically  to  reproduce.  This  is  shown. to  be  the  case  in  a 
series  of  calculations  for  particular  combinations  of  transducer  and  delay  line 
materials. 

The  Impedances  and  the  mismatch  of  impedances  for  LiNbOj  and  for  quartz 
are  given  in  the  matrix  of  Table  6.  Lithium  niobate  is  chosen  as  an  example 
because  of  the  available  piezoelectrics  it  most  closely  matches  the  impedances 
of  the  low -loss  delay  media  for  shear  modes.  Quartz  is  used  as  an  example 
because  of  its  availability,  frequent  use,  and  poor  impedance  match.  The  mis- 
match to  LiNbOr,  can  be  reduced  if  the  level  of  effort  entailed  is  warranted. 

The  impedance  can  he  lowered  somewhat  by  carefully  exploring  rotated  cuts, 
with  due  regard  to  the  effective  coupling  coefficient  and  the  generation  of  spu- 
rious modes.  The  impedance  can  be  raised  by  growing  mixed  crystals  withLiTaO,. 
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TIME  OF  OUTPUT  FROM  DELAY  LINE  IN  UNITS  OF  DELAY  TIME, 
T,  OF  THE  TRANSDUCER 


Figure  20.  Amplitude  Versus  Time  of  Output  for  100  Trips 

with  0.1,  1,  2,  and  5 Percent  impedance’  Mismatch 


AMPLITUDE 


In  Table  6 the  acoustic  impedances  were  taken  from  Tables  1 and  3.  For  YHi15 
the  Impedance  value  used  was  20,2  x 10*  gm/sec  cm8. 

Figure  22  shows  a plot  of  amplitude  versus  time  for  LiNbOn  on  MgAl,04. 
The  number  of  pulses  has  been  chosen  to  provide  NI  products  of  1,  2,  and  3 
for  comparison  with  Figures  20  and  21.  One  finds  close  similarity  in  the  shape 
of  the  curves  for  constant  NI  product,  verifying  the  usefulness  of  Equations  52 
and  53.  Figure  23  shows  similar  results  for  LiNbO,  on  Y3Al,01#i  Note  that 

Table  6 

PERCENT  IMPEDANCE  MISMATCH 
OF  TRANSDUCERS  AND  DELAY  LINES 


I.INI>0„  X -mil 

KIOt,  Y-CUt 

MkAI.,11, 

'1.  4 t 

7I>< 

•>  ,M.t  i, 

1 . 'A  i 

711  i 

1 1 i 

(I7< 

IN  UNITS  OF  T 

Figure  22.  Amplitude  Versus  Time  of  Output  for  LiNbO,  Transducer 
on  MgAljO^  for  3,  30,  f>0,  and  00  Trips 
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AMPLITUDE 


TIME  OF  OUTPUT  FROM  DELAY  LINE 
IN  UNITS  OF  T 


Figure  23.  Amplitude  Versus  Time  of  Output  for  LiNbO,,  on  YAOS 
for  75,  150,  and  300  Trips 
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the  close  Impedance  match  provides  N = 75  trips  before  the  condition  NI  = 1 
Is  reached.  In  Figure  24  a similar  plot  1b  given  for  LiNbOs  on  In  this 

case  the  NI  ■ 1 condition  is  reached  for  only  nine  trips  of  the  signal  through 
the  delay  line.  In  Figure  25  the  cases  of  2,  3,  and  10  trips  are  shown  for 
SiOa  on  MgAla04,  In  spite  of  the  good  insertion  loss  and  bandpass  characterls 
tics  of  this  combination,  the  poor  impedance  matoh  makes  it  a poor  choioe  for 
memory  of  coded  signals  for  any  number  of  trips  through  the  delay  line. 


TIME  OF  OUTPUT  FROM  DELAY  LINE 
IN  UNITS  OF  T 


Figure  24.  Amplitude  Versus  Time  of  Output  for  I.iNbO,  on  Vllw 
for  0,  1H,  and  36  Trips 
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TIME  OF  OUTPUT  FROM  DELAY  LINE  IN  UNITS  OF  T 


Figure  25. 


Amplitude  Versus  Time  of  Output  for  SiO„  on 
for  2,  3,  and  10  Trips 


MrA1,0( 


Summary  and  Conclusion 


Impedance  mulching  of  transducer  and  doluy  media  is  Important  In  order  to  min- 
imize distortion  of  stored  coded  signals.  An  expression  has  been  provided  in 
Equation  48  i’or  calculating  the  effects  of  the  impedance  mismatch  in  detail. 

Two  simple  criteria,  Equations  52  and  53,  arc  given  for  estimating  the  de- 
gree of  impedance  matching  required  to  meet  n system  specification.  The 
allowable  distortion  on  the  output  of  the  matched  filter  must  be  specified  for 
the  system  of  Interest  if  the  degree  of  impedance  matching  is  to  be  calculated 
for  the  memory  element. 

The  constraints  of  impodance  matching  do  not  create  problems  in  the  in- 
sertion loss  characteristics  of  the  delay  lines.  The  diameter  of  the  trans- 
ducers must  be  kept  largo  in  order  to  minimize  diffraction  losses  and  provide 
long  delay  times.  The  large  area  and  high  dielectric  constant  result  in  a very 
low  electrical  impedance  for  the  transducer.  In  addition,  the  conversion  ef- 
ficiency of  LiNbO,  is  too  high  for  long  memory  times.  Figure  26  shows  cal- 
culated values  of  Insertion  loss  versus  frequency  for  a delay  line  of  MgAl#04 


with  Source  Impedance  and  Series  C'npaeltance  Varied 


H4 


with  an  X-cut  LiNbOa  transducer  centered  at  200  MH*  bonded  to  it,  with  1 pm 
of  gold  and  0,1 -pm -thick  gold  top  electrode.  'When  driven  from  a low  imped- 
ance, 0.8-ohm  source  this  combination  has  a very  low  insertion  loss,  less 
than  1 decibel,  end  varies  in  insertion  loss  by  leas  than  4 decibels  over  a 
100-MH*  bandwidth.  When  driven  from  a 50-ohm  source  the  insertion  loss 
variation  with  frequency  does  not  change  appreciably,  but  the  minimum  in- 
sertion loss  is  raised  substantially,  to  30  decibels. 

Further  control  of  the  insertion  loss  can  be  attained  by  driving  the  trans- 
ducer via  a aeries  capacitor.  With  6 pfd  of  aeries  capacitance  the  insertion 
loss  is  increased  to  about  50  decibels,  without  appreciable  changes  in  the  fre- 
quency dependence  other  than  a alight  increase  in  the  frequency  of  minimum 
! insertion  loss.  The  use  of  series  capacitance  has  been  reported  previously 

! by  Larson  and  Winslow  (Ref.  2),  who  showed  that  a marked  resonant  effect 

on  bandpass  characteristics  resulted.  Their  study  involved  X-cut  LlNbOB 
i bonded  to  rutile  with  indium.  Unlike  the  present  case,  the  transducer,  bond, 

and  delay  line  were  badly  matched  acoustically.  Resonance  was  observed  when 
i the  delay  line  was  mismatched  both  electrically  Bnd  acoustically.  The  acoustic 

j impedance  matching  required  for  memory  of  coded  signals  eliminates  such 

| resonant  effects.  Insertion  loss  characteristics  maybe  tailored  to  needs  by 

j controlling  the  source  impedance  and  by  the  addition  of  a series  capacitance. 
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Section  3 


DELAY  LINE  FABRICATION 


Fabrication  of  delay  rods  and  transducer,  deposition  of  metals  for  bonds 
and  electrodes , bonding  of  the  transducer  to  the  delay  line,  and  thinning  of  the 
bonded  transducers  have  been  explored  experimentally.  The  selection  of  metals 
as  bond  and  electrode  materials  is  considered  in  this  section.  The  relevant 
phenomena  are  indicated,  a brief  literature  review  and  status  of  knowledge  is 
given,  and  a practical  system  is  described  for  simultaneously  minimizing 
attenuation  in  the  bond  and  matching  Impedances.  A description  is  given  of 
the  equipment  used  and  the  general  methods  employed. 

A.  MATERIALS  FOB  METAL  BONOS 

There  are  several  criteria  to  apply  to  metals  which  are  to  be  UBed  for 
bonds j 


• Ease  of  thcrmocomprcssion  bonding 

• Matching  acoustic  impedance 

• Minimizing  Hound  attenuation 

Ease  of  forming  a bond  has  been  found  (Ref.  21)  to  correlate  with  yield 
stress  of  the  material.  Thus  the  ductility  of  the  crystal  lattice  is  an  important 
factor,  and  this  decreases  as  one  goes  through  the  sequence  FCC:  HCC ; HCP. 
This  criterion  points  to  the  use  of  inetals  with  the  FCC  copper  structure. 

The  acoustic  impedances  of  several  of  the  materials  of  interest  are  listed 
in  Table  7.  The  delay  media  and  the  transducer  material  have  impedances  in 
the  range  of  20.2  to  23.4  xlOr  g cm/sec.  The  metals  listed  have  the  desired 

Table  7 

SHEAR  MODE  IMPEDANCES 


■M*. 



Zno 

1 — 

Material 

z 

^ m 

^100 

T's* 

Gold 

24.0 

22.4 

24. 8 

20.7 

1H.  (i 

Pall  ndium 

24 . 0 

22.il 

2fi.fi 

20.3 

in.  l 

Silver 

17.7 

- “ 

“ ~ 

MgAlpq 

. ~ 

“ - 

~ - 

23.4 

■ •" 

YAG 

-- 

22.il 

-- 

UNIjO, 

-- 

-- 

-- 

22.  r, 

^ ^68 

-- 

-- 

20.2 

— 

tlnita  of  Impedance  vuIuih  10e  gin/cm"  sec. 


FCC  copper  wtructurc  and  form  solid  solutions  over  a large  composition 
range.  The  columns  labeled  Zm,  ZLl0,  and  Z100  contain  impedances  for  shear 
waves  propagating  in  the  [111],  [110],  and  [100]  crystallographic  directions, 
respectively.  The  column  headed  ”2  gives  the  value  of  impedance  for  a poly- 
crystalline sample  with  random  orientation  of  grains.  From  Table  7 one  sees 
that  the  Impedances  of  all  of  the  oxide  crystals  can  be  met  with  alloys  of  the 
Pd  Au-Aff  system  if  the  alloy  has  a random  grain  structure. 

In  the  work  of  Larson  and  Winslow  (Ref.  2)  with  gold  films,  the  films  were 
found  to  have  a [111]  texture  and  an  impedance  of  22  percent  less  than  the  value 
for  bulk  gold.  The  exact  causes  of  the  discrepancy  are  not  understood  and  may 
involve  such  factors  as  the  extent  of  interdif fusion  during  ultrasonic  welding, 
the  density  of  the  deposited  metal,  and  the  details  of  extracting  the  velocity  and 
attenuation  of  sound  in  gold  from  the  insertion  loss  data.  The  effect  of  a 22- 
percent  impedance  reduction  Z is  shown  in  the  column  labeled  Zm.  One  would 
expect  a somewhat  reduced  impedance  in  comparing  a deposited  film  with  a 
single  crystal.  Whether  a 22 -percent  reduction  in  impedance  for  a thin  metal 
film  is  realistic  and  a limit  to  what  can  be  achieved  is  under  investigation. 

Even  if  this  is  the  case,  the  Pd-Au  system  offers  several  advantages. 

The  deposition  of  an  alloy  makes  it  easier  to  obtain  a film  with  random 
orientatio_n  of  grains,  with  an  increase  in  impedance  as  a result  of  a random 
texture,  Z,  rather  than  a Till]  film  texture,  Zm.  Because  alloys  of  gold 
and  palladium  follow  Vegard's  law  to  within  a few  percent  (Ref. 22),  density 
varies  nenrly  linearly  with  alloy  composition.  However,  the  torsion  modulus 
does  not  vary  linearly  with  composition  (Ref.  23):  as  a result,  the  acoustic  im- 
pedance of  alloys  is  greater  than  that  of  either  component,  as  shown  in  Table  8. 
Other  elements  may  be  added  to  the  Pd-Au  system.  Additions  of  silver  can  be 
made  over  wide  composition  ranges  to  reduce  impedance.  Over  restricted 
composition  ranges,  additions  of  platinum,  rhodium,  and  iridium  can  be  made, 
and  the  latter  should  be  useful  in  increasing  film  impedance. 


Table  8 

SHEAR  MODE  IMPEDANCES,  OF  Au-Pd  ALLOYS 


Atom  Percent  Palladium  \ 

Z 

0 

24. 0 

8.80 

24.0 

17.0 

20,0 

HI . (i 

20.0 

41.2 

27.2 

27.8 

lull 

24.0 

Units  of  impedance  values  10*  gm/cm8  sec. 
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The  third  criterion  for  the  selection  of  a metal  system  to  be  used  as  a 
bond  is  minimization  of  sound  attenuation  in  the  bond.  Attenuation  in  such 
metals  arises  from  four  causes:  1)  scattering  of  sound  by  the  crystal  grains, 

2)  dislocation  attenuation  of  the  sound  wave,  3)  attenuation  of  the  sound  wave 
by  electrons,  and  4)  attenuation  of  the  sound  wave  by  phonons.  It  is  shown  here 
that  these  effects  are  minimized  by  use  of  an  alloy.  Attenuation  by  scattering 
from  crystal  grains  is  small  for  the  metal  thicknesses  contemplated,  and  is 
minimized  by  deposition  procedures  which  can  readily  produce  films  of  metal 
with  grains  much  smaller  than  an  acoustic  wavelength.  Attenuation  by  phonons 
ia  not  an  appreciable  fraction  of  the  attenuation  in  the  delay  medium.  However, 
attenuation  due  to  dislocations  and  to  electrons  can  be  appreciable. 

Mason  and  Rosenberg  have  studied  attenuation  in  aluminum  (Ref.  24), 
cadmium  and  copper  (Ref.  25),  and  lead  (Ref.  26).  They  show  that  the  attenu- 
ation can  be  separated  into  a portion,  rt,  which  is  caused  by  electrons  and  varies 
quadratically  with  frequency,  and  a portion,  r„,  which  is  caused  by  dislocations 
and  follows  a normalized  loss  curve.  Typically,  dislocation  losses  dominate  at 
frequencies  up  to  a fewhundred  MHz,  and  electron  damping  dominates  in  the  GHz 
frequency  range.  The  electron  related  attenuation,  rt,  was  shown  to  be  pro- 
portional to 


TWc/pVa 


where 


u = angular  frequency 
o ■ electrical  conductivity 
p B density 
V ■ sound  velocity 

The  electrical  resistivity  of  gold -palladium  alloys  has  been  measured  (Refs. 

23-27),  and  values  (Ref.  23)  are  quoted  in  Tabic  9. 

Table  9 

CONDUCTIVITY  OF  AU-PD  ALLOYS  AT  0°C 


Atom  Percent  Palladium 


4 . fiii 
fi.Zl 
P.70 


The  extant  data  for  shear  wave  attenuation  in  gold  films  is  that  of  Larson 
and  Winslow  (Ref.  2)  for  films  of  [111]  texture.  They  find  a value  of  1000  dB/cm 
at  1 GHz,  and  at  that  frequency  the  attenuation  is  expected  to  be  dominated  by 
interactions  with  electrons.  From  Table  9,  then,  it  is  apparent  that  re  for 
Au-Pd  alloys  can  be  as  low  as  100  dB/cm.  This  factor  of  ten  lowering  of  at- 
tenuation reduces  to  a negligible  amount  the  contribution  of  the  bond  to  total 
delay  line  attenuation.  The  contribution  of  a gold  layer  to  total  attenuation  can 
be  comparable  to  that  of  a low -loss  delay  medium  such  as  MgAl#04. 

The  attenuation  due  to  dislocations,  rB,  can  be  appreciable.  The  available 
experimental  results  (Refs.  24-26)  separate  these  losses  from  electron  damping 
and  provide  values  for  the  dislocation  drag  coefficients  in  the  metals  studied. 
They  do  not  provide  useful  means  for  predicting  dislocation  damping  in  other 
materials.  However,  Cain  and  Thompson  (Ref.  28)  have  studied  sound  velocity 
in  Cu-Al  alloys.  In  that  study  they  note  that  the  dislocation  attenuation  in  the 
alloys  is  substantially  lower  than  in  pure  copper,  as  a result  of  pinning  of  dis- 
locations. This  result  is  expected  from  known  properties  of  dislocations.  One 
expects  reduction  in  rB  by  orders  of  magnitude  in  going  from  pure  gold  to  Pd-Au, 
Data  are  not  available  for  making  a quantitative  estimate  of  the  effect  of  alloying 
on  rD. 

The  advantages  of  the  Au-Pd-Ag  system  for  acoustic  bonds  are  summarized 
as  follows.  This  system  permits  acoustic  impedance  matching  to  the  Iow-Iobs 
crystals  of  interest.  Metal  films  can  be  deposited  with  fine,  randomly  oriented 
grains.  Attenuation  due  to  electrons  can  be  reduced  to  a tenth  that  of  gold,  and 
attenuation  due  to  dislocations  can  be  markedly  reduced.  Quantitative  estimates 
of  the  latter  are  not  available. 

B.  FABRICATION  EQUIPMENT 

Metal  deposition  is  carried  out  by  evaporation  or  sputtering.  The  appar- 
atus used  for  diode  or  triode  sputtering,  shown  in  Figure  27,  consists  of  two 
stainless  steel  stations  and  control  circuitry.  The  vacuum  stations  use  copper 
gaskets  and  graded  metal-to-glass  seals  exclusively.  Each  station  is  pumped 
by  a mercury  diffusion  pump,  which  is  isolated  from  the  station  by  two  glass 
liquid  nitrogen  traps.  The  mechanical  roughing  and  backing  pump  is  isolated 
from  both  the  diffusion  pump  and  the  system  by  a glass  liquid  nitrogen  trap 
filled  with  copper.  The  stations  routinely  reach  10“T  torr  prior  to  backfilling 
with  gas  for  sputtering. 


Tne  bonding  of  a transducer  to  a delay  rod  has  been  carried  out  in  a number 
of  ways,  such  as  the  use  of  low -melting -point  solders  (Ref. 29),  cured  varnishes 
(Ref.  30)  or  epoxies  (Ref.  31),  thermocompression  welding  (Ref.  2),  and  ultra- 
sonic welding  (Ref.  2).  These  methods  have  been  reviewed  and  extended  by 
Larson  and  Winslow  (Ref.  2).  Electric  field  assisted  bonding  (Refs.  32  and  33) 
may  be  applicable  but  has  not  been  explored.  Preliminary  experiments  with 


60 


Figure  27.  Sputtering  System,  D-C  Diode  or  Triode  Operation 

ultrasonic  bonding  indicated  that  this  method  would  not  be  practical  for  large  - 
diameter  transducers  unless  thick  metal  films,  an  appreciable  fraction  of  an 
acoustic  wavelength  in  thickness,  were  used.  With  thin  films,  excessive 
pressures  were  required  to  form  a bond  by  the  ultrasonic  technique.  Thormo- 
compression  bonding  also  requires  high  pressures  (Rof,  34)  when  performed 
in  air.  When  carried  out  In  vacuum,  substantially  lower  pressures  are  ade- 
quate (Ref.  21). 

Therefore  a thermocompression  bonding  system  has  been  built  which  op- 
erates inside  a vacuum  system.  In  this  system  the  top  surface  of  the  delay 
line  rod  and  what  will  eventually  be  the  bottom  surface  of  the  lithium  niobate 
transducer  are  held  in  a jig,  in  a coplanar  position,  underneath  a metal  evap- 
oration source.  The  jig  is  held  in  a motor-driven  rotating  dUk.  The  motor 
is  shown  in  the  upper  portion  of  Figure  28.  The  jig  is  shown  as  two  square 
holes  just  to  the  left  of  center  in  Figure  28;  one  square  hole  holds  the  rod,  the 
other  holds  the  transducer. 

In  the  first  position  of  the  disk,  a thin  chromium  film  is  evaporated  onto 
both  the  rod  and  the  transducer.  The  disk  is  then  rotated  to  a second  position 
in  the  bell  jar,  where  a second  thin  metallic  film,  perhaps  silver,  is  evaporated 
onto  both  the  rod  and  the  transducer.  In  the  third  position  the  transducer  plnte 
is  inverted  by  means  of  a hinge  built  into  the  jig  so  that  the  metal- coated  face 
of  the  transducer  now  lies  on  top  of  the  delay  line  rod.  It  is  also  accurately 
centered  over  the  rod. 
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Figure  28.  Substrate  Carrier  for  Depositing  Metals  and  Bonding 
Transducers  in  Vacuum 


In  the  fourth  position  the  transducer- metal  film- rod  sandwich  is  placed 
between  the  jaws  of  a vise  which  is  operated  by  an  hydraulic  ram.  This  ram 
is  actuated  from  outside  the  vacuum  chamber  by  a piston  pump  which  pres- 
surizes the  oil  inside  a stainless  steel  line  leading  into  the  vacuum  chamber 
and  then  to  the  ram. 


A collection  of  power  feedthroughs,  hydraulic  cylinder,  ram,  and  so  forth 
for  use  with  the  vacuum  system  is  shown  in  Figure  29.  Experiments  have 
shown  that  close  alignment  of  the  transducer,  delay  rod,  and  ram  surfaces 
are  essential  if  forces  are  to  be  evenly  distributed  so  that  the  transducer  will 
not  fracture.  This  alignment  requires  both  a rotation  and  translation  of  the 
ends  of  the  delay  rod.  Thin  shims  of  teflon  have  been  used  successfully  to 
provide  the  required  motion. 


This  present  bonding  is  all  done  at  room  temperature.  Once  the  sand- 
wich has  been  pressurized  and  the  metal  bond  is  secure,  the  pressure  is  re- 
leased and  the  transducer  is  ready  for  the  next  step.  This  step  is  thinning  of 
the  transducer  so  that  it  will  resonate  at  the  appropriate  frequency.  This 
thinning  must  be  done  after  the  bonding,  because  a transducer  of  the  resonant 
thickness  would  crack  easily  under  the  bonding  pressure.  The  vacuum  sta- 
tion, as  shown  in  Figure  30,  is  equipped  for  RF  sputtering  of  insulators,  and 
it  can  be  used  in  this  mode  for  thinning  of  bonded  transducers.  The  bonded 
transducer  can  also  be  thinned  by  a mechanical  process  of  grinding  and  polishing. 
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Figure  29,  Parts  for  Vacuum- Bonding  Station 


Figure  30.  Vacuum  Station  for  Bonding  or  RF  Sputtering 
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A mechanical  system  For  thinning  transducers  ia  being  developed  using 
techniques  similar  to  those  of  Twyman  (Ref.  35).  In  this  apparatus,  however, 
metal  laps  replace  the  pitch  laps  of  Twyman.  The  machine  is  built  from  a 
precision  computer  disk  memory  drive  unit  and  is  shown  in  Figure  31.  Dia- 
mond is  embedded  in  the  metal  laps  and  a rotating  sponge  collects  loose  par- 
ticles from  the  lap.  Figure  31  also  shows  the  box  which  houses  the  laps  for 
different  grades  of  diamond,  one  of  the  large  polishing  blocks,  and  a large- 
diameter  optica1,  flat  used  to  simultaneously  evaluate  the  condition  of  the  lap 
and  the  sample.  Aluminum  laps  have  proven  to  be  satisfactory  for  diamond 
of  3-pm  size  and  above,  but  they  induce  scratching  when  used  with  finer  grades 
of  diamond.  Copper  laps  have  been  prepared  for  use  with  1-pm  and  1/2-pm 
grades  of  diamond. 


Figure  31.  Polishing  Machine 

C.  UNDOPED  SPINEL  CRYSTALS  WITH  LITHIUM  NjOjATE  TRANSDUCERS 


The  most  promising  material  to  date  for  fabrication  of  low -loss  acoustic 
delay  lines  at  a frequency  of  200  to  300  MHz  is  spinel  (MgAl804),  The  design 
of  the  transducer -bond -crystal  sandwich  for  MgAl#04  using  a LiNbO,  piezoelec- 
tric transducer  has  been  described  in  Section  2 ("Delay  Dine  Design").  A 
nearly  optimum  design  is  the  use  of  an  X-cut  LiNbO,  crystal  with  a thickness 
of  1/2  wavelength  to  generate  shear  waves  at  200  MHz;  this  will  be  12  microns 
thick.  It  is  to  be  bonded  to  one  end  of  a spinel  crystal  rod,  which  has  a [100] 
axis  along  the  length  of  the  rod,  by  the  use  of  a pure  silver  or  gold  bond.  The 
metal  bond  thickness  should  be  considerably  less  than  1/2  wavelength  thick; 
this  means  preferably  of  the  order  of  0.  5 micron  thick. 


Two  boules  of  Czochralski  grown  spinel  were  purchased  from  Crystal 
Technology,  of  Mountain  View,  California.  One  was  grown  with  a [111]  axis 
and  was  10  cm  long  and  3 cm  in  diameter.  The  other  was  grown  along  a [100] 
axis,  and  was  5 cm  long  and  3 cm  in  diameter.  These  boules  were  scanned 
with  a visible  light  beam  when  immersed  in  an  index -matching  fluid  in  order 
to  select  regions  with  a minimum  number  of  light  scattering  centers.  A total 
of  six  separate  crystal  rods  were  cut  from  the  boules,  two  from  the  [100] 
boule  and  four  from  the  [lll]boule.  The  rods  were  about  1 cm  in  diameter  and 
0.8  to  2.5  cm  long;  they  were  all  cut  with  a [100]  crystallographic  direction 
along  their  length.  These  rods  were  centerless -ground  and  then  sent  to  the 
Adolph  Meller  Company,  in  Providence,  Rhode  Island,  where  the  end  faces 
were  carefully  ground  and  polished  to  be  optically  flat  and  parallel. 

Shear  wave  transducers  of  X-cut  LiNbOa  0.95  cm  in  diameter  and  125  pm 
thick  were  purchased  from  the  Valpey-Fisher  Corporation,  in  Hopkinton, 
Massachusetts.  The  procedure  is  to  compress  ion -bond  these  transducers  to 
the  spinel  crystals  in  vacuum.,  using  evaporated  metal  films  on  the  two  op- 
posing faces.  Once  bonded,  the  exposed  face  of  the  transducer  is  to  be  ground 
and  polished  by  removing  about  113  pm  of  material,  thus  attaining  the  desired 
1/2  wavelength  thickness.  This  polished  face  will  then  be  covered  with  an  evap- 
orated metal  contact,  to  act  as  a second  electrode.  The  special  Twyman-type 
grinding  and  polishing  machine  was  used  to  remove  the  requisite  amount  of 
LiNb03  and  keep  the  polished  face  flat  and  parallel.  This  was  installed  in  a 
laminar -airflow  hood  in  order  to  keep  room  dust  out  oT  the  system. 

Some  grinding  experiments  on  fused  quartz  rods  and  plates  were  run  with 
this  machine.  The  grinding  was  done  with  6-pm-diametcr  diamond  powder 
embedded  into  the  surface  of  an  aluminum  wheel.  The  wheel  was  30  cm  in 
diameter  and  2.  5 cm  thick,  and  rotated  at  a speed  of  30  rpm.  At  a vertical 
pressure  of  2 x 10B  dyncs/cma  the  quartz  removal  rate  was  1.  1 pm/mln,  Finer 
grinding  was  carried  out  with  a similar  copper  wheel,  using  embedded  1-pm- 
diametcr  diamond  powder.  The  cutting  rate,  at  a pressure  of  2 x 10s  dynes/ 
cm3,  was  0.  12  pm/min.  Somewhat  similar  grinding  rates  can  be  anticipated 
for  LiNbOij.  The  grinding  rates  were  measured  using  comparison  standards 
of  metal  gage  blocks  and  a precision  dial  indicator  made  by  the  Federal  Pro- 
ducts Corporation  (Providence,  Rhode  Island),  The  rates  wore  found  to  be 
reproducible.  One  problem  encountered  in  the  reproducibility  waB  to  find  a 
mounting  wax  for  the  spinel  or  quartz  rods  that  would  not  cold-flow  during  the 
grinding  and  measuring  operation,  but  that  also  had  a low  enough  melting  point 
to  enable  the  rods  to  be  mounted  in  the  brass  polishing  jig.  The  best  was  found 
to  be  .(ewelor's  Special  Wax  made  by  the  Princeton  (Illinois)  Sealing  Wax  Company. 

D.  METAL  BONDING 

The  evaporator  for  the  metal  films  in  Figure  30  was  evaluated  with  fused 
quartz  test  blanks.  Chromium  and  then  silver  films  were  evaporated  onto  the 
quartz  irom  tungsten  wire  filaments  (Materials  Research  Corporation,  Orange- 
burg, New  York)  in  a vacuum  of  about  2 x 10""'  tnrr.  The  time  and  temperature 
were  empirically  adjusted  so  that  in  15  minute  s a 700?-lhick  layer  of  chromium 


was  deposited.  The  sample  was  then  rotated  into  the  second  section  of  the  cham- 
ber, where  4000  of  silver  wore  deposited  in  about  45  minutes  time.  The 
chromium  layer  thickness  was  measured,  on  a dummy  substrate  which  wds 
left  unsilvered,  by  means  of  a Dektak"*  profilometer.  The  silver  layer  thick- 
ness was  measured  by  weighing  a second  dummy  substrate,  before  and  after 
evaporation,  on  a sensitive  Mettler  balance  manufactured  by  Will  Scientific, 

Inc.,  Rochester,  New  York. 

After  two  fused  quartz  rods  were  coated  with  the  double  layer  of  metal 
film  they  were  rotated  to  the  third,  and  then  the  fourth,  section  of  the  vacuum 
chamber,  whore  they  were  automatically  positioned  between  the  jaws  of  the 
hydraulic  ram  (Figure  29).  The  maximum  pressure  that  could  be  tolerated 
by  the  0.  95-cmrdiameter  rods  was  about  2,2  x 10s  dynes/cm*  (32,  000  lbs/in8); 
at  this  pressure  they  invariably  shattered.  On  smaller- diameter  quartz  rods 
pressure  was  occasionally  run  up  to  3.  G x 10®  dynes/ cm8  before  shattering. 

As  a rule,  the  compression  bonding  of  quartz  was  attempted  at  about  2.0  x 
10®  dynes/cma.  Of  the  three  tries  that  were  made,  none  was  successful.  The 
metal  layers  pulled  away  from  the  quartz,  even  though  the  quartz  had  been 
very  carefully  cleaned  prior  to  the  evaporation.  The  reasons  for  the  lack  of 
success  are  unclear.  More  work  is  needed  to  perfect  this  bonding  technique. 


1M  Trademark  of  the  Sloan  Technology  Corporation,  Santa  Barbara,  California. 
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Section  4 


GROWTH  AND  CHARACTERIZATION  OF  YB  ^ CRYSTALS 
A BACKGROUND 


The  Research  and  Development  Center  of  the  General  Electric  Company 
has  supported  efforts  on  understanding  the  theory  of  sound  attenuation  and  on 
preparing,  evaluating,  and  measuring  the  properties  of  crystals  to  develop 
useful  microwave  memory  crystals  and  to  test  theory.  This  is  borne  out  in 
the  early  reported  work  of  Jacobsen  (Ref.  36),  Woodruff  and  Ehrenreieh 
(Ref.  37),  and  other  publications  (Refs.  38  and  39). 

Under  U.S.  Air  Force  Contract  No.  AF33(657)-8935,  a systematic  cor- 
relation (Ref.  38)  was  developed  between  microwave  frequency  sound  atten- 
uation in  insulators  and  thermal  properties  of  these  cyrstals;  particularly 
thermal  conductivity  and  Debye  temperature.  This  work  made  it  possible  to 
rank  crystals  in  the  order  of  their  expected  high-frequency  acoustic  attenuation 
using  a parameter,  designated  R,  which  can  be  calculated  from  thermal  prop- 
erties such  as  thermal  conductivity  and  Debye  temperature.  This  ranking  not 
only  served  to  systematize  available  data  but  also  afforded  a means  for  pre- 
dicting the  acoustic  loss  of  crystals  whose  acoustic  properties  have  not  been 
measured.  Crystals  which  have  lew  sound  attenuation  at  microwave  frequencies 
and  room  temperature  were  shown  to  be  those  with  light  mass  atoms,  complex 
crystal  structures,  high  Debye  temperatures,  and  low  thermal  conductivities. 

A basis  was  provided  for; 

• Selecting  crystals  for  low  sound  attenuation  at  room  temperature 
from  thermal  measurements. 

• Selecting  host  crystals  in  which  reduced  acoustic  attenuation,  par- 
ticularly near  room  temperature,  might  be  achieved. 

t Selectively  doping  or  alloying  the  crystals  to  achieve  a reduction 
of  sound  attenuation  at  room  temperature. 

Further  support  through  U.S.  Air  Force  Contract  No.  AF33(615)-2590 
made  it  possible  to  refine  the  theory  of  microwave  acoustic  attenuation  in 
crystals,  to  study  the  reduction  of  room -temperature  sound  attenuation  by  the 
presence  of  paramagnetic  ions  in  spinel  containing  Fe8+,  to  rank  borides  from 
thermal  measurements,  and  to  grow  yttrium -aluminum  garnet  as  a host  crystal 
for  further  studies  to  reduce  acoustic  attenuation.  Sapphire,  yttrium -aluminum 
garnet,  and  magnesium  aluminate  spinel  were  identified  as  crystals  with  acoustic 
losses  near  the  lower  limit  for  oxides.  Icosahedral  boron  compounds  were  iden- 
tified as  a system  with  potentially  lower  losses  than  the  oxides,  and  YB^,  was 
selected  as  one  of  the  best  candidates  among  borides  because  of  its  unusual 
structure.  Under  U.S.  Air  Force  Contract  No.  F3361 5-67 -C -1399  (Ref.  40), 
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comparison  of  attenuation  was  made  in  MgAl804  spinel  (Ref.  41)  and  in  yttrium 
and  rare  earth  garnets,  and  preliminary  steps  were  taken  to  grow  YB^  single 
crystals.  Flux  grown  spinel,  MgAip4(Ref.  42)  was  shown  to  be  the  material 
with  lowest  acoustic  attenuation  at  room  temperature. 

Further  studies,  under  Contract  No.  F33815-69-C-1286  (Ref.  43)  provided 
the  first  intercomparison  of  sound  attenuation  in  very -low -loss  materials  grown 
by  several  different  techniques.  Comparisons  were  made  among  sapphire 
crystals  grown  by  Verneuil,  Czochralski,  and  chemical  vapor  deposition  tech- 
niques, and  MgAla04  spinel  grown  by  flux  evaporation  and  Czochralski  tech- 
niques. Spinel,  once  again,  was  shown  to  be  the  lowest-loss  material,  but 
the  Czochralski  grown  crystals  available  at  that  time  were  found  to  be  unsatis- 
factory because  of  the  high  density  of  precipitates  within  them.  The  available 
flux  grown  crystals  have  been  too  small  for  many  practical  delay  line  require- 
ments. Recent  developments  in  the  growth  of  spinel  as  a substrate  for  silicon 
epitaxy  have  made  it  worthwhile  to  reevaluate  Czochralski  spinel  as  an  acoustic 
delay  line  material. 


Under  Contract  No.  F33615-69-C-1286  (Ref.  43)  methods  were  developed 
for  synthesizing  and  growing  single  crystals  of  YB^  and  the  compound  was 
shown  to  melt  congruently.  Attenuation  measurements  showed  that  a low -tem- 
perature relaxation  process  dominates  the  sound  attenuation.  No  estimate  of 
the  thermal  phonon  related  losses  was  possible  because  of  the  relaxation  attenua- 
tion. The  crystal  stoichiometry,  uniformity  of  composition  in  individual  crystals, 
and  purity  were  not  under  sufficient  control  to  permit  identification  of  the  re- 
laxation mechanism.  It  was  necessary  to  identify  and  eliminate  the  cause  of 
the  relaxation  absorption  so  that  the  potential  of  YBg,  as  an  acoustic  material 
could  be  realized , 


During  the  present  contract  the  stoichiometry  and  purity  of  YB^  crystal 
has  been  examined.  The  compound  has  been  shown  to  melt  congruently  at  a 
composition  near  B/Y  = 81.7.  Small  crystals  of  a number  of  compositions 
have  been  grown  and  characterized,  including  the  congruently  melting  com- 
position. Precipitates  of  other  Y -B  phases  have  been  identified  in  crystals 
of  Borne  compositions,  and  growth  techniques  have  been  adapted  to  eliminate 
them.  The  effort  to  grow  larger  crystals  was  warranted,  and  equipment  modi- 
fications were  made  to  scale  up  the  length  of  the  crystals  grown.  These  results 
and  characterization  of  the  crystals  are  described  in  this  section  of  the  report. 


In  Section  2,  the  design  study  of  delay  lines  showed  that  MgAl804  spinel 
was  superior  to  other  oxides  as  a delay  line  medium.  To  realize  its  potential 
it  is  necessary  to  (1)  obtain  large,  high-quality  boules,  (2)  provide  shear 
mode  transducers  which  can  be  manufactured  reproducibly,  and  (3)  to  use 
electrodes  with  lower  acoustic  attenuation  than  the  metals  presently  used. 
Since  YB^  is  potentially  a lower -loss  material  than  MgAla04,  the  three  con- 
clusions above  apply  equally  to  YB^  once  the  relaxation  losses  are  eliminated. 
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B GROWTH  METHODS 


The  synthesis  and  growth  of  single -crystal  YB^  has  been  reported  pre- 
viously (Kefs.  4 3-45).  Since  the  techniques  used  are  not  widely  applied  they 
will  be  described  here  briefly  to  make  the  remainder  of  Section  4 more  in- 
telligible. Alterations  of  the  growth  procedures  are  described  which  will  per- 
mit growth  of  crystals  with  improved  uniformity  with  a composition  very  close 
to  the  congruent  melting  composition. 

Synthesis 

YB^  crystals  are  synthesized  directly  from  the  elements  yttrium  and 
boron.  The  reactivity  of  the  components  eliminates  the  use  of  heated  ceramic 
or  metal  crucibles.  Synthesis  is  complicated  by  the  high  heat  of  formation  of 
an  intermediate  compound  YB4  which  forms  with  explosive  violence  if  the  re- 
action rate  is  not  limited  by  the  synthesis  technique.  The  material  is  synthesized 
in  a water-cooled  copper  boat  (Kefs.  44  and  46),  Heating  is  accomplished  with 
RF  induction,  and  coupling  to  the  charge  is  initiated  by  preheating  the  charge  to 
about  1000°  C with  a removable  susceptor.  Chemical  interaction  with  the  crucible 
is  negligible  (Refs.  44  and  47). 

The  exothermic  formation  of  YB4  is  controlled  by  placing  pieces  of  yttrium 
metal  in  appropriate  quantity  on  blocks  of  boron,  rather  than  by  mixing  the  ele- 
ments intimately.  Upon  heating,  a molten  pool  of  eutectic  is  formed . As  the 
temperature  is  increased,  YB4  forms  at  the  interface  of  the  solid  boron  and  the 
eutectic  pool:  its  formation  is  readily  observed  from  indandescence  at  the  edge 
of  the  pool.  In  the  geometry  employed,  the  reaction  rate  is  limited  by  diffusion 
of  boron  and  yttrium  and  by  the  large  heat  flow  to  the  water-cooled  boat.  Manual 
power  adjustment  is  sufficient  to  control  the  reaction  under  these  conditions. 

Directionally  Solidified  Ingots 

The  reacted  YBff  from  the  cold  boat  is  polycrystalline,  and  it  fractures 
on  reheating.  A pulling  arrangement  was  added  to  the  cold  boat  so  that  seeded 
ingots  could  be  grown  by  the  Czochralski  method.  These  ingots  are  typically 
single -crystal  but  have  appreciable  mosaic  structure  due  to  the  asymmetric 
temperature  distribution  of  the  cold  metal  boat  system. 

Pedestal  Technique 


Crystals  for  study  are  prepared  by  using  the  ingots  pulled  from  the  cold 
boat  as  pedeRtals  (Kefs.  4 5,  48,  and  49).  The  ingots  are  supported  from  the 
lower  end  with  a pyrolytic  BN  support.  RF  heating  of  the  upper  end  of  the 
ingot  is  initiated  with  a removable  susceptor,  the  molten  pool  is  seeded  with  a 
piece  of  single  crystal,  and  the  crystal  is  grown  by  passing  the  liquid  region  through 
the  ingot.  This  procedure  is  equivalent  to  a single  zone  refining  operation  and 
has  been  advantageous  in  the  studies  of  stoichiometry  reported  below. 


The  basic  characteristics  of  the  pedestal  method  can  be  understood  by 
reference  to  Figure  32,  a simple  phase  diagram  for  a two-component  system 
which  forms  eutectics  at  compositions  I and  III  and  solidifies  congruently  -- 
liquid  and  solid  have  the  same  composition  at  composition  n.  The  congruent 
melting  composition  need  not  correspond  to  ideal  stoichiometry  of  the  solid, 
though  it  is  frequently  close  to  the  ideal  composition.  The  structure  of  YB^ 
is  sufficiently  complex  and  the  site  occupation  is  partially  statistical  so  that 
an  ideal  stoichiometry  has  never  been  previously  defined.  It  has  therefore 
been  desirable  to  establish  the  congruent  melting  composition,  to  prepare 
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Figure  32.  Binary  Phase  Diagram  with  Two  Eutectics  at  I and  III 
and  a Congruently  Melting  Compound  of  Composition 
AB„  at  II 

crystals  of  uniform  composition,  and  to  establish  the  relation  of  the  congruent 
melting  composition  to  ideal  stoichiometry  and  ultrasonic  relaxation  phenomena. 

If  one  starts  with  a pedestal  of  composition  (C)  as  shown  on  the  ordinate 
of  Figure  32,  melting  occurs  at  point  Lt.  However,  solidified  crystal  forms 
at  composition  Sj  as  the  molten  region  passes  through  the  rod.  Since  the  solid 
being  formed  is  richer  in  component  A than  is  the  pedestal  being  melted,  the 


liquid  is  depleted  of  component  A an  pedestal  growth  continues.  At  some  later 
time  the  liquid  will  be  at  point  La  and  the  solid  at  point  Sa,  the  same  composition 
as  the  pedestal.  This  condition  will  remain  stable  until  the  pedestal  is  exhausted 
and  there  is  no  more  pedesta’  to  melt.  The  very  last  material  to  freeze  will  be 
at  the  eutectic?  this  will  be  a mixture  of  the  compound  at  Sf  and  constituent  B 
at  Sf. 


A single  solidification  with  the  pedestal  technique  then  leads  to  nonuniform 
crystals,  unless  one  is  using  a pedestal  of  the  congruent  melting  composition. 
This  fact  is  extremely  useful  in  establishing  the  congruent  melting  composition. 
Comparison  of  the  first  material  to  freeze  with  later  freezing  material  gives  an 
extremely  sensitive  test  of  the  departure  of  the  pedestal  composition  from  the 
congruent  melting  composition.  Metallographic  examination  of  the  final  solid- 
ified portion  shows  unambiguously  which  eutectic  has  been  reached. 

Repetitive  zone  refining  (multiple  pedestal  growth)  produces  a crystal  of 
which  a part  can  be  very  uniform  and  very  nearly  the  congruent  melting  com- 
position. If  one  begins  with  a long  pedestal  of  composition  Z:  (Figure  32)  and 
performs  multiple  zone  passes,  the  liquid  composition  at  the  upper  end  is  of 
composition  Z1#  Zfl,  Z3  at  the  start  of  the  first,  second,  and  third  zone  passes 
and  the  first  material  to  freeze  is  of  compositions  ZSj,  ZSa,  and  ZS3.  After 
a large  number  of  zone  passes,  the  upper  end  of  the  crystal  will  be  near  the 
congruent  melting  composition  and  the  lower  end  will  depart  radically  from 
this  composition.  If  the  crystal  is  made  sufficiently  long  bo  that  the  upper  20  per- 
cent of  the  crystal  is  adequate  in  length  for  measurements,  then  that  sample 
portion  can  be  made  very  uniform  and  near  the  congruent  melting  composition 
by  multiple  pass  zone  refining.  The  pedestal  apparatus  in  use  is  suitable  for 
zone  refining;  one  need  only  have  a pedestal  of  the  order  of  four  to  six  inches 
in  length. 

Crystals  of  uniform  composition  other  than  the  congruent  melting  com- 
position can  be  obtained  by  Czochralski  growth  of  a small  crystal  from  a 
large  volume  of  melt,  or  by  zone  leveling  a pedestal  which,  initially,  was  nonuni- 
form in  composition. 

Other  Growth  Techniques 

Most  other  growth  techniques  either  are  not  applicable  to  YB^  growth, 
because  of  the  reactivity  of  the  melt,  or  are  not  advantageous.  One  additional 
method  is  useful:  Czochralski  growth  from  a Hukin  cold  crucible  (Ref.  50). 
This  crucible  will  levitate  a large  volume  of  melt,  and  crystal  growth  takes 
place  in  a symmetric  temperature  distribution.  Such  a system  should  provide 
long  pedestals  for  zone  refining  and  zone  leveling.  It  has  not  been  used  in  this 
work  because  it  requires  a large  volume  of  melt  to  function  properly,  and  such 
a quantity  of  melt  would  have  been  excessively  expensive.  The  Hukin  crucible 
has  not  been  used  previously  for  materials  which  must  be  preheated  to  initiate 
coupling,  and  this  function  would  have  to  be  provided. 
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Production  of  Crystals 


Crystal  growth  has  been  continuously  improved  by  the  use  of  charges  of 
the  congruently  melting  composition  and  yttrium  of  higher  purity  than  previously 
available.  The  yttrium  metal  was  distilled  for  use  on  this  project  in  an  ultra- 
high-vacuum  system,  by  G.  W.  Demarko,  General  Electric  Neutron  Devices 
Department,  through  arrangements  made  with  the  help  of  R.  L,  Hickmott, 
(Wright-Patterson  Air  Force  Base)  and  the  U.  S.  Atomic  Energy  Commission. 
An  analysis  of  the  starting  material  and  the  distilled  metal  is  given  in  Table  10. 
Vacuum  fusion  measurements  were  used  for  oxygen,  nitrogen,  and  hydrogen. 
Other  impurity  concentrations  were  determined  by  spark  source  mass  spectro- 
scopy. A list  of  crystal  growth  runs  is  given  in  Table  11. 

Table  10 

IMPURITIES  IN  YTTRIUM  METAL 


Impurity 

WelKht  Percent 
Starting  Material 
(ppm) 

Wrujlil  Percent  In 
IJlHtllle.l  Yttrium 
(ppm) 

c 

4 

4 

O 

532 

174 

N, 

4(1 

2 

H, 

30 

10 

Tn 

7H0 

47 

Cu 

550 

13 

Ni 

4 

1H 

l‘Y 

24 

12 

Cr 

!i 

2 

•ri 

180 

40 

Sc 

110 

«-7 

Ca 

7 

4 

Cl 

ni 

2 

SI 

27 

23 

A1 

5 

* :i 

v 1 

1 

X 

< I 

1 

. i 

Table  11 

CRYSTAL  GROWTH  RUNS 


Run 

Comporltlon  11/ Y 

Apparatus 

P„ 

01.7/1 

L'oUl  l>« < nt 

P,. 

• 

>% 

1’,. 

PcHrstol 

P-. 

P-, 

■ 

Cnltl  boat 

P.p 

■ 

72 


The  stoichiometric  crystals  grown  from  high  purity  yttrium  were  used  in 
sound  attenuation  measurements  to  establish  the  diffraction  losses  and  intrinsic 
losses  in  YB4#(see  Section  6). 

D.  CRYSTAL  CWWSIW  AND  WTERSTTT1AL  METAL  ATOMS 

There  have  been  several  papers  in  the  literature  on  YB^  and  related  com- 
pounds. The  data  of  Spear  and  Solovyev  (Ref.  51)  and  of  Schwetz  et  al.  (Ref.  52) 
show  that  all  of  the  rare  earths  from  neodymium  to  lutetium  (except  europium) 
make  YBM  phases.  The  larger  rare  earths  (lathanlum.  cerium,  and  praseo- 
dymium) do  not.  The  work  of  Carlsson  and  Lundstrom  (Ref.  53)  shows  that 
scandium,  manganese,  iron,  and  other  elements  do  not.  These  latter  elements 
can  be  incorporated  into  the  /3-boron  structure  and  cause  a lattice  expansion, 
but  do  not  change  the  crystal  structure. 

It  appears  to  take  r,  certain  minimum  size  Interstitial  atom  to  trigger  the 
crystal  structure  change  from  the  rhombohedral  structure  of  0-boron  to  the 
cubic  structure  of  YBM.  This  minimum  size  is  somewhere  between  that  of 
scandium  and  lutetium,  or  between  1.63  and  1.72  angstroms.  If  1.68  angstroms 
is  taken  as  the  minimum  radius,  then  from  Figure  3 of  Reference  52  the  mini- 
mum lattice  constant  of  this  phase  is  23.370  angstroms.  The  maximum  radius 
appears  to  be  thorium  at  1.81  angstrom,  and  the  maximum  lattice  constant 
(Ref.  54)  is  23.53  angstroms. 

Assuming  that  all  of  these  cubic  crystals  have  the  same  number  of  boron 
atoms  per  unit  cell  of  N,  * 1632,  the  average  volume  per  boron  atom  ranges 
from  7.821  (X)3  to  7.983  (l)3.  In  0-boron  with  105  boron  atoms  per  unit  cell, 
the  volumes  range  (Ref.  53)  from  7.814  (X)3  at  pure  0-boron  to  7.969(X)3  at 
manganese  saturated  0-boron.  Thus  the  density  of  packing  of  boron  atoms  in 
the  two  structures  is  almost  the  same,  and  the  ranges  of  existence  are  com- 
parable. For  comparison,  the  very  dense  o-boron  phase  (Ref.  55)  has  a 
volume  per  boron  atom  of  7.267  (J)3,  which  is  noticeably  less  than  for  either 
0-boron  or  the  YBae-type  crystals. 

Whether  a particular  metal-boron  composition  will  occur  as  the  0-boron 
or  the  cubic  YBM  structures  is  determined  by  the  size  of  the  interstitial  metal 
atom.  The  0-boron  range  (Ref.  54)  is  from  1.  63  angstrom  radius  at  scandium 
down  to  or  below  1.  24  angstroms  at  nickel.  To  obtain  the  cubic  crystal  struc- 
ture for  ultrasonic  delay  line  applications,  it  is  necessary  to  put  in  the  larger 
metal  atoms  --  yttrium,  thorium,  plutonium,  or  the  rare  earths. 

E.  CRYSTAL  STQM^mX 

Uniformity  of  composition  and  crystal  stoichiometry  are  important  in  the 
acoustic  applications  of  the  material.  The  composition  of  the  compound  des- 
ignated YBgs  has  been  reported  as  YB100  (Ref.  56),  YBW  (Ref.  57),  and  YB^ 

(Ref.  58).  The  latter  results  are  the  most  accurate  and  are  the  result  of  an 
x-ray  structure  determination.  The  complexity  of  the  structure  and  the  statis- 
tical occupation  of  a number  of  the  sites  precluded  an  ideal  stoichiometry  from 
the  structure  data,  and  it  was  considered  to  be  uncertain  to  within  about  66  t 2. 
(Ref.  59). 
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The  specimen  studied  in  the  x-ray  structure  determination  was  grown  by 
Seybolt  (Ref.  56)  and  formed  in  a (3-boron- YB,a  eutectic.  It  must  be  considered 
to  be  rich  in  boron.  The  structure  suggests  that  there  may  be  an  appreciable 
solubility  of  yttrium  in  the  crystal.  X-ruy  lattice  parameter  measurements 
showed  a range  of  values  from  23.  446  * 0,  005  angstroms  to  23.  482  f 0.  003 
angstroms,  and  the  x-ray  density  (Hof.  58)  2.482  g/cm3  differed  appreciably 
from  a pycnometric  density  of  2.  568  t 0.  005  (Refs.  43  and  60). 

The  congruent  melting  point  has  been  located  by  pedestal  crystal  growth, 
precise  density  comparison,  and  metallographic  examination.  (The  utility  of 
pedestal  growth  and  metallography  is  explained  above.)  Crystals  have  also 
been  examined  by  making  precise  x-ray  lattice  parameter  determinations  and 
analytical  wet  chemical  analysis  of  composition.  Only  the  density  measurements 
had  the  requisite  precision  to  lead  to  the  congruent  melting  composition.  Sim- 
ilar experience  has  been  reported  by  other  authors,  (Refs.  61-G3),  who  have 
had  to  rely  upon  parameters  such  as  ferroelectric  transition  temperature, 

NMR  line  width,  and  density,  which  are  very  sensitive  to  small  changes  in 
stoichiometry. 

A simple  density  apparatus,  shown  in  Figure  33,  was  constructed  to  measure 
the  relative  density  of  samples  sectioned  from  a crystal  boule.  A constant-tem- 
perature water  bath  was  made  by  fitting  a glass  tank  with  an  inner  glass  cylinder 
such  that  the  areas  inside  and  outside  the  inner  cylinder  were  equal.  With  this 
arrangement  a simple  magnetic  stirring  motor  and  impeller  provided  rapid 
water  circulation.  Temperature  control  was  provided  by  inserting  a thermistor 
with  millidegree  sensitivity  through  the  plastic  foam  top.  The  bath  temperature 
was  maintained  constant  by  a proportional  controller  which  drove  a heating  ele- 
ment made  from  a Ni-Cr  heater  (Amperex  Electronic  Corporation,  Itauppnuge, 
New  York)  with  an  Inconel  jacket  0.058  inch  outside  diameter.  Temperature 
was  monitored  with  a mercury-in-glass  thermometer  calibrated  in  hundredths 
of  a degree  centigrade. 

The  sample  tube  was  a test  tube  with  a ground,  tapered  upper  rnd.  A 
tapered  glass  stopper  was  fitted  to  a capillary  tube  and  inserted  in  the  test  tube. 
The  test  tube  was  loaded  with  the  samples  to  be  measured,  filled  with  a thallium 
formate-water  solution,  and  closed  with  the  glass  stopper.  Density  was  mea- 
sured by  adjusting  the  temperature  until  samples  floated  at  a specific  height  in 
the  tube.  The  sample  height,  temperature,  and  height  of  the  solution  in  the 
capillary  were  monitored  with  a cathetomcter  to  eliminate  parallax.  I'he  volume 
of  the  tube  and  capillary  were  calibrated  by  careful  weighing  and  tpmpernturc 
scanning  with  distilled  water. 


Absolute  densities  were  determined  from  the  weight  of  solution  added  to  the 
tube  and  capillary  and  from  tie  observed  volume.  These  values  were  consistent 
with  those  found  for  the  solution  density  using  a pycnometer  bottle.  Greater  ac- 
curacy was  obtained  in  the  relative  density  of  samples  observed  during  a single 
run  by  calculating  the  density  change  of  the  fluid  from  the  column  height  in  the 
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capillary  tube.  It  was  necessary  to  carefully  wet  the  samples  with  fluid  and  to 
expose  them  to  ultrasonic  agitation  to  ensure  that  no  gas  bubbles  adhered  to  them 
and  caused  erroneous  results.  The  absolute  density  values  are  accurate  to  ap- 
proximately 0.2  percent  and  the  relative  value  to  0.01  percent. 

During  the  course  of  this  study  approximately  twelve  synthesis  runs,  29 
crystal  growth  runs,  and  128  density  determinations  have  been  carried  cut. 
Selected  data  are  summarized  in  Table  12  with  x-ray  lattice  parameters  and 
wet  chemical  analysis. 

Inspection  of  Table  12  shows  that  the  analytical  measurements  of  com- 
position are  too  inaccurate  to  assist  in  locating  the  congruent  melting  com- 
position. There  is  little  correspondence  between  the  chemical  formula  cal- 
culated from  these  results  and  the  starting  composition  of  the  melt  used  in  the 
initial  synthesis  of  the  compound  and  listed  in  the  first  column,  "Composition. " 

Table  12 

DENSITY,  X-RAY,  CHEMICAL  ANALYSIS  RESULTS 


Composition 

Density 

*0 

Analytical  Composition 
Weight  Percent  V and  Equivalent  CompoHitiou 

Yh*>" 

• - 

23.  482  t 0.  003 

• • 

vn40 

-- 

23.  471  + 0.  004 

-- 

YfW 

2.  004  1 

23.455  i 0.005 

2.  5782 

-- 

YtW 

2.  57  f>  1 

23.443  t 0.004 

YB.* 

2.  577  0 

-- 

YB*„ 

2.  5730 

23.451  l 0.  005 

y«*t 

2.  5706 

-- 

Y«W 

2.  5737 

-- 

Y'W 

2.  5747 

-- 

Yll„, 

2 5712 

-- 

-- 

' B«.» 

2.  5700 

-- 

-- 

2.  5685 

23.  445  t 0.  005 

-- 

YB* 

2.  5064 

-- 

-- 

YB* 

2.  5700 

23.  442  i 0.  005 

-- 

YV 

2.  5622 

-- 

YB*f 

2.  5609 

23.448  1 0.005 

-- 

^ H/vtfti) 

2 5676 

-- 

-- 

YIV. 

2.  5576 

23.  443  * 0.  005 

•- 

YIW 

2.  5686 

23.  443  1 0.  004 

-- 

' ! V' 

- 

23.440  ! 0 006 

v«V.* 

2.  5749 

-- 

12.  24  * 0 t V n„-(. 

YB*,4 

2 5782 

12  13  • 0.  1 YM„  . 

\ H .r' 

2.  5736 

12  18  • 0.  1 Ylt.#i, 

Sample  contained  srcimri  phase 
Sample  contained  -econd  phase  (1-lioron. 

♦Two  200-mg  samples  sectioned  along  most  of  the  length  of  the 
crystal  am!  matched  In  density. 
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The  densities  of  these  samples  correspond  satisfactorily  with  the  melt  com- 
positions. The  x-ray  lattice  parameters  show  a definite  increase  for  the 
samples  richest  in  yttrium.  However,  the  change  in  lattice  parameter  from 
the  most  boron-rich  sample  through  the  congruent  melting  composition  is  of 
the  order  of  or  smaller  than  the  standard  deviation  of  the  me  astir  ements.  A 
factor  of  six  improvement  in  accuracy  of  the  lattice  parameter  is  required  to 
make  it  comparable  to  the  density  data.  Such  an  improvement  would  be  ad- 
vantageous in  allowing  one  to  assess  the  relative  roles  of  lattice  parameter 
change  and  composition  change  in  affecting  density. 

The  density  data  may  be  understood  more  completely  by  recalling  the  de- 
scription of  the  pedestal  growth  method  (page  69  ).  The  first  material  to  freeze 
is  closer  to  the  congruent  melting  composition  than  the  parent  liquid;  the  com- 
position of  the  melt  during  the  synthesis  process  is  listed  in  the  first  column 
of  Table  12.  The  pedestal  grown  crystals  were  grown  from  melted  pedestal, 
and  that  liquid  is  expected  to  be  somewhat  closer  to  the  congruent  melting  com- 
position than  the  synthesis  melt  was.  because  of  segregation. 

All  of  the  samples  labeled  T (Table  12)  were  sectioned  from  the  very  top 
of  the  crystals.  They  are  expected  to  show  only  a gradual  increase  in  density 
as  the  quantity  of  yttrium  in  the  melt  increases  and,  correspondingly,  the 
amount  of  yttrium  in  the  crystal  is  increased.  This  dependence  is  shown  by 
the  solid  circles  in  Figure  34.  The  samples  labeled  B (Table  12)  were  sec- 
tioned from  the  bottoms  of  the  crystals.  The  pedestal  growth  process  causes 
the  very  bottom  of  the  crystal  to  reach  the  eutectic  point  on  the  same  side  of 
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Figure  34.  Sample  Density  Versus  Composition  of  Melt  During  Synthesis 
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the  congruent  melting  composition  as  the  pedestal  composition.  If  the  pedestal 
was  boron-rich  the  bottom  of  the  crystal  should  reach  the  boron-rich  eutectic 
and  have  a low  density.  If  the  pedestal  was  yttrium-rich  it  should  reach  the 
YBjg-YBjg  eutectic  and  have  a high  density.  Exactly  this  behavior  is  shown 
by  the  open  squares  in  Figure  34.  Bottom  samples  from  melts  very  rich  in 
yttrium  showed  a second  phase  of  YB19  and  bottom  samples  from  melts  very 
rich  in  boron  showed  0-boron  as  a second  phase.  Clearly  the  congruent  melting 
composition  is  near  but  slightly  more  boron-rich  than 

The  congruent  melting  composition  can  be  estimated  more  accurately  by 
plotting  symmetrically  the  bath  temperature  difference  for  floating  top  and 
bottom  samples  versus  the  composition,  as  is  done  in  Figure  35.  The  con- 
gruent melting  composition  is  estimated  to  be  YB81t7Sdfco#1.  The  remaining  un- 
certainty in  composition  will  strain  the  density  method  and  may  be  influenced 
by  the  segregation  of  impurities . The  most  reliable  method  of  producing  a 
uniform  crystal  of  CM  composition  is  to  apply  multiple  zone  refining  to  a ped- 
estal grown  from  a melt  of  composition  61.75.  A crystal  of  composition  61.5 
was  grown  with  the  intent  of  fabricating  an  acoustic  sample  for  attenuation 
measurement.  Subsequent  examination  showed  appreciable  mosaic  structure 
to  be  present,  so  that  the  sample  has  been  used  instead  for  electrical  resistivity 
measurements. 


Figure  35.  Temperature  Difference  for  Floating  Samples  from  the  Top  and 
Bottom  of  a Crystal  (Plotted  Symmetrically)  Versus  the  Com- 
position of  the  Melt  Duri  .•  Synthesis 
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Metallographic  examination  of  the  yttrium- rich  samples  initially  produced 
a puzzle.  A second  phase  was  observed  as  expected,  but  its  color  (Ref.  45) 
suggested  that  it  was  YEa  rather  than  YBi8  as  expected.  Because  the  second 
phase  grains  were  of  micron  and  submicron  size  an  electron  beam  microprobe 
was  applied,  with  the  results  shown  in  Table  13.  The  precipitate  is  clearly 
YBia  in  agreement  with  the  phase  diagram,  and  the  subjective  color  identifica- 
tion was  erroneous. 


Table  13 

ELECTRON  BEAM  MICROPROBE  EXAMINATION 
OF  YTTRIUM -RICH  SECOND  PHASE 


Katio  V C ounts/  11  Counts 
standard  Sample 

Hatio  Y Counts/ B Counts 
Mptallofiruphlr  Sample 

YBn  0.373  i 0.0045 

Matrix  0.320  t 0.0030 

YU„  3.012  t 0.0205 

Precipitate  2.344  t 0.024H 

Y»r  4.274  t 0.0402 

Yttrium  ratio 

CountH  rnt»c» 

60/ 1 2 - 5.5 

7.4  metallographic  sample 

5.4  standard 

66/ 6 = 11 

11.4  Mtnndard 

F.  PENSfTY.  ANP  ..CPMPQSrnpN  - A MQPEl 


The  progress  in  the  growth  of  single  crystals  of  YB^,  enables  one  to  draw 
some  conclusions  about  the  crystal  structure.  The  congruent  melting  composi- 
tion was  found  to  be  YBr,  with  n = 61.75.  Samples  had  a density  of  2.5687 
±0.005  g/cm3  and  an  x-ray  lattice  parameter  of  a0  = 23.445  ± 0.005  ang- 
stroms. The  density  can  be  calculated  from 


0s 


m,  • NB  + mv  • Nv 
” (ao“>4  Na 


(55) 


where 

mB  = atomic  mass  of  boron 
m¥  = atomic  mass  of  yttrium 
Ne  = number  of  boron  atoms  in  a unit  cell 
Ny  = number  of  yttrium  atoms  in  a unit  cell 
N*  = Avogadros’  number 
This  can  be  rewritten  as 


0 = 


™aN8 

(a0)3N, 


J + niyNy 

mBNB 


(56) 
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Using  Nb/Ny  = 61.75  and  the  measured  values  of  p and  Sq,  the  calculated 
value  of  Nb  is 

Nb  = 1628  ± 4 (57) 

This  value  is  higher  than  that  given  by  Richards  and  Kasper  (Refs.  58  and  59), 
who  postulated 


Na  = 1584  ± 48 

The  difference  can  be  explained  if  it  is  assumed  for  the  present  crystals,  which 
are  grown  slowly  and  close  to  equilibrium  conditions,  that  in  a unit  cell  of  vol- 
ume (ao)®  there  are  8 units  of  (Bia)u  icosahedra  plus  8 units  of  configuration  II 
in  the  channels  (Ref,  58).  These  units  have  48  boron  atoms  each.  Thus  NB  * 

8 (156+48)  * 1632.  This  model  avoids  the  necessity  for  putting  units  of  both 
configuration  I (with  36  atoms)  and  configuration  II  in  the  channels  statistically. 

Now  the  calculated  density,  with  N,  = 1632  and  a0  = 23.445  angstroms, 
comes  out  to  be  (x-ray)  = 2.575  g/cm®,  at  n = 61.75,  which  is  in  good  agree- 
ment with  the  measured  value  of  2.569  ± 0.005  g/mm®.  The  next  step  is  to 
compute  NYi 

N>  = 6Trr5  = 26-4  (58) 

Since  there  are  48  possible  yttrium  sites  per  unit  cell,  then  only  55.0  percent 
of  them  are  occupied  at  n = 61.75.  If  exactly  half  of  them  were  occupied,  the 
chemical  formula  would  be  YBgg.  Therefore  YBe8  probably  should  represent 
the  stoichiometric  compound.  The  congruent  melting  point  at  61.75  is  slightly 
shifted  away  from  this  value  of  n. 

If  Nb  = 1632  and  the  possible  range  of  Nv  is  0 < Ny  <48,  and  if  the  boron 
lattice  parameter  is  independent  of  the  yttrium  concentration  in  this  range, 
then  the  range  of  densities  versus  Ny  would  be 

p=  2.2724  [l+(Ny/201. 7)]  g/cm*  (59) 

The  data  on  a0  versus  density  (see  Table  12)  show  that  a0  is  essentially  in- 
dependent of  yttrium  concentration  from  densities  of  2.558  to  2.575  g/cm®. 

This  density  range  convert::  to  a range  of25.3<NY<26.9  atoms  in  a unit  cell. 

The  model  of  Richards  and  Kasper  states  that  the  maximum  value  of  Ny  is 
24;  that  is,  only  50  percent  of  the  sites  are  filled.  For  NY>24  the  yttrium 
atoms  in  adjacent  octants  may  begin  to  interfere  with  one  another.  The  inter- 
ference for  25  <Ny  <27  is  assumed  to  have  only  a small  effect  on  ao  for  the 
present  purposes:  this  is  justified  in  a later  section. 

The  present  model  of  YBn  shows  that  there  is  very  little  disorder  in  the 
boron  lattice  of  Np  = 1632  atoms  per  cell.  All  of  the  variation  seems  to  lie  in  Ny. 
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The  model  of  Richards  and  Kasper  shows  that  for  each  yttrium  atom  present 
in  a lattice  site  there  is  another  equally  possible  vacant  site  located  2. 55  ang- 
stroms away  along  a [100]  direction.  The  Y radius  in  yttrium  metal  is  1.78 
angstroms,  comparable  to  the  value  found  for  Y in  YBm.  Thus  it  is  not  pos- 
sible for  yttrium  atoms  to  occupy  both  adjacent  sites  without  some  considerable 
local  distortion  and  expansion  of  the  boron  lattice.  This  may  be  the  reason  for 
the  increase  in  ao  seen  for  Ny  approaching  48  in  Table  12. 

The  two  adjacent  yttrium  sites  are  separated  by  an  octagonal  ring,  or  win- 
dow, of  eight  boron  atoms.  The  radius  of  this  window  is  1.45  angstroms, 
slightly  less  than  the  radius  of  an  yttrium  atom  (see  Figure  36).  Under  stress, 

BORON 


Figure  36.  Cavity  Occupied  by  an  Yttrium  Atom  in  YB^.  The  Yttrium 
is  Shown  as  a Sphere,  the  Boron  as  a Solid  Mass;  the  Ap- 
propriate Neutral  Atomic  Radii  were  Used  in  this  Construction 

however,  an  yttrium  atom  might  well  move  through  this  window  in  the  [100] 
plane  into  its  mirror-image  vacant  site.  Such  stress-induced  motion  may  be 
responsible  for  some  of  the  acoustic  losses  reported  below. 

If,  as  suggested,  the  motion  of  yttrium  atoms  within  the  dumbbell-shaped 
holes  in  the  lattice  is  responsible  for  the  high  acoustic  losses,  then  it  might 
be  possible  to  fill  both  sides  of  all  of  the  dumbbells  at  NY  =48.  The  maximum 
a0  of  the  cubic  YBgg  phase  occurs  (Ref.  54)  for  thorium  where  a0  = 23.53  ang- 
stroms. The  structural  study  (Ref.  54)  of  this  crystal  showed  that  the  largest 
metal  atom  that  could  fit,  two  at  a time  into  the  dumbbell  cavity  would  have  a 
radius  of  1,67  angstroms.  Thus  it  might  be  possible  to  make  LuaBee  = I,uB34. 
If  one  could  choose  two  atoms  where  the  average  radius  was  1.67  angstroms, 
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then  a possibility  of  scandium  plus  lutetium  would  do.  One  might  try  to  grow 
ScLuB6#.  In  this  manner  it  might  be  possible  to  grow  a crystal  with  low  acoustic 
losses.  The  losses  in  YB^  are  too  high  to  be  useful. 

G.  SEGREGATION  COEFFICIENT 

The  authors  have  shown  that  YBn  for  n * 61.7  is  a congruently  melting 
compound  with  a melting  point  of  about  2100°C.  From  Figure  32  it  can  be  seen 
that,  for  a melt  composition  where  B/Y >61.7,  B will  be  rejected  into  the  melt 
from  the  growing  solid  crystal.  If  B/Y  <61.7,  then  the  solid  will  be  richer  in 
boron  than  the  melt.  If  one  defines  a segregation  coefficient  as  the  ratio 

v * deviation  of  boron  concentration  in  solid  from  n « 61.7 
118  deviation  of  boron  concentration  in  liquid  from  n = 61.7 

then  kg  is  approximately  0.3,  according  to  the  data  of  Figures  34  and  37.  This 
value  is  obtained  from  the  following  analysis. 


50  55  SO  65  7 0 

B/Y  RATIO  IN  MELT 

Figure  37.  Comparison  of  the  B/Y  H ratio  in  Crystal  Houles  Crown  from 
the  Melt  to  the  Starting  Composition  of  the  Ingot.  Samples 
Have  Been  Taken  frqm  Both  Top  rand  Bottom  of  the  Boulefi; 
Congruent  Melting  i-'oint  is  at  B/\  = 61,75 
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Assume  that  the  lattice  parameter  of  YBn  is  independent  of  n near  n - 61.7, 
This  assumption  is  borne  out  by  the  data,  as  discussed  previously.  Assume 
that  the  crystal  density  variation  is  only  caused  by  the  variation  in  Ny  < see 
equation  59).  Then  the  density  gives  a unique  value  for  the  B/Y  ratio.  Fig- 
ure 37,  which  is  a replot  of  Figure  34  using  equation  59  to  convert  density  to 
B/Y  ratio,  shows  the  results  for  crystals  that  have  been  regrown  twice  from 
melts  of  the  stated  composition.  The  top  part  of  the  boule  grew  first.  After 
two  growths  the  top  part  of  the  boule  shows  a deviation  from  the  congruent 
melting  composition  of  about  0.08  of  that  of  the  melt.  For  example,  a melt 
of  B/Y  * 50  produces  a crystal  with  B/Y  = 60.7,  or  a deviation  ofl.OinB/Y. 
Thus  ka  e 1.0/ (61.7-50.0)  * 0.085.  The  square  on  k comes  from  two  successive 
growth  processes  from  the  starting  melt  composition.  Thus,  k~  0.3. 

With  k this  far  from  unity  it  was  rather  difficult  to  grow  large,  homogeneous 
crystals  of  YB,  at  any  composition  other  than  n “ 61.7.  For  example,  if  one  de- 
sires to  grow  in  a single  growth  step  by  a seeded  pulling  technique  (as  was  used 
here)  a crystal  of  YB#ff(with  NY  ■ 24),  the  starting  melt  composition  would  have 
to  have  been  about  B/Y  - 83. 

The  maximum  yttrium  concentration  in  any  YB^,  crystal  grown  from  the 
melt  would  correspond  (Ref.  43,  p.  113)  to  a YB^-YBia  eutectic  melt  at 
B/\  « 30.  With  k = 0.3  the  resultant  crystal  in  one  growth-step  would  have 
B/Y  *>  52  or  Ny  » 31.  Here  only  65  percent  of  the  yttrium  sites  would  be  filled. 
Testing  to  see  if  this  predicted  value  of  NY  « 31  agrees  with  other  experiments 
one  sees  from  Table  12  that  a melt  of  composition  B/Y  * 30  would  have  given 
a crystal  with  a0  * 23.476  angstroms.  If  it  is  assumed  that  a YfB^,  crystal 
would  have  ao  ■ 23.55  angstroms  (slightly  larger  than  ThBrt),  that  YBM  has 
®o  * 23.44  5 angstroms,  and  that  a0  increases  linearly  with  Y content,  then  from 
a0  “ 23.476  angstroms  one  calculates  * 31.  The  linear  increase  corre- 
sponds to  assuming  that  for  YB„  with  34  s ns 68  one  has 

a0  ° (23.445  + 4.4  x 10"® NY ) angstroms  (61) 

The  close  agreement  at  NY  e 31  is  fortuitous,  but  the  model  makes  some 
sense.  The  total  length  of  a dumbbell  of  two  yttrium  atoms  is  7.12  angstroms. 
The  length  of  the  cavity  in  ThBffi  is  6,68  angstroms.  Thus  the  a0  of  Y^ 
should  be  just  slightly  greater  than  that  of  ThB«,  at  23.53  angstroms. 

H-  ELASTIC  CONSTANTS 

The  cubic  crystal  YB^g  has  three  independent  elastic  constants.  These 
have  been  measured  by  the  pulse-echo  technique.  The  sound  velocity  of  long- 
itudinal waves  in  the  [100]  direction  gives  cu  = p V?  and  the  [100  j shear  wave 
velocity  gives  cw  = pVJ,  where  p is  the  crystal  density  of  2.568  g/emB.  The 
values  of  Vu  and  V,  (Ref.  43,  p,  36)  are  Vt  » (1.21  -t  0,05)  x 10p  rm/spc  and 
v.  = <0.70  t 0.02)  x 10*  cm/sec.  The  longitudinal  velocity  in  a [111]  axis  rod 
labeled  P-32P,  which  was  0.27  cm  in  diameter  and  1.27  cm  long,  was  men- 
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sured  at  60  MHz  at  300° K.  The  velocity  was  1.18  x 10*  cm/ sec.  Using  the 
equation  3 p v8  = c„  + 2cl0  + 4c*.  one  calculates  that  >0.13  scia  s + 0.85  in 
units  of  1019  dynes/cm8. 

Thus,  for  YBgg! 

cn  = 3.  80  ± 0.3  x 101#  dyne/ cm* 

c4*  “ 1.80  ± 0.08  x 10ia  dyne/cm8 

cw  * 0.4  ±0.5  xlO1*  dyne/cm 8 (62) 


These  values  give  for  the  elastic  anisotropy  constant,  A a 2c44/(cn  -cia),  a 
range  of  A 3 0.96  t 0.14.  Thus  YBgg  is  nearly  isotropic,  since  A is  close 
to  unity.  The  acoustic  diffraction  constant,  b,  (see  equation  8),  is  related  to 
A by 


b 


i ( 1-a\  „ 

3 V A JC«  cn(cn.C) 


1.38  [A-1-!] 


(63) 


Thus  the  value  of  b for  YB,^  is 


b » +0.05  0.19 


(64) 


This  means  that  the  diffraction  loss  from  the  acoustic  beam  in  YB^  may  be 
small  (see  Figure  4).  A better  value  of  cl8  is  needed  in  order  to  calculate  b 
with  certainty. 

i.  vzmjmsmite 

From  the  known  elastic  constants  at  0°K  one  can  calculate  the  Debye  tem- 
perature at  absolute  zero.  Since  the  Debye  temperature  is  in  the  range  of 
1000° K,  the  present  approximation  uses  the  room -temperature  elastic  constants. 
ThiB  approximation  is  supported  by  the  observation  that  the  shear  wave  veloc- 
ity, which  determines  c^,  changes  very  little  (Ref.  43,  p.  36)  between  300"  K 
and  4°K.  The  calculated  Debye  temperature  ie  thus 

fl0  a 1300  t 50" K (65) 

Some  preliminary  specific  heat  capacity  measurements  on  YBM  mnde  at  Stan- 
ford University  (Ref.  64)  are  in  general  agreement  with  this  value.  Samples 
of  YBW  and  YB6l<,  both  gave  0O  * 1200°  K. 

J.  ACOUSTIC  ATTENUATION 

The  only  previous  data  (Ref.  43.  p.  37)  on  the  acoustic  attenuation  In  YB*, 
were  taken  at  1.15  GHz  for  shear  waves  in  a [1001  direction  between  2 ' K and 
10°K.  The  attenuation  was  found  to  be  quite  large;  and  it  was  suggested  (Ref.  4 3) 
that  measurements  be  made  at  low  frequencies.  This  has  now  been  done,  and 
the  data  are  shown  in  Figures  38  through  40. 
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TEMPERATURE.  *K 

Figure  38.  Acoustic  Attenuation  Versus  Temperature  at  0.060  GHz  for  a 
Melt-Grown  Sample  of  YB,,;  Direction  of  Propagation  is  [100] 


TEMPERATURE,  ‘K 

Figure  39.  Acoustic  Attenuation  of  Longitudinal  Phonons  Versus  Tempera- 
ture at  0.065  GHz  and  at  its  Third  Harmonic  for  a Melt -Grown 
Sample  of  YB„it7;  Direction  of  Propagation  is  [100], 
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Figure  40.  Acoustic  Attenuation  of  Transverse  Phonon  Versus  Temperature 
at  0.06  GHz  in  a Melt-Grown  Sample  of  YBili7  (solid  curve)  and 
at  1.15  GHz  in  a Sample  of  YB*  (dashed  curve);  Direction  of 
Propagation  is  [100] 


Sound  attenuation  measurements  using  the  pulse -echo  technique  and  the 
exponential  comparator  have  been  described  in  a previous  report  (Bef.  65). 

In  order  to  make  sound  attenuation  measurements  at  lower  frequencies  and 
with  narrow  pulse  widths  the  instrumentation  was  modified  as  shown  in  Fig- 
ure 41.  Narrow  pulses  were  needed  because  the  specimens  were  short  and 
the  acoustic  velocities  were  high.  Zero  bias  into  the  cascaded  ring  modulator 
cuts  off  the  CW  signal  by  approximately  80  decibels.  Pulse  widths  as  short 
as  0.1  .us  can  be  generated  by  this  technique  and  amplified  up  to  the  10-watt 
level  by  the  wideband  power  amplifier.  The  coaxial  sample  holder  and  the  re- 
ceiver part  of  the  instrumentation  are  the  same  as  were  used  previously. 

The  transducers  were  quartz  with  overtone -polished,  gold  plated  faces, 
and  their  thickness  was  one-half  acoustic  wavelength  at  65  MHz.  Transducers  of 
of  X-cut  quartz  were  used  for  longitudinal  waves;  AC -cut  quartz  was  used  fcr 
shear  wave  generation.  Considerable  difficulty  was  experienced  in  making  re- 
liable bonds  between  the  quartz  transducer  and  the  boride  sample.  The  bonds 
needed  to  operate  from  room  temperature  down  to  2°K.  This  bonding  problem 
is  probably  caused  by  the  large  difference  in  thermal  expansion  between  quartz 
and  the  boride  samples.  Only  two  bonding  materials  were*  found,  out  of  the 
many  materials  tried,  which,  with  care,  made  reliable  bonds.  These  were 
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Figure  41.  Block  Diagram  of  Instrumentation  Used  for  Low-Frequency 
Narrow -Pulse -Width  Acoustic  Attenuation  Measurements 


NONAQ  Stopcock  Grease  (Fisher  Scientific  Company,  Fair  Lawn,  New  Jersey) 
and  200  silicon  fluid  (Dow  Corning  Corporation,  Midland,  Michigan),  with  a 
viscosity  of  20,000  centistokes. 

In  order  to  operate  at  the  third  harmonic  of  the  quartz  transducer  it  was 
necessary  to  change  the  resonant  circuit  and  the  coupling  of  the  transducer- 
coaxial  line  system.  The  two  methods  used  are  shown  in  Figure  42.  In  Fig- 
ure 42A  the  length  of  the  stainless  steel  coaxial  line  is  greater  than  l/4  wave- 
length at  65  MHz,  so  that  the  entire  transmission  line  can  be  series -resonated 
with  variable  capacitor  Cj  and  the  coupling  controlled  by  capacitor  C„.  This 
method  worked  well  when  the  Arenberg  pulse  generator  was  used  as  a high 
power  source  at  65  MHz.  The  circuit  also  resonates  at  the  3rd,  5th,  7th, 
etc.  harmonics  without  retuning,  but  the  losses  increase  with  frequency.  In 
addition,  the  3rd  harmonic  source,  having  less  output  power,  produced  only 
marginally  detected  echoes. 
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TO  HYBRID 


Figure  4 2.  Cross -sectional  View  of  the  Two  Methods  Used  to  Resonate 
and  Couple  the  Transducer-Coaxial  Line  System 

Resonating  the  transducer  with  Lt  and  C*  and  coupling  with  Ca  as  Bhown 
in  Figure  42B  improved  the  conversion  efficiency  so  that  echoes  were  de- 
tected at  IBS  MHz,  The  tuning  method  shown  in  Figure  42B  worked  equally 
well  at  65  or  195  MHzj  however,  it  was  necessary  to  change  Lx  and  Ct  when 
the  frequency  was  changed.  Neither  of  these  tuning  methods  had  the  conversion 
efficiency  of  the  reentrant  cavity  and  ZnO  transducer  combination  previously 
used  at  1 GHz. 

Results 


The  acoustic  attenuation  results  are  shown  in  Figures  38  through  40.  The 
interpretation  of  these  results  is,  at  best,  tentative.  There  are  two  features, 
however,  which  are  present  for  all  samples.  There  Is  an  attenuation  peak  at 
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20  to  30°  K that  covers  a large  temperature  range  and  rises  in  magnitude  aB 
the  frequency  increases;  above  200° K the  attenuation  rises  rapidly  with  increas- 
ing temperature,  with  no  apparent  plateau.  It  is  conjectured  that  the  25°  K peak 
is  caused  by  some  disorder  in  the  lattice,  perhapB  by  yttrium  atoms  moving  be- 
tween their  two  possible  siteB  within  the  dumbbell -shaped  cavity  under  the  ap- 
plied sonic  stress.  The  rise  in  attenuation  at  temperatures  above  200° K may 
be  the  intrinsic  Akhieser  loss.  If  so,  the  loss  is  much  higher  than  in  0-boron 
at  these  temperatures,  as  found  by  Young,  Oliver,  and  Slack  (Ref.  86);  however, 
the  temperature  dependence  is  similar. 

The  comparison  of  samples  of  YB—  with  YB,itV  does  not  show  any  appreci- 
able difference  in  acoustic  loss.  The  primary  difference  in  the  samples  is  prob- 
ably one  of  homogeneity  of  composition,  the  YE#.  , samples  being  better  because 
they  were  grown  at  the  congruent  melting  point.  However,  as  indicated  In  the 
discussion  of  the  crystal  chemistry  given  earlier  in  thlB  section,  the  is 

even  further  than  YBg,  from  the  stoichiometric  composition  of  YBM.  Neither 
is  very  close  to  the  possible  low -loss  composition  at  YB34  where  all  the  Y lat- 
tice sites  might  be  filled. 

More  work  is  needed  for  an  understanding  of  the  sources  of  acoustic  loss 
in  YBa.  Clearly  it  is  not  a good  material  for  delay  lines  as  it  is  now  made. 

For  any  immediate  applications  0-boron  or  spinel  would  be  better. 

K.  ELECTRICAL  PROPERTIES 

Electrical  properties  of  crystal  samples  can  be  measured  quickly,  and  they 
give  useful  information  with  regard  to  impurities  and  defects.  Data  for  YB* 
are  shown  in  Figure  43,  along  with  results  for  0-boron  from  Delta  and  Hermann 
(Ref.  67)  and  for  a -boron  from  F.  H.  Horn  (Ref.  60),  Four  alloy  contacts 
were  made  to  the  sample  with  platinum  wires  by  heating  the  sample  near  the 
Pt-YBM  eutectic  and  then  pulsing  current  through  the  contacts.  The  eutectic 
temperature  of  approximately  700° C limited  the  temperature  to  which  measure- 
ments could  be  made. 


The  temperature  dependence  of  the  conductivity  shows  the  material  to  be 
a semiconductor.  The  thermal  activation  energy  calculated  from  the  slope  of 
the  p versus  l/T  curve  is  0,88  eV.  This  value  does  not  correlate  well  with 
optical  transmission  which  begins  near  3000cm"1  or  0,36  eV.  The  resistivity 
data  and  the  optical  data  may  be  controlled  by  the  impurities  present  rather 
than  the  bandgap  of  the  material.  All  samples  were  found  to  be  p-type,  in- 
dependent of  stoichiometry  when  tested  for  polarity  of  Seebeck  EMF.  Four 
point  probe  measurements  of  sample  resistivity  showed  nearly  identical  re- 
sistivity at  room  temperature  for  all  samples.  Attempts  were  made  to  measure 
Hall  mobility  with  six  contact  samples  between  300°  K and  100° K,  but  ITall  volt  - 
ages were  undetectable.  All  of  these  results  suggest  impurity-dominated  con- 
duction and  the  desirability  of  preparing  samples  with  impurity  concentrations 
reduced  by  a factor  of  ten  or  more. 
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Figure  43.  Electrical  Resistivity  Versus  Reciprocal  Temperature 
for  YR6a#s,  a-boron,  and  /3-boron 
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The  purity  of  YB—  crystals  has  been  described  previously  (Refs,  43  and  45), 
and  the  yttrium  metal  used  was  shown  to  be  the  principal  source  of  Impurities 
(see  Table  10),  Sample  purity  nan  be  improved  by  a combination  of  zone  refihing 
and  use  of  higher -purity  yttrium  metal.  The  zone  refining  methods  which  are 
desirable  for  obtaining  uniform  samples  of  the  congruent  melting  composition 
will  also  remove  Impurities  with  segregation  coefficients  less  than  unity.  Seg  - 
regation coefficients  for  impurities  in  YB—  are  unknown;  approximate  values 
will  be  a byproduct  of  the  crystal  growth  by  zone  refining. 

The  purity  of  th<  yttrium  metal  used  is  very  important.  Purification  of 
rare  earth  metals  has  been  reviewed  by  Hufflne  and  Williams  (Ref.  69).  Zone 
refining  of  yttrium  has  little  effect  on  a number  of  the  most  prevalent  impurities 
(Ref.  45),  including  oxygen,  nitrogen,  copper,  nickel,  iron,  and  titanium.  Ex- 
periments in  zone  refining  of  YB,,  will  have  to  be  conducted  in  order  to  learn 
whether  any  of  these  elements  have  favorable  segregation  coefficients  in  YBW. 
Three  techniques  are  applicable  to  the  purification  of  yttrium  metal;  distillation 
under  high  vacuum,  filtering  through  porous  tungsten,  and  electrolysis.  The 
most  promising  is  distillation.  Initial  attempts  to  procure  distilled  yttrium 
fromthemetal  suppliers  proved  disappointing;  it  was  less  pure  (Ref.  45)  than 
vacuum-melted  metal.  N.H.  Parsons  (Ref.  70)  haB  reported  quite  good  re- 
sults in  distillation  of  erbium  for  the  Atomic  Energy  Commission.  Arrange- 
ments were  made  through  the  AEC  to  have  a quantity  of  yttrium  purified  by 
this  process  and  to  have  an  Impurity  analysis  provided  of  both  the  starting  ma- 
terial and  the  distillation  product  (see  Table  10). 

Purification  by  filtering  through  porous  tungsten  removes  some  impurities 
and  contributes  others.  It  could  not  be  intelligently  applied  to  this  problem  un- 
til the  effects  of  the  other  purification  procedures  are  known  for  a number  of 
impurities.  Purification  by  electrolysis  is  a potentially  valuable  purification 
method  and  could  be  applied  following  distillation,  provided  the  elemental  anal- 
yses of  other  purification  procedures  show  that  it  is  warranted.  As  in  the  other 
procedures,  the  impurity  segregation  achieved  (Refs.  69  and  71)  is  different 
for  each  impurity. 
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Section  5 


IRON-DOPED  SPINEL  CRYSTALS 

The  possibility  of  using  Fe8+  for  doping  of  spinel  crystals  in  order  to  re- 
duce the  acoustic  attenuation  at  300s  K has  been  explained  in  a previous  report 
(Ref.  43).  One  of  the  previous  limitations  on  this  possibility  was  the  fact  that 
only  natural  spinel  crystals  containing  iron  were  available,  and  these  crystals 
contained  other,  unwanted  impurities  as  well.  With  the  successful  growth  of 
spinel  by  flux  techniques  (Ref.  72)  it  has  been  possible  to  produce  synthetic 
crystals  containing  only  iron  as  an  impurity.  A number  of  these  crystals  were 
made  available  to  us  from  RCA,  through  the  Wright -Patterson  Air  Force  Base 
and  Mr.  Robert  Hickmott.  Three  crystals  were  selected  for  optical  study,  as 
listed  in  Table  14. 


Table  14 

SYNTHETIC  IRON -DOPED  SPINET,  CRYSTALS 


Crystal 

Label 

Weight 

Percent 

Iron 

Sample 

Thickness 

(mm) 

. cm-1 

Fep+/  Fct-u 

RCA  1 

0.15: 

5.71 

4.7 

0.28 

Novak 

0.40  * 

3.15 

8.1 

0.18 

RCA  5 

1.34: 

2.34 

37.0 

0.25 

*Iron  concentration  as  stated  on  crystals  as  received 
•^'•lron  concentration  aa  measured  by  standard  chemical  methods 


The  optical  absorption  spectra  of  these  crystals  were  measured  over  the 
wavenumber  range  500  cm-1  <v<  20,000cm"1  at  300°K.  The  broad  absorption 
band  characteristic  (Refs.  17  and  18)  of  Fea+in  tetrahedral  sites  was  clearly 
seen.  This  band  peaks  at  v = 5130  cm"1  and  extends  from  about  v - 3000  cm"1 
to  v = 6000  cm"1.  The  optical  absorption  coefficient  at  the  peak  is  given  in 
column  4 of  Table  14.  From  this  value  and  the  assumption  that  all  of  the  iron 
in  natural  spinel  sample  R97  of  Reference  72  occurs  as  Fe8"1- , one  can  cal- 
culate the  absolute  Fes+  concentration  in  tetrahedral  sites  in  the  synthetic 
spinels.  This  concentration  turns  out  to  be  18  to  28  percent  of  the  total  amount 
of  Iron  in  the  sample  as  measured  chemically. 

The  remainder  of  the  iron  may  be  present  as  metallic  inclusions,  Fe3+ 
in  octahedral  sites,  or  as  Fea  + . In  view  of  the  lack  of  control  of  the  partial 
pressure  of  oxygen  over  the  flux  during  growth  (ReT.  42),  it  is  not  surprising 
that  not  all  of  the  iron  is  Fe34-  ions  that  should  be  effective  in  reducing  the 
microwave  acoustic  losses  at  300°  K,  more  control  during  growth  is  needed  to 
assure  this  condition.  In  view  of  this  problem  with  the  iron  content  and  the 
small  size  (0.6  cm  maximum)  of  the  crystals,  no  acoustic  attenuation  measure- 
ments were  made  on  these  samples. 
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Section  6 


SOUND  ATTENUATION  AND  DIFFRACTION  LOSS 


Sound  attenuation  measurements  have  been  conducted  on  samples  of  YBm, 
MgAla04,  and  A1803  from  4®K  to  300°  K.  The  thermal  phonon  losses  have  been 
evaluated  previously  (Ref.  43)  for  MgAla04,  YAG,  and  AlaOa  at  frequencies  of 
500  MHz  and  above.  Similar  measurements  on  YB^  at  500  MHz  were  unsuc- 
cessful except  at  very  low  temperatures  (near  4°K).  The  present  attenuation 
measurements  were  conducted  at  low  frequency,  60  MHz,  using  bonded  quartz 
transducers.  This  report  gives  experimental  values  for  diffraction  losses  in 
all  of  these  materials,  so  that  one  is  no  longer  forced  to  rely  on  calculation 
(see  Section  2). 

The  attenuation  at  60  MHz  for  temperatures  below  about  100°  K is  dominated 
by  diffraction  losses.  At  higher  temperatures,  thermal  phonon  losses  (Ref.  43) 
may  enter.  For  example,  in  Figure  38,  for  YBj*  using  longitudinal  waves  the 
attenuation  below  100° K is  constant  at  0.79  dB/^is,  which  is  close  to  the  cal- 
culated loss  of  0.75  dB//is.  The  very  close  agreement  may  be  fortuitous  in 
view  of  the  uncertainty  in  b for  YBtf  (see  Equation  64).  The  results  of  the  com- 
pleted measurements  for  MgAla04,  Ala03,  and  YB^  are  given  in  Table  15. 

The  agreement  between  theory  and  experiment  is  good. 

Table  15 

SOUND  ATTENUATION  AND  DIFFRACTION 
LOSSES  AT  60  MHz 


Sample 

Mode 

Velocity 
xlO6  cm/sec 

Anisotropy 

b 

Measured 

Attenuation 

(dB/ps) 

Calculated 
Diffraction 
Loss  (dB/ps) 

MgAl804 

Longitudinal 

r1001 

8.56 

-1.1 

1.7 

1.5 

MgAla04 
RRE  2 

Shear 

noon 

6.63 

+ 1.5 

0.15 

0.17 

ai8o3 

A13213 

Longitudinal 

c-axis 

11.1 

+ 1.5 

0.15 

0.31 

Longitudinal 

noo] 

12.2 

+ 0.05 

0.79 

0.75 
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Section  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  and  related  recommendations  resulted  from  the 
studies  reported  heres 

1.  Magnesium  aluminate  spinel  is  shown  to  be  the  best  commercially 
available  material  for  acoustic  delay  lines,  in  both  attenuation  and 
beam  steering  characteristics.  Substantial  improvements  in  delay 
line  performance  can  be  achieved  with  spinel  and  appropriate  trans- 
ducers. 

2.  Appropriate  transducers  can  be  made  with  bonded  quartz,  or  LlNbQ,, 
or  deposited  ZnO  films.  In  all  cases  new  electrode  metals  must  be 
used.  Design  criteria  reported  include  transducer  characteristics, 
attenuation,  diffraction,  and  beam  steering. 

3.  A method  has  been  developed  to  calculate  the  effect  of  acoustic  im- 
pedance mismatch  between  an  acoustic  delay  line  and  its  transducers 
on  the  performance  of  a matched  filter  processing  a stored,  coded 
signal. 

4.  For  long  time  memory  of  coded  signals,  an  impedance  mismatch 
greater  than  a few  percent  can  lead  to  substantial  distortion  of  the 
output  from  a matched  filter.  With  the  allowable  distortion  speci- 
fied, the  required  degree  of  impedance  matching  can  be  calculated. 

5.  The  metal  system  Au-Pd-Ag  is  identified  for  bonding  transducers  to 
delay  lines.  Its  advantages  include:  the  ability  to  match  impedances, 
and  the  reduction  of  attenuation  due  to  both  electrons  and  dislocations. 

6.  The  efforts  on  spinel  directed  toward  making  a shear  wave  transducer 
bond  to  single  crystals  at  300  MHz  were  only  partially  successful. 

No  fundamental  limitations  were  encountered,  but  much  more  work 
and  attention  to  detail  will  be  needed  to  perfect  the  metal -film  trans- 
ducer bonding  technique. 

7.  The  work  on  YB#  demonstrated  that  the  anomalous  acoustic  attenuation 
in  either  YB#U,,  or  YB^,  is  quite  large  at  60  MHz  and  above.  The  rea- 
sons for  this  loss  are  not  clear,  but  one  of  them  may  be  the  presence 
of  the  dumbbell -shaped  cavities  that  are  present  at  the  lattice  sites  of 
the  yttrium  atoms.  These  might  be  filled  by  using  crystals  such  as 
LuB34. 


Another  solution  would  be  the  use  of  pure  /1-rhombohedral  boron, 
where  no  interstitial  metal  atoms  are  needed  to  stabilize  the  structure. 
Another  possibility  might  be  to  use  crystals  of  tetragonal  boron  which 
are  stabilized  by  carbon  or  nitrogen  atoms.  Clearly  a lot  of  work 
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would  still  be  involved  to  find  a material  better  than  spinel  for  room- 
temperature  microwave  acoustic  delay  lines. 

Sound  attenuation  measurements  at  low  temperatures  and  lou r fre- 
quencies have  confirmed  the  previously  calculated  values  of  diffrac 
tion  loss  in  MgAlaOv  A1203>  *nd  YB*  crystals. 
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